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Die stationiren Bewegungen des selbsterregten Kreisels 
und ihre Stabilitat. 


Von R. Grammel. 


1, Einleitung. Als selbsterregt soll ein Kreisel (d.h. ein Kérper, der um einen in ihm und 


im Raum festen Punkt O unbeschrankt frei drehbar ist) bezeichnet werden, wenn der Moment- 


vektor J{ seines Antriebs ein im Kreisel fester Vektor ist — vorbehiltlich einer spateren 
Erweiterung dieser Definition —, etwa verwirklicht durch einen oder mehrere korperfeste 
RiickstoBantriebe (wie sie bei einigen Anwendungen vorkommen). Seine Bewegung gehorcht 
den bekannten Eulerschen Gleichungen, die in einem rechtshandigen, kérperfesten Haupt- 
achsenkreuz (x, y, z) mit den Haupttragheitsmomenten A, B, C, den Komponenten @,, Wy, @; 
des Drehvektors 0 und den in diesem ‘Falle festen Komponenten M. M,, M, des Vektors 


x? 
lauten: : 


Air, +(C— B) 0,0, = Mu | 
Boy, 5 (A a C) O,0, = M,. (1) 
Co, (3. — A) 0,0, = M,. | 


Wir wollen hier aber nicht die Integration dieser Gleichungen versuchen, die meines Wissens 
bisher nur fiir den symmetrischen Kreisel (B = C) vollstandig gelungen ist!, sondern die Frage 
beantworten, ob ein derartiger Kreisel stationdre Bewegungen, und zwar zunachst solche in 
der Form von permanenten gleichférmigen Drehungen um sowohl kérper- wie 
raumfeste Achsen vollziehen kann und, wenn ja, ob diese stabil sind, — Bewegungen also, 
bei denen @, =, =@, =0 ist. Dabei wollen wir aber von vornherein die zu ft =0 gehérenden 
permanenten Drehungen um Hauptachsen aufer Betracht lassen, da deren Stabilitatsver- 
halten langst bekannt ist, und also JJt--0 voraussetzen. Wir werden finden, daB es dann beim 


| symmetrischen Kreisel eine (iibrigens auch schon bekannte, aber noch nicht auf ihre Stabilitat 
) hin nachgepriifte) co®-fache Schar von solchen permanenten Drehungen gibt, welche ,,fast 
| alle‘‘ stabil sind, beim unsymmetrischen Kreisel dagegen eine 0o*-fache Schar solcher perma- 
) nenten Drehungen, welche ,,fast alle«« labil sind. Weiterhin werden wir noch untersuchen, 
| wie man diese Labilitat beseitigen kann. Und schlieBlich werden wir beweisen, da es auber 
| diesen permanenten Drehungen keine allgemeineren stationaren Bewegungen des selbst- 
| erregten unsymmetrischen Kreisels gibt. 


2. Der symmetrische Kreisel. Es sei B= C; dann verlangen? die Gleichungen (1) mit 


| @, =O, =, =0 offensichtlich M, =0. Weil also der kérperfeste Drehmomentvektor Jt in 


der Aquatorebene des Kreisels durch 0 liegt, diirfen wir ohne Einschrankung M, =0 setzen, 
indem wir die y-Achse in die Richtung des Vektors St fallen lassen. Dann folgt aus (1), daB 


bei konstanten Drehkomponenten w, und w, notwendigerweise w, verschwinden muB: 
Wj, OK, wy =9, OSC, (2) 
und da8 die Konstanten a und c mit M; =|3)t| verkniipft sind durch die Beziehung 
M, =(A— B) ae. (3) 


Dies ist eine althekannte stationare Bewegung, namlich eine gleichférmige Drehung des 
Kreisels um eine Achse g durch O in einer Meridianebene, auf welcher der Drehmomentvektor Jt 


vom Betrag (3) senkrecht steht, mit der Drehgeschwindigkeit w = j@ + c? im einen oder im 


anderen Drehsinn um die Achse g (da nur das Vorzeichen des Produktes ac festliegt), also 


1 U. T. Bédewadt, Math. Z. 55 (1952) S. 310. ia 
2 Den belanglosen Fall 4 = B= C des sogenannten Kugelkreisels wollen wir hier und im Fol- 


) genden stets ausschlieBen. 


11 
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eine regulare Prazession ohne Eigendrehung. Die permanenten Drehungen um die Symmetrie 
achse (c =0) oder um die Aquatorachsen (a =0) zahlen dabei verabredungsgemaB nicht mit 
da fiir sie auch M, — 0 wird. | 

Um die Stabilitit dieser stationaren Bewegungen unter der Wirkung des Momentes (3) 
festzustellen, fiigen wir dem Drehvektor 0(a, 0, c) (2) eine Kleine Stérung hinzu, setzen alsog 
statt (2) 

@,=a+t ee, wy= ee, w,= eee’, (4 
wobei €1, 9, €, als klein von erster Ordnung gelten sollen, und prifen nach, ob die Stérungy 
unter der Wirkung des Drehmomentes (3) anwachst oder nicht. Fihrt man die Werte (4) 
in die Gleichungen (1) mit B=C, M, = M,=0 ein und vernachlassigt Produkte der kleinent 
GréBen ¢;, so kommt wegen (3) 


Ave, =e | | 
Boe, + (A— B)(ce, + aes) =0, (5) 
Boe, — (A— B)ae,=0. | 


Damit ¢,, €,, €,; von Null verschieden sind (also eine Stérung vorhanden ist), mufs die Deter 
minante der Koeffizienten von ¢é, &, €; der drei Gleichungen (5) verschwinden: | 


Ao 0 4 0 
(A —.B)ec Bo (A— B)a|=0, 
0 —(A— B)a Bo 
und dies liefert die drei Werte 
A—B Z 
a=9, @ .,=+ j wetaths (6)) 


Mithin gibt cs fir a + 0 keine Stérung mit positivem Realteil von 0, also keine Stérung, diet 
vom Drehmoment (3) angefacht werden kénnte; vielmehr gilt: 
Die permanenten Drehungen des selbsterregten symmetrischen Kreiselss 
sind, sowie ihre Komponente a+ 0 ist, alle stabil, oder besser ausgedrickt: 
parastabil! (weil die Wurzel 9 = 0 anzeigt, daB die Stabilitat nicht vollstandig in dem 
-Sinne ist, daB jede Stérung in Schwingungen um die urspriingliche Bewegung bestiinde). 
Von der permanenten Drehung um die Symmetrieachse (c = 0) ist bekannt, daB sie voll- 
standig stabil ist, von den permanenten Drehungen um alle Aquatorachsen (a = 0), daB sie# 
labil sind (was auch daraus hervorgeht, da®B dann nach (6) 0, = 9, = 0 wiirde, und das hat 
zusammen mit 0, =—0 bekanntlich Lésungen von der Form ¢ + ¢’t + ¢’’t? zur Folge). 
Man kann sagen, daf die Parastabilitat der anderen Achsen vom ,,Pol‘* nach dem ,,Aquator® 
hin von vollstandiger Stabilitat bis zur Labilitat abnimmt und mit Annaherung an den Aquator 
praktisch nur noch einen sehr geringen Stabilitatsgrad besitzt. | 


3. Der unsymmetrische Kreisel. Wir diirfen von jetzt an ohne Einschrankung / 
ie BC (7 


voraussetzen. Aus (1) folgt dann mit den Festwerten 


(Og Oe Oe (8) 

zunachst | 
Leip oe eae 

(A— €)eca = M,, O 

(A — B)ab = — | 


Aus (9) schlieBt man, daB keine der drei Komponenten M,, M,, M, des Drehmomentes Jt} 
fiir sich allein verschwinden kann: 
Drehmomentvektoren, die in einer der drei Hauptebenen von O liegen) 


(aber in keiner Hauptachse), vermigen keine permanenten Drehungen zu. 
unterhalten. | 


1 Fiir die Definition der Parastabilitat vgl. R. Grammel, Der Kreisel Bd. I, S$. 57, S. 214, S. 240, | 
2. Aufl. Berlin-Géttingen-Heidelberg 1950. 


j 


= 
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Dagegen gibt es nach (9) permanente Drehungen, wenn zwei von den drei Komponenten 

des Vektors St verschwinden: 
2 Drehmomentvektoren, die in einer der drei Hauptachsen liegen, vermigen 
eine permanente Drehung zu unterhalten, deren Drehvektor in der zu jener 
-Hauptachse senkrechten Hauptebene liegt. 
Zwischen den Drehkomponenten b,¢ bzw. c,a bzw. a,b und dem zugehérigen Dreh- 
“moment | Jt} = M, baw. = M, baw. = M, gilt die erste bzw. zweite baw. dritte Gleichung 

(9), und offenbar sind diese Drehungen die Verallgemeinerungen der permanenten Drehungen 
des symmetrischen Kreisels von Abschnitt 2. 

| Sind schlieBlich alle drei Komponenten des Vektors Jt von Null verschieden, so gibt (9) 
| aufgelést die folgenden Quadrate der Komponenten der permanenten Drehungen (8): 


1 Heo “MM, 


| Srp ied 6 My M; 
i (4— 8) (4 — C) Mz, 
jays A—C M:M. | 
a BR ec) i, (10) 
wees Vig? Mz My | 


Die Vorzeichen von a, b, c sind dabei durch (9) bestimmt, und dies liefert die Zuordnung der 
-nebenstehenden Tabelle, aus der man folgendes entnimmt: 


Zu jedem vorgegebenen permanenten | i SonteMe [Tal nib: ; | 
Drehvektor o(a,6,c) gibt es ein kérperfestes 
Drehmoment Jt(M,, M,, M,), das diese Drehung 
zu unterhalten vermag; aber nicht jedes vor- 
gegebene Drehmoment Jt kann eine reelle 
permanente Drehung 0 unterhalten, sondern 
nur solche Drehmomente (von dannallerdings 
beliebiger GréBe), deren Vektoren Jt in einem 
der vier Oktanten des Hauptachsenkreuzes 
liegen, in welchen die gerichteten Achsen 
des gréBten, des mittleren und des kleinsten 
Haupttragheitsmomentes in dieser Reihen- 
folge ein rechtshandiges System bilden, und 
der zugehérige Drehvektor 0 liegt dann jedes- 
mal in demjenigen Oktanten, der aus dem 
Oktanten von Yt durch eine Viertelsdrehung 
um die zugehérige gerichtete Achse des 
gro8ten Haupttraigheitsmomentes A im Sinne einer Rechtsschraube hervor- 
geht, oder auch in dem beziiglich des Nullpunktes O dazu symmetrischen 
Oktanten. (Zu jedem erlaubten Drehmoment 3) gehéren also auch beim unsymmetrischen 
Kreisel zwei permanente Drehungen von entgegengesetztem Drehsinn, aber gleichem Betrag.) 

Multipliziert man die drei Gleichungen (9) der Reihe nach mit den Faktoren —a, --b, —e, 
so kommt aM, + bM, + ¢M, = 0 und dies besagt: 

Der kérperfeste Momentvektor J steht senkrecht auf dem Vektor 0 der 
von ihm unterhaltenen permanenten Drehung. 

Weil das Drehmoment Jt demgem48 bei der Bewegung keine Leistung vollbringt, so geht 
hieraus auch noch hervor, daB dieses kérperfeste Drehmoment zwar die permanente Drehung 0 
‘unterhalten, aber nicht hervorbringen kann. Um sie einzuleiten, bedarf es eines anfanglichen 
| DrehstoBes mit den Komponenten Aa, Bb, Cc. 


| 

| 

| 
hax 


| 
4 


| 
| 
l 
| 


| 
| 
| 
I 


| 
I 
| 
I 
| 


ee ee Se 
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4, Die Stabilitait der permanenten Drehungen. Wir untersuchen auch hier die Stabilitat 
idieser Drehungen, indem wir zu a, b, c eine kleine Stérung hinzufiigen, also mit den Ansatzen 


| O, =a +e, e*, @y, =b + ee", O, = c¢ +6é,ee* (11) 


‘in die Gleichungen (1) eingehen, dann Produkte der ¢; vernachlassigen und auf (9) achten. 
So entstehen drei lineare homogene Gleichungen fiir die ¢;, deren Koeffizientendeterminante 
die folgende Form hat und verschwinden muB: 

ie 
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Ao —(B—C)c —(B—C)b 
(A — C)e Bo (A— C)a|=0 (12 
— (A— B)b —(A— B)a Co 
oder geordnet | 
oo? +a,0? +a,e +a, = 9 (13) 
mit den Beiwerten | 
Op == 24 DG | 
| = 0, asi 


= A(A— B) (A= Oe — B(A—B) (BC) P + C(A = 0) (B= Oe 
poe 2) (AB) (A Cy 6B nt oe 


Wir nehmen zuerst den Sonderfall a = 0 (Drehachse in der (y,z)-Hauptebene). Hier hay 
die Gleichung (13) auBer einer Wurzel 0, = 0 zwei entgegengesetzt gleiche 

reelle c ~ B(A — B) 

mare | Warzeln 0o,3, falls BS CASO)’ (15) 

im ersten Falle also eine reelle positive Wurzel 9,, die zu Labilitat der Bewegung fihrt, und 

dies trifft auch noch auf den Grenzfall des Gleichheitszeichens in (15) zu. Im Falle imaginarey 
Wurzeln Q,,, wird man wegen der Wurzel 0, = 0 wieder nur von Parastabilitat sprechen. 

In gleicher Weise kommen fiir b = 0 (Drehachse in der (x, z)-Hauptebene) auBer 0, = 

in jedem Falle zwei imaginaére Wurzeln 0,,,, was Parastabilitat bedeutet. 

Endlich fiir ¢ = 0 (Drehachse in der (x, y)-Hauptebene) hat man aufer 9, = 0 zwei ent- 

gegengesetzt gleiche 


reelle \ 


imagindre j 


2 

Wurzeln 0,,,, falls e $ Tao ' (16) 
im ersten Falle also Labilitat, einschlieBlich des Gleichheitszeichens in (16), im zweiten Fallet 
Parastabilitat. Wir fassen zusammen: 

Abgesehen von den villig stabilen permanenten Drehungen um die Achse} 
des gréBten und um die Achse des kleinsten Tragheitsmomentes (Jt =O 
sind alle permanenten Drehungen des selbsterregten unsymmetrischen Krei- 
sels um die Nullpunktsachsen der Hauptebene des gréften und kleinsten| 
Tragheitsmomentes parastabil, ebenso um diejenigen Nullpunktsachsen der 
Hauptebene des kleinsten und mittleren Tragheitsmomentes, die mit de» 
mittleren Hauptachse einen Winkel : 


p> = arctg /24—#) (17)| 

bilden, und ebenso um diejenigen Nullpunktsachsen der Hauptebene des; 
gréBten und mittleren Tragheitsmomentes, die mit der mittleren Hauptachse: 
einen Winkel : 
P > Pp = arc tg joao (18) 

bilden (in den letzten beiden Fallen also hinreichend weit von der mittleren Hauptachse 
abstehen, zu welcher bekanntlich eine labile permanente Drehung mit Jt = 0 gehdrt). Die 
permanenten Drehungen um die itbrigen Nullpunktsachsen in den Haupt- 
ebenen des kleinsten und mittleren sowie des gré8ten und mittleren Trag-| 
heitsmomentes (y <q, und ym <4q,) sind labil. 

Und nun werden wir vollends beweisen: 

Die permanenten Drehungen um alle anderen Nullpunktsachsen des selbst- 
erregten unsymmetrischen Kreisels, die in keiner der drei Hauptebenen. 
liegen, sind labil. | 

Dies folgt in der Tat rasch aus (13) und (14) an Hand der Hermite-Hurwitzschen Bedin- 
gungen? dafiir, dafs eine Gleichung dritten Grades keine Wurzeln mit positivem Realteil be: 


* Sie werden in der Regel nach Hurwitz oder nach Routh benannt, finden sich aber im wesentlichen 
cee bei Ch, Hermite, Crelles Journ. 52 (1850) S. 39, auf welchen sich A. Hurwitz ausdriicklich be- 
zient, 
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sitzt. Sie lauten 


Ope OE Oe 0G a0 Tapas a, ag 0% (19) 
In unserem Falle ist a > 0, a, =0 undo, +0 fir a+0,b = 0, c +0, und somit wider- 
sprechen sich die dritte und vierte Bedingung (19) fiir a +0, b +0,c +0, womit der Be- 
weis erbracht ist. 

_ Das Ergebnis, daB die von einem kérperfesten Drehmoment Jt unterhaltenen perma- 
nenten Drehungen eines unsymmetrischen Kreisels ,,fast alle‘ labil sind, ist erstaunlich, zumal 
da es einem anderen bekannten Ergebnis der Kreiseltheorie zunichst schroff zu widersprechen 
‘scheint. Wie Staude entdeckt hat, besitzt nimlich ein unsymmetrischer Kreisel, dessen Schwer- 
punkt nicht mit seinem Stiitzpunkt O zusammenfallt, unendlich viele Achsen durch 0, um 
welche er, wenn sie lotrecht gestellt sind, unter dem Einflu® der Schwerkraft permanente 
‘Drehungen vollziehen kann. Diese Achsen liegen auf einem Kegel zweiter Ordnung, der durch 
die drei Hauptachsen geht und also stets auch unendlich viele Achsen enthalt, die in keiner 
der drei Hauptachsen liegen. Die Stabilitatsverhaltnisse der Staudeschen Drehachsen sind 
genau bekannt!: der Staudesche Kegel enthalt im allgemeinen neben unendlich vielen labilen 
auch unendlich viele stabile Drehachsen, also auch solche, die wir soeben als labil nachge- 
wiesen haben. Nun werden aber doch auch die Staudeschen permanenten Drehungen durch 
ein kérperfestes Drehmoment unterhalten, namlich durch das Moment Mt, des Kreiselge- 
wichts beziiglich O, und dieses Moment St, kénnte man offenbar anstatt durch die Schwer- 
kraft ebensogut durch einen kérperfesten RiickstoBantrieb vom Moment Jt = Mt, erzeugen. 
‘Der scheinbare Widerspruch, der darin besteht, da®B dann eine bisher (unter der Wirkung 
‘von ‘)t,) stabile Drehachse plétzlich labil wird, lést sich indessen ganz einfach durch die Er- 
kenntnis auf, da®B das Moment Me der Schwere zwar bei der permanenten Staudeschen Drehung 
kérperfest ist, bei einer Stérung aber — infolge der Anderung der Schwerpunktslage gegen- 
iiber dem Stiitzpunkt — sich im Kérper (und im Raum) andert, wogegen ein seiner Natur 
nach kérperfestes Drehmoment, wie wir es hier immer voraussetzen, auch bei einer Stérung 
der permanenten Drehung im Kreisel fest bleibt. Die Stabilitat der wirklich stabilen Staude- 
schen Drehachsen wird eben dadurch erzeugt, da das Schweremoment Jt, bei einer Stérung 


‘sich selbst im stabilisierenden Sinne dndert. 
i 
5. Die Stabilisierung der labilen permanenten Drehungen. Die letzte Bemerkung fiihrt nun 


auch vollends zur Beantwortung der praktisch wichtigen Frage, mit welchen Mitteln man die 
als labil erkannten permanenten Drehungen eines selbsterregten unsymmetrischen Kreisels 
‘stabilisieren kann. Man muB offenbar das Drehmoment 3Jt in einem stabilisierenden Sinne 
‘steuern. Soweit dies vom Kreisel selbst aus méglich ist, wollen wir ihn immer noch selbst- 
erregt nennen, obwohl bei einer solchen Steuerung das Drehmoment St nicht mehr ein im 
kérperfesten Bezugssystem konstanter Vektor sein wird. 

Da man den Drehvektor 0 eines Kérpers grundsatzlich durch Hilfsmittel im Kérper selbst 
‘registrieren kann (die mechanischen Versuche zum Nachweis der Erddrehung leisten dies mehr 
‘oder weniger vollkommen) und dann auch in seine Komponenten @,, wy, @, im kérperfesten 
/System zu zerlegen vermag, so kommen vor allem gesteuerte Drehmomente 9{ in Frage, 
: deren Komponenten M,, M,, M, von @,, Wy, Wz abhangen. Der einfachste Ansatz, der, wie 
wir best 4tigen werden, auch zum Ziele fihrt, lautet 
[ 

} 

| 


M,, = My — M, Wy » 
Vee Ty Oy, 
et (ea eee 
Dabei sollen M;, My, M;, m,, m,, m, Festwerte sein. Die Komponenten a, b,c des zu- 
_gehdrigen permanenten Drehvektors 0 gehorchen nach (1) und (20) den Gleichungen 

M2 =m,a—(B—C)bc, | 

My =m, b--(A—C) ¢ @, 

M? =m,c—(A— B)ab. | 
‘Diese kann man zwar nicht in geschlossener Form nach a, b, c auflésen; man kann aber aus 
ihnen zu jedem Wertesatz a, b,c bei frei gewahlten ,, Steuerungskonstanten m,, m,, m, die 


,,Hauptkomponenten* Mz, My, M; von Nt berechnen. 


(20) 


(21) 


1 R. Grammel, Math. Z. 6 (1920) S. 124, oder Der Kreisel, Bd. I, S. 178, 
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Geht man jetzt mit dem Stérungsansatz (11) in (1) samt (20) ein und beachtet (21), sq 
kommt statt (12) 


Ao +m, —(B—C)c —(B—C)b ! 
(A— C)e Bo +m, (A— C)a|=0 (22) 
= Ap —(A— B)e Co+m, 
oder aufgelést, die Gleichung (13) mit den neuen Beiwerten 


aera AGBIGs | 
% = BCm,+CAm, +A Bm,, 
a = A(A— B)(A—.C) & — B(A— B)(B—C)P4C(A—C)(B—O)e2? + 
+Am,m, + Bm,m, + Cm, m, , 
& = 2(A— B)(A— C)(B— C)abe +m,m,m, + 
+(A— B)(A— C) @ m, — (A — B)(B— C) Bm, +(A— C)(B—C)em,. | 
Weil nach den Hermite-Hurwitzschen Stabilitatsbedingungen (19) neben o> 0 insbe- 


sondere auch a, und a, positiv sein miissen, so wird man die Steuerungskonstanten m, und m, 
positiv und m, zunichst einmal gleich Null wahlen. Mit m, = 0 findet man aus (23) vollendss 


Oy Oe — Oy Wg = B(A—.C)[(A— B) Aa? m, —2(A— B)(B— C)ACabe + | (24) 
+ (B— C) C? c2 m,] + B? [m, m,— (A— B)(B— C)b?](Cm,+Am,) (m,=0). | 
Jetzt geniigt es offenbar zur Stabilitat fiir alle Drehungen, bei denen b = 0 ist und min-! 


destens eine der Drehkomponenten a und c ebenfalls von Null verschieden ist, wenn man den} 
Steuerungskonstanten die Werte 


e| 


ie (AB) bi, my, = 0, m, = (B— C) |b! (25) 
erteilt; denn dann erhalt man nach (23) und (24) jeweils fiir bz0 
Og =A BC> OF 
OB A =) bil Ss 00y on 
a, = (A— C) [(A— B) a+ (B— C)c]*|b|> 0, (26 
Oy Bg — %) 0, = B(A— B)(A— C)(B—C)(AaF Ce)?|b|> 0. 


Somit sind alle Stabilitatsbedingungen (19) erfiillt und demgem4f{ alle Drehachsen, fiir) 
welche a+ 0, 6+ 0, c+ 0 ist, aber auch diejenigen, fiir welche eine der beiden Drek-} 
komponenten @ und c verschwindet, stabil geworden, also auch die bisher nur parastabilen 
Drehachsen in den drei Hauptebenen. 

Weil die Stabilitat der Hauptachse des gréften Tragheitsmoments (b = 0, c = 0) und} 
derjenigen des kleinsten Tragheitsmoments (a =0, 6 =0) von. vornherein feststeht, se 
bleibt jetzt nur noch die Stabilisierung der (bekanntlich von sich aus labilen) Hauptachse des 
mittleren Tragheitsmoments (a =0, b+ 0, ¢ = 0) ibrig, fiir welche natiirlich von vorn-} 
herein My = 0, M; = 0 und My, =m, b sein muB, da es sich um eine permanente Drelaaa 
mit Jt = 0 handelt. Hier folgt aus (23) mit a=0,c¢ =0 

a; = [m, m, — (A — B)(B— C) b?|m,, | 

oy Og — ay a, = B® [m, m,— (A — B)(B—C)¥](Cm, + Am,) + 

+2ABCm,m,m, + A Cm}(Cm, + A m,) + Bm, (C? m2 + A? m?) , | 


z 


‘— 


und man erkennt, daB mit 


Tit (As) Be m,> 0, m, = (B— C)|8| (28) 
alle Stabilitatsbedingungen (19) befriedigt werden, so daB nun auch die mittlere Hauptachse 
stabil geworden ist. 

Man stellt aber dariber hinaus leicht fest, da8 mit (28) auch in allen friiheren Fallen die 
Stabilitat ebenso verbiirgt ist wie mit (25), und hat somit das Ergebnis: | 

Mit Steuerungskonstanten m,, my, m,, welche die Vorschrift (28) erfillen, 
lassen sich alle permanenten Drehungen des selbsterregten unsymmetrischen 
Kreisels véllig stabil machen. 

Endlich wollen wir nur noch darauf hinweisen, da®B man die Stabilisierung statt mit dem 
Ansatz (20) in gleicher Weise auch mit dem viel allgemeineren Ansatz 


M, = F,(@,) , M, = F, (@,) , M, = F, (@,) (29) 
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} versuchen und untersuchen kann, worin F,, F,, F, gegebene Funktionen sind, die an der 
{ Stelle O, =a bzw. Wy = b bzw. w, =c Ableitungen besitzen. Hierbei ist a, b, c ein reelles Lé- 
sungstrippel der drei Gleichungen 


Pee ee 8 CAAT 
(4— C) ea — -+-F, (bh), (30) 
z (i 3R) ce F.40) il 
) Fir den Stérungsansatz (11) hat man jetzt die Werte 
Py (@z) = F, (a) + Fi (a) e, e*, | 
F, (wy) = F, (b) + F; (8) e, 8, : (31) 
FP; (,) = Fs (c) + F; (c) e, e*, 


fin (1) zusammen mit (29) und (30) einzufiihren. Man erkennt, da dann die bisherigen Steue- 
jrungskonstanten m,, m,, m, einfach durch die negativen Ableitungen — TG) — 1 (his 
>— F’, (c) zu ersetzen sind. Sonst dndert sich an unseren Rechnungen und Ergebnissen nichts. 


6. Beweis fiir die Nichtexistenz weiterer stationirer Bewegungen des selbsterregten unsym- 
# metrischen Kreisels. Wahrend der symmetrische Kreisel wenigstens in dem ausgearteten Fall 
) M =O (als kraftefreier Kreisel) auBer den in Abschnitt 2 erwahnten permanenten Drehungen 
s auch noch allgemeinere stationére Bewegungen vollziehen kann, namlich seine bekannten 
»regularen Prazessionen w, um eine raumfeste Achse nebst einer Eigendrehung w, um seine 
; Symmetrieachse, so gibt es bekanntlich beim unsymmetrischen kr4ftefreien Kreisel auBer 
‘seinen permanenten Drehungen um seine drei Hauptachsen keine weiteren stationdren Be- 
) wegungen. Hier erhebt sich nun die naheliegende Frage, ob nicht der selbsterregte unsymme- 
i trische Kreisel (J)(+-0) auBer den in Abschnitt 3 aufgezihlten permanenten Drehungen (um 
}raum- und zugleich kérperfeste Achsen) noch allgemeinerer stationérer Bewegungen fahig ist, 
)etwa einer Eigendrehung w,+0 von festem Betrag um seine (kérperfeste) Hauptachse und 
)zugleich einer Prazessionsdrehung w,-+0 von ebenfalls festem Betrag um eine raumfeste 
) Achse, die den festen Winkel 6+ 0 mit jener Hauptachse bildet. DaB eine solche Bewegung 
* beim selbsterregten unsymmetrischen Kreisel nicht méglich ist, lieBe sich leicht an Hand der 
| Gleichungen (1) zeigen. Wir wollen aber sogleich die viel allgemeinere Behauptung beweisen: 
~ Die permanenten Drehungen (10) sind die einzigen stationdren Bewegungen, 
/ deren der selbsterregte unsymmetrische Kreisel (bei festem MM) fahig ist. 

Die allgemeinste stationire Bewegung, welche ein starrer Kérper, der sich um einen 
‘festen Punkt O frei drehen kann, kinematisch zu vollziehen vermag, setzt sich zusammen 
aus einer Eigendrehung 0, von festem Betrag w, um eine kérperfeste Achse durch O, die aber 
» keine Hauptachse sein muB& (es jedoch sein darf), und einer Prazessionsdrehung 0, von festem 
Betrag w, um eine raumfeste Achse durch O, wobei die Vektoren 0, und 0, den festen Winkel 6 
_miteinander bilden. Legt man nun eine kérperfeste x*-Achse in die Richtung des Eigendreh- 
) vektors 0,, so ist die dortige Drehkomponente 


Ws = W, + W,cosd. (32) 
| Die kérperfesten y*- und z*-Achsen, die die «*-Achse zu einem rechtshandigen, rechtwinkligen 
Koordinatensystem erginzen, laufen in einer zum Vektor 9, senkrechten Ebene durch O mit 


| der Drehgeschwindigkeit w, um, und man darf dann mit einem geeignet definierten Nullpunkt 
) des Drehwinkels y der Drehung w, die zugehérigen Drehkomponenten in der Form 


: ot =o,sindsing, Wf=M)sindcosp mit Pp =O (33) 
| ansetzen. Durch (32) und (33) ist mit festen Werten ,, w,, 0 die allgemeinste stationare 
| Bewegung des Kreisels kinematisch dargestellt. 

Um sie auch dynamisch untersuchen zu kénnen, miissen wir die Eulerschen Gleichungen (1) 
» von dem Hauptachsenkreuz (x, y, z) auf unser jetziges Achsenkreuz (x*, y*, z*) transfor- 
»mieren. Dies geschieht am schnellsten folgendermafen. Die Tragheitsmomente um die x, 
\ y*-, 2*-Achse seien A*, B*, C*. Da sie im allgemeinen keine Haupttragheitsmomente sind, 
, so hat der Drehimpulsvektor D’, der zu einer Drehkomponente w* gehért, auBer einer Kom- 


 ponente Di, = A*@® in der x*-Achse im allgemeinen auch noch Komponenten in der y*- und 
z*-Achse, namlich 


le * * 
Dy = We Ons Dy = WA Ws > 


1 
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wenn man die sogenannten Kippmomente? 

We = — [y*o*dm, We = — [e* x* dm, We = — [x* y*¥ dm (34) 
des Kérpers einfithrt, wo jeweils tber die Massenelemente dm des Kreisels zu integrieren ist 
Somit hat man fiir die drei Komponenten D*, D*, Df des gesamten Drehimpulses im (x >| 
y*, 2*)-System | 

Ds =-A* wh Wi wy og Os (35) | 
und zwei daraus zyklisch hervorgehende Ausdriicke fiir D¥ und DZ. Zerlegt man jetzt die 
vektoriell geschriebene Eulersche Gleichung (1) 


dD 


re + [oD] = M 
in ihre Komponenten im (x*, y*, z*)-System, so kommt 
+ (C* — B*) ow} oF + Wolo, — Wore: = Mz 


sowie zwei daraus zyklisch hervorgehende weitere Gleichungen®. Hierbei sind natiirlich Mf, | 
M*, M# die festen Komponenten des Vektors {Jt im (x*, y*, 2*)-System. 
' Gehen wir jetzt mit (32) und (33) in die Gleichung (36) ein, so lautet sie 


(W* sin p — WX cos p) w, sind cos 6 + 2 (C*— B*) sin 29 — W cos 29] w, sin? d = MZ. 


Weil rechts das feste Moment M* steht, so muB auch die linke Seite einen festen Wert haben. | 
Hierfiir gibt es aber nur folgende Méglichkeiten: 
1. w, = 0, also gy = konst., 


2. w, + 0, dann aber entweder 
a) ®p = 0 (und also MZ = 0), oder 
b) 6 = 0 oder a (und also Mf = 0), oder 
c\Ou= 2. BC, We —0. (undsaleo si s—"0), oder 
d) WE. = WA = WS = 0,4) B* = C*) (und,alsot MZ .0)= 
Die Falle 1. sowie 2a) und 2b) sind offenbar nur wieder permanente Drehungen, wie wir sie 
von Abschnitt 3 her kennen. Die Falle 2c) bzw. 2d) lassen die erste der beiden aus (36) zy- 


klisch entstehenden Gleichungen tbergehen in 


A*o,0, cos » +-5 (Wt sin 2p + Wt cos 2g) 03 — W3 (o: =+ 0%) = M? (37) 


bzw. in 
[A*w, + (A* — B*)@, cos 6], sin 6 cos p = M3. (38) | 
Gleichung (37) verlangt entweder m, =0, also py =konst. oder w, =0 und bedeutet also eben- | 


falls eine permanente Drehung. Gleichung (38) verlangt entweder w,—0 oder 6 =0 oder Z, | 
was wieder permanente Drehungen bedeutet, oder aber das Erfiilltsein der Bedingung 


A@!--(A = 1B) aeos) = 0% (39) 


wobei wir die Sterne fortlassen durften, da es sich ja wegen 2d) um ein Hauptachsenkreuz | 
handelt. Man bestatigt, daB fiir diesen selben Fall 2d) die letzte der beiden aus (36) zyklisch 
entstehenden Gleichungen ebenfalls auf permanente Drehungen oder aber auf die Bedingung 
(39) fithrt, und diese Bedingung (39) ist ja einfach diejenige fiir die kraftefreie Prazession des 
symmetrischen Kreisels, welche hierdurch als die einzige mégliche stationare Bewegung auBer 
den permanenten Drehungen nachgewiesen ist. Damit ist der Beweis geschlossen. 


(Eingegangen am 1. Oktober 1952.) 
Anschrift des Verfassers: Prof. Dr. R. Grammel. Stuttgart-N, Robert-Bosch-Str. 101. 


+ Vel. etwa R. Grammel, Der Kreisel, Bd. I, S. 21. 
* Diese Gleichungen sind im deutschen Schrifttum wenig bekannt, aber nicht neu; schon L. Euler 


hat sie gefunden und veréffentlicht in seiner Arbeit: Découverte d’un nouveau principe de mécanique, 
Mém. Acad. Royale Berlin, 6 (1750) S. 213. 
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Anwendung der Matrizenrechnung zur Erfassung von Betriebsablaufen. 
Von O. Pichler. 


1, Einleitung. In einer fritheren Arbeit! wurde gezeigt, wie sich die Matrizenrechnung bei 
‘9 betriebswirtschaftlichen Aufgaben in einem aus vielen Teilbetrieben bestehenden Industrie- 
} werk mit Vorteil anwenden lat. Das angewandte Verfahren beruhte darauf, daB die Betriebs- 
§ vorgange eines einzelnen Teilbetriebes mit praktisch hinreichender Genauigkeit durch eine 
| Matrix beschrieben werden kénnen. Aus diesen einzelnen Teilmatrizen lat sich dann unter 
} Beriicksichtigung der zwischen den einzelnen Teilbetrieben bestehenden Kopplung und unter 
9 weitgehender Benutzung des Matrizenkalkiils eine das gesamte Betriebssystem kennzeichnende 
4 Strukturmatrix aufbauen. An zwei einfachen Beispielen, der unverzweigten Betriebskette 
fund der Betriebskette mit Riickfithrung wurde dies ausfihrlich erlautert. Im folgenden soll 
nun zunachst gezeigt werden, wie sich das Verfahren bei allgemeineren Kopplungsbedin- 
} gungen anwenden 1aft; sodann soll auf die weiteren Méglichkeiten eingegangen werden, die 
oder Matrizenkalkiil bietet, um die Matrizen in eine fiir betriebswirtschaftliche Aufgaben 
ii méglichst geeignete Form zu bringen. 


j 2. Konforme Kopplungsbedingungen. LEin Teilbetrieb i eines aus m Teilbetrieben be- 
» stehenden Betriebssystemes sei durch seinen als Vektor darzustellenden Durchsatz 0; gekenn- 
» zeichnet (s. Abb. 1). Der Durchsatz D; gibt die dem Teilbetrieb i zur Verarbeitung zuge- 
i fihrte Stoffmenge (man denke etwa an ein Gas mit n Komponenten) an, aus der neben der 
4 Ausbeute a;, die als Zwischenprodukt anzusprechen ist und anderen Teilbetrieben zur Weiter- 
verarbeitung zugeleitet wird, auch noch ein Fertigprodukt {; entsteht. Es werde angenommen, 
s daB sich der Zusammenhang zwischen 0;, a; und {; durch die Vektorengleichungen 

a; = U;d,, fi = Bid; (1) 
| mit praktisch hinreichender Genauigkeit beschreiben 148t. Um nun auch die Kopplung des 
) vorliegenden Betriebssystems zu beschreiben, mu® angegeben werden, in welcher Weise die 
| Ausbeuten zur Weiterverarbeitung auf die einzelnen Teilbetriebe verteilt % 

5 werden. In einem Ansatz, der eine solche Verteilung beschreibt, soll aber 
* zum Ausdruck kommen, da8 bei einer Aufteilung der Ausbeute a, in m Teile 

* gemaf3 


Betricb tv 


ese Onkee nk 
# samtliche Vektoren a,; dieselbe Zusammensetzung haben, so daB also a fi 


Abbildung 1. 
Qik = Wik 


/ mit skalarem w;; gesetzt werden mu8. Eine Kopplung dieser Art, die besonders in der che- 
4 mischen Industrie vorkommt, soll als ,,konforme Kopplung* bezeichnet werden zum Unter- 
 schied von der ,,durchsatzbedingten Kopplung“, bei der diese Zusatzbedingung nicht besteht 
) und die im Abschnitt 3d behandelt wird. Um zu einem endgiiltigen Ansatz fiir die konformen 
) Kopplungsbedingungen zu kommen, ist noch zu beriicksichtigen, da die Durchsatze D; nicht 
| allein durch die verschiedenen Q;,; gebildet werden, sondern da dazu auch noch von auSen 
an das Betriebssystem gelieferte Rohprodukte t; treten kénnen, hinsichtlich deren Verteilung 
s| dasselbe wie bei den Ausbeuten , gelten soll. Dann 148t sich eine Kopplung sehr allgemeiner 
#) Art durch den Ansatz 


m 1 m I 
D> Oy, DY = DS ie De, G (2) 
pea y= pa jet 


_beschreiben, der angibt, wie sich der Durchsatz 0; aus den Ausbeuten Q; der m Betriebe und 
/ aus den I in das betrachtete Betriebssystem eingehenden Rohprodukten zusammensetzt. Die 
) Verteilungskoeffizienten w;, und wj; miissen den Bedingungen 
| VS wy =, 0 suf, <1 
und 

i=l i=l 


| 10, Pichler, Ing.-Archiv 21 (1953) S$. 119, im folgenden zitiert mit if 
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geniigen. Die Ausdriicke: 
‘= (1 — a Wik) Ay und re =f (i! = wi) ty (3) 


geben die im Betriebssystem nicht weiter verarbeiteten und (mit Einschlu8 etwaiger Lei- 
tungsverluste) nach auBen abgegebenen Ausbeuten und Rohprodukte an und sollen als ,,Zwi- 
schenproduktabgaben“ bezeichnet werden. Dabei ist itbrigens auch der Fall, daB 3 vor- 
kommt, praktisch durchaus méglich. Die Wahl skalarer Verteilungsgréfen w;, und wr 5 ist, | 
wie bereits erwahnt, notwendig, um zum Ausdruck zu bringen, daB simtliche Komponenten } 
von a, und t; in gleicher Weise verteilt werden. Man kann nun ohne Einschrankung der All-| 
gemeinheit — wenn nétig durch Einfiigen weiterer Komponenten vom Betrag Null — an- 
nehmen, daf simtliche Vektoren D;, a, tj sowie die Fertigproduktvektoren {; von der Di-+ 
mension 7 sind; dann |aBt sich mit den Skalarmatrizen 


Wir = wir Cr und WF, = ui; C, (©, n-reihige Einheitsmatrix) 


die Betriebskopplung durch den Ansatz darstellen 


De % = (iz) | + (RBH) | tr | = (Wie) a + (WF) v, 
a W | 
| : (4) 
Repose ; 
oder mit Kinfiihrung der ,,Kopplungsmatrizen“ 
YW = (YW: ») (mn Zeilen und Spalten) , 
A= (207, (mn Zeilen und Jn Spalten): 
d= Wa + W* r. 
Fiihrt man noch die Matrizen . 
n=/(uN, 0 0 |} und S=(B 0 - O 
0 Uy, 0 Oinbouk ae” 
0 Ov ombihiils Own Oi aubss Bm | 
mit je mn Zeilen und Spalten ein, so ist das gesamte Betriebssystem durch folgende Gleichungen 
Ga Ds] : 
PSS (Umsatzgleichungen) (5) 
D= Wa + Wr (Kopplungsgleichungen) 


gekennzeichnet. Setzt man die erste Umsatzgleichung in die Kopplungsgleichung ein, so folgt | 
D=WUdD WF rs 
; (CF — WN) d = Wr. 
Ist (© — WW N) nicht singular, so wird 
D= (C— BU)! WF r, | 
a= U(C— WU) WF r, 
f= 8(EC— Bw Br 
oder zusammengefaBt | 
P= (Di) © | (C— BN) Wr Fr, 
a U 
haya) 
wo nunmehr der ,,Betriebsablauf* { als Funktion des Rohproduktvektors r dargestellt ist. 


Die Beschrankung, die in dem Ansatz (1) insofern liegt, als fiir jeden Teilbetrieb nur eine 
Ausbeute und nur ein Fertigprodukt vorgesehen ist, 148t sich durch den erweiterten Ansatz 


Qi ] =n, fia = %; 0D; 
Ais | fro (1a) 


a; sj | fi o; 


(6) 


(6a) 
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jaufheben. Man hat dann s; Ausbeuten und o; Fertigprodukte im Betrieb i, fiir das gesamte 


, Betriebssystem aber s = 5») s; Ausheuten und o = » 9; Fertigprodukte. In den Gleichungen 
‘tf i=1 ri 

(| (2) und (3) ist dann lediglich der obere Summenzeiger m durch s zu ersetzen, desgleichen in 
» Gleichung (4) a,, durch a, und die Matrix YW hat nicht mn, sondern sn Spalten. Sonst andert 
) sich an der ganzen Ableitung nichts. 


Eine Anwendung der Gleichungen (6) mége fiir die Kopplungsmatrizen 


} W=|( 0 » haa ; 0. An &, | und W* = E 
, AS, 90 . 0 0 (©, n-reihige Einheitsmatrix) 
Cee 0 oO 
re aed 
| ‘ok ails Po An pn PO 
© gezeigt werden. Bei dieser Wahl von ¥ folgt aus ) — Wa + YW* ¢ fiir A; = 1 
D: = a4 +7; 


ub (unter zyklischer Auffassung der Folge 1, 2,...m, also mit Q,_; = Qj) = Qn): 
Dieser Fall wurde in einer fritheren Abhandlung! als Betriebskette mit Riickfithrung be- 
‘weichnet. Fir J; = 1(i <m),/, =0, kommt 


Ds = a_, +1; 0, — the (l<tigm) 


ein Fall, der als unverzweigte Betriebskette? bezeichnet werden soll. Nun wird 


eee 0 0 LAG Oe 0, OE fo 0 0c 0M Ag tty) 

! PG 0 fac 0. 0 Ome srens 104) On| eel; 1,072) 0 0 
OSE a soon Ort 4G | OMS etaaty 6 Oily 0° 0 
SBC OAL E, 0 OIRO 90> tA Dig OAU AAO 9): 


| Zur Vereinfachung der Schreibweise werde voriibergehend 1; an Stelle von A; Ul; gesetzt, 
_ferner werden die Abkizungen 


pe Ui ely 
Bro; Ue | 
t..= (©, = ToS a 
BVin= Vr 


| unter zyklischer Auffassung der Folge 1, 2,...m eingefiihrt. 


Gh 182-8 tn) (7) 


Dann gilt die Gleichung: 


4 CS, 0 : 0 — ite | | Th DER vas 4 {De i < € 
z : | = i 7 (€ mn- 
1 ——- %, GC, ; 0 0 | Vn Vn Use : Uy m Un m | reihige 
,) : : 7 ‘ ‘ oF at , “_ |. Einheits- © 
0 0 - ine C,, | | ee. 1 Uns a5 Us x is m J matrix). 


(© — W UN) 
| Beim Nachweis dieser Gleichung ist lediglich zu beachten, daf 
Zeile i des Produktes = Zeile i des Rechtsfaktors — U,;_, - Zeile (i — 1) des Rechtsfaktors 


i 


ist. Somit gilt, wenn das Matrizenprodukt voriibergehend mit OS) = (Dix) bezeichnet wird, 
i fir i +k: Oy —= Wyk en a Cee Wax) = Wi_1,% Une — Wi. Wes = Ds 
® fir i — Kk: Osp—= Wer — Ves (Va_on Ven) = Ue — Uae Une = (©, — Ue) er 
= (© — Vn_1,%) (C, — Yea.) = &, 
| Damit ist gezeigt, da® der Rechtsfaktor in dem obigen Produkt die Kehrmatrix von (G — YW U) 
darstellt. Fiir diese 1a8t sich somit schreiben, wenn man fiir 1; nun wieder 4; U1; setzt und die 


1 [, Seite 125 ff. 
2 1, Seite 123 ff. 
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(zyklisch aufzufassenden) Abkiirzungen 
Mie As Asa. ++ Me» 
I= Tae e/g iiss (z, i 1, 2p kee m) (Tb) 


i=l 


Wir (A) — (C, —A Ven ye 


einfiihrt, 
(E— Buy =| 1:(A) Aina Uma Uoo(A) + Anim Ura m Urn m(A) | 
| Au Un My (A) U,,(A) ; Aim Miss Un m(A) 
es Ss ; ed ” 
= Veet Uy (A) Am—1,2 as U,,(A) ee Vem) ( ) 
ay Or, == Oe tt bate Ih 
= (ik C, —+ Amat Wa U, (A) ( 5 | : i i) 5 
Dabei wurde fiir die letzte Umformung die Identitat 
Ty(A) = ©, ae Ups Tan(A)= Cn + Anan Ven(A) Use (9) 


benutzt, die aus ©, =A Uy_4,. + (©, —A VU_1,x) folgt, wenn man beide Seiten von rechts 
bzw. von links mit (G, —A Ws.) = Uee(A) multipliziert und beriicksichtigt, daB gemaB 
(Tb) A = Ay_1,, ist. Setzt man noch Y* = & und beschrankt sich im Interesse der einfacheren 
Schreibweise auf m = 3, so erhalt man schlieBlich fiir den Betriebsablauf { aus (6a) mit den 
weiteren Abkiirzungen %;, gemaB (7) 


i d| a7 | U(A) Asx Uso Ug9(A) Asg Uss Tl, 9() : 
i Ay Uni Ui (A) U,(A) Arg Uys Us3(A) 
Asy Us Ui(A) Ase Uso Up9(A) U13,(A) 
| My, UW, (A) Ase Wyo Us9(A) Azz Uys Us3(A) 
Ayy Us) Wy (A) My5 Ug9(A) Ay3 Ugs Us0(A) (10) 

Asx Us, Uy (A) Age Uso W20(A) Us3 Ug(A) 

Bi Thy) don Bie Uaa(®) ag Bs Usa) | 

Ay Bar y(A) Boo Ug9(A) Ais Bos U53(A) | 

| As B31 Vi (A) Age B32 Ug2(A) B33 Ug3(A) J 


Unter Verwendung des Kroneckersymbols 6;; 1aBt sich allgemein fir eine Kette mit m Teil- 
betrieben schreiben 


{=/{bd)= Fe C, +A: ak Wee Uye(A) tie 
a Oin Ui + Aran Wie  Unn(A) (10a) 
f | Oink Bi t+Aan Bin UWser(A) 
Fiir Zwischenproduktabgaben folgt aus (3) unter Vertauschung der Zeiger i und k 
gi = (1 — A) a: (3a) 
Man erkennt ferner, daB man aus (10) bzw. (10a) fiir 2; = 1 die Betriebskette mit Riick- 
fithrung! erhalt (es wird dann 3; = (1—A,) a; = 0), wahrend mit 4; = 1 fiir i = 1, 2,... m —1, 
Am = 0 die Gleichungen fiir die unverzweigte Betriebskette? herauskommen. (Es wird dann 
si = 0 fir 1 = 1,2,...m—1 und 3, =a,,). Bei dem Nachweis ist zu beachten, daB in 
diesem Fall, da A,, = 0, sowohl 
AS da 4 ov be 0 A undedamity ti (€,, —A U,_1,4) = ©, , 
als auch fiir i << k 
Mine Within A A Ana. le und damit auch hik Ui. = 0 
werden. Allgemein gilt, daB, wenn r GréBen A; verschwinden, das Betriebssystem in r un- 
verzweigte Betriebsketten zerfallt. 


1 Siehe S. 128. 
2 Siehe S. 124, 
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Fiir die Aufstellung der Lésungsmatrix 9 sind m reziproke, n-reihige Matrizen Un (A) 
zu berechnen; diese lassen sich, wie jetzt gezeigt werden soll, samtlich auf eine von ihnen 
beispielsweise 11;;(A) zuriickfithren. Dazu wird die Identitat 


A(C, —A BA = (EG, —A AB) 
d bendtigt, die sich aus der Identitiat 
(©, —A AB) A = W(E, — a B A) 
* ergibt, indem man beide Seiten von links mit (&, — 4 2 8)! und von rechts mit (E, — 2B Y)-! 
* multipliziert. Nun ist 
Usain Wne(A) = Uae (C, —A Uae) = Wan (G —A Ue-se Us aae 
oder unter Benutzung der soeben gezeigten Identitat 
Use Vn.) = (©, — A Wa Ve.) Wei = (G, — A Wa) Wes) = Wis(A) Wise - 
Somit gilt 
Vian Wan) = Wis() Uae - (11) 


| Unter Benutzung von (9) und (7b) 148t sich weiter schreiben 
et.(4) = ©, +A Wan Une(A) = © +A Vea: Use Uea(A) = © +A Ves Wis(A) Vp. (12) 


G Damit ist die Berechnung der Kehrmatrix Ye» (A) auf die von ll,; (A) zuriickgefiihrt ; fiir i = 1 
» ergibt sich insbesondere 


Une) = Cn +A We aa Ua(A) Une « (12a) 
© Es geniigt somit, zur Aufstellung der Lésungsmatrix ©, in (10a) lediglich die Kehrmatrix 11,,(A) 
» zu berechnen; die iibrigen Elemente von © lassen sich dann, wie hier jedoch im einzelnen nicht 
| weiter gezeigt werden soll, von 1, (A) ausgehend durch je eine Matrizenmultiplikation er- 
s rechnen. 


3. Umformung der Matrizen fiir verschiedene Verwendungszwecke. Es wurde bis jetzt 
3) gezeigt, wie man, ausgehend von den Umsatzmatrizen lI;, die sich fiir die einzelnen Teilbetriebe 
eines Betriebssystems aufstellen lassen, auch den Betriebsablauf des ganzen Betriebssystems 
unter Beriicksichtigung der Kopplungsbedingungen durch eine Matrix darstellen kann. Wenn 
' eine solche Matrix nun auch einen weitgehenden Einblick in die Betriebsstruktur vermittelt, 
» so ist sie doch in dieser Form fiir manche praktischen Aufgaben nicht zweckmafig, weil sie 
’ einmal viele nicht weiter interessierende Einzelheiten bringt, andererseits aber auch aus ihr 
| ableitbare GréBen, die fiir die Beurteilung des Betriebsablaufes von Bedeutung sind, nicht 
/ unmittelbar enthalt. Da nun die Brauchbarkeit der Matrizenrechnung zum Teil gerade darauf 
beruht, daB méglichst viele Rechnungen vorbereitend erledigt werden kénnen, ist es notwendig, 
| die bis jetzt erhaltenen Matrizen weiter umzuformen. 

a) Bewertungen und Transformationen. Durch die Gleichungen (10) bzw. (10a) ist 
! der Betriebsablauf, d.h. Durchsatze, Ausbeuten und Fertigprodukte, als lineare Vektor- 
) funktion der Rohproduktvektoren t, gegeben. Nun ist es haufg so, dafs diese, wenn es sich 
\ beispielsweise um Gase handelt, wohl hinsichtlich ihrer Menge, aber nicht oder nur in geringem 
) Umfang hinsichtlich ihrer Zusammensetzung variiert werden kénnen. Kann man nun mit 
/ einer festen Zusammensetzung rechnen, so ist es wiinschenswert, diese mit in die Matrix 
) hineinzunehmen. Zu diesem Zweck fihrt man eine ,,Bewertung: der Rohprodukte gemaB 
R, a Dee Ts 

- ein, wo bi; einen Zeilenvektor, den ,,Bewertungsvektor“, darstellt. Setzt man beispielsweise 
Ge Oe oes 

- so stellt, wenn t, etwa ein Gas mit n Komponenten bedeutet, 


ee leds tl) ra, ea, ley ie * Fak 


To 1, 


Thk 
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einfach die Menge des Gases 1, (etwa in Nm?) dar und man kann von einer ,, Mengenbewertung“ 
sprechen, wahrend durch 

bi, b=1i(0 fel, wet) 
eine ,,Auswahlbewertung“ (etwa nach der GO -+- H ninth: was bei der Wasserstofferzeugung 
gelegentlich zweckmaBig ist) dargestellt wird. Von einer ,,Bewertung nach dem Heizwert* 
kann man bei dem Bewertungsvektor 


ad (itt Need hh 


sprechen, wo hy, h,,...h, die Heizwerte der einzelnen Gaskomponenten bedeuten. SchlieBlich 
kann man auch, z. B. bei wirtschaftlichen Rechnungen, eine gleichmafige Bewertung nach 
DM einfithren. Es sei jedoch darauf hingewiesen, da sich die hier eingefiihrte Bewertung 
nicht mit der sonst in der Mathematik itblichen deckt. Man kann nun weiter gemaB 


a 
w 
= 


+ — 6 
ie 


oe b>] 
r e 


,normierte** Vektoren fiir die Rohprodukte einfihren. Es gilt dann by, T,= 1 und es lassen 
sich dann gemafh 


Fatics To ONT R eT ; 
CS 0, Oe er ee i (13) 
Ty a Oe eee eile a) eee 


die in die Vektorgleichung eingehenden Werte 1, durch die Werte R, mit Hilfe einer im all- 
gemeinen nicht quadratischen Matrix © ausdriicken. Bemerkt sei noch, da man die Gleichung 


t;, = Th R, 


auch als Umsatzgleichung auffassen kann, wenn man unter R, etwa einen bestimmten Brenn- 
stoff versteht, der in einem bestimmten Betrieb zur Gaserzeugung eingesetzt wird. 1, bedeutet 
dann die aus der Gewichtseinheit dieses Brennstoffes unter den vorherrschenden Betriebs- 
bedingungen gewinnbare Gasmenge (als n-dimensionaler Vektor aufgefaBt). SchlieBlich kann © 
man auch noch die Matrix © mit der Matrix ¥$* fiir die Rohproduktverteilung (Ab- 

schnitt 2) verbinden, und erhalt dann unter Erweiterung der Definition fir © 


t, | =| Bh Wh - Tl Tie Oe 0] Re | ay apie tay R,| =6 [Ry | 

Vp WI WI. - War 0 4s 0) 4 a1 Ten * Vad R, | R, (13a) 
4 é 6 a 5 5 e 9 . 5 0 4 A | 3 

tm FO tee Bek 0 Le ai R, Ymitme ° Umi R; | R | ; 


wo in © = (t;,) die Elemente 1;, Spaltenvektoren 


—_—_—_——_-. 
= 
“mk 
= 
~ 
3 
= 


bedeuten. Auch fir die den Betriebsablauf zum Ausdruck bringenden Vektoren dD, a, f empfiehlt 
es sich vielfach, gemaB ) 


Dee Dy, dD; ° A= Da, A; ? f= by, fi 


Bewertungen einzufiithren. Dann wird 
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WeH(be > WAnw. ee | eee |) + 0) (b,)=2(d 

: : ; : : ; ; : “ 

Dn, 0 ‘ Don 0 i is 0 ‘ 0 8) f 
A, | Ca ee ree One ih ee gt Ig, 

; ‘ : ; : : : : : ‘ : 14 
An | 0 ° 0 0 : be .0 : 0 | ies a 
F, | 0 . 0 0 . 0 be, : 4 ts 
F,, | 0 130 ds 16.6 = ch Oripee © ~ Bin}. {fo | 


4 Es kann vorkommen, daf fiir eine GrifBe, etwa D;, mehr als eine Bewertung bendtigt wird; 
vedann setzt man, emer iecice bei zwei Hewertadeen, an Stelle von 


| 
i D;= bb, D; 
: 
5 
i 


peal en 
Dt) Na 


Sund erweitert damit die Matrix fT hinsichtlich der Zeilenzahl. 
Schreibt man die Vektorgleichung (10a), welche den Betriebsablauf in Abhangigkeit vom 
 Rohstoffeinsatz zum Ausdruck bringt, kurz 7 


y=Oz 


und den Ubergang auf die Bewertungen durch die »» Lransformationsgleichungen“ 


LCF bzw. p=, 


“einfach 


D 


g9=20 607 = Gr. 
© Geht man von der Vorstellung aus, daf durch die Transformationsgleichungen der Rohstoff- 
: einsatz ¢ bzw. der Betriebsablauf {) auf andere »,Koordinatensysteme“ mit im allgemeinen 
7 verringerten Dimensionszahlen bezogen werden, so kann man den Ubergang von der Matrix 9) 


«zur Matrix $$ gemaB 


B= T6 


in verallgemeinertem Sinne als Transformation mit den ,,Bewertungsmatrizen“* © und £& be- 
jzeichnen. Hierbei sind allerdings im allgemeinen © und & singular; die Matrix © schrumpft 
i bei der Transformation auf die Matrix S$ zusammen. Diese enthalt zwar nicht mehr so viele 
«| Einzelheiten wie ©, sie beschrankt sich vielmehr auf die fiir die Betriebsbeurteilung wesent- 
Hlichen GréBen und ist fiir die praktische Zahlenrechnung bequemer zu handhaben. Auch hier 
J zeigt sich wieder, wie zwanglos sich der Matrizenkalkiil auf die Untersuchung von Betriebs- 
vorgangen anwenden 1]aBt. 

In den Tabellen 1 und 2 ist die Transformation der in einer fritheren Abhandlung! gebrachten 
© Matrizen gezeigt. Dabei wurden folgende Gasanalysen zugrundegelegt?: 


Volumteile 
Einsatz in 


t 
t 
| 


co, co H, N CH 


Wassergas aus Koks Konvertierung | 0,07 | 0,40 | 0,50 | 0,025 | 0,005 
| Wassergas aus Grude | Konvertierung | 0,22 | 0,34 | 0,40 | 0,025 | 0,015 


Abgas aus einem Druckwasser- 
anderen Betriebs- wasche 
| system 0,30 | 0,07 | 0,60 | 0,02 | 0,01 
i Die gewahlten Bewertungen sind aus der Tabelle selbst ersichtlich. 


| b) Umkehrproblem. Durch die (unter Umstanden durch eine Bewertungstransformation 
Bob cewandelte) Gleichung (10a), ist der Betriebsablauf { als Funktion der zum Einsatz ge- 
ilangenden Rohprodukte 1 dargestellt gemaB 


fe St. 


1 T, Tabelle 2 u. 3. 
2 Siehe etwa Krénig, Die katalytische Druckhydrierung von Kohlen, Teeren und Mineralélen, 5.184, 


+ Berlin- Géttingen-Heidelberg 1950. 


| 
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Wassergas aus Abfa 
Koks Grude Abele e Koks Grude fellgas 
R R, R, R, R, R, 
Nm® Nm? Nm Nm Nn’ Nm? 
co+ H, co+ H, co+ H, co +H, Co+ H, co+ H, 
Durchsatz der Konvyertierung co, 0,0778 0,2973 — 0,0987 0,3186 0,0222 
co 0,4444 0,4594 _— 0,4953 0.5103 0,1053 
Dd, H, 0,5556 0,5406 — 0,5753 0,5603 0,0197 
(Nm?) N, 0,0278 0,0338 _— 0,0284 0,0345 0,0006 
CH, 0,0055 0,0203 — 0.0059 0,0217 0,0010 
Durchsatz der Druckwasserwasche co, 0,4722 0,7067 0,4478 0,5405 0,7754 0,5691 
co 0,0500 0,0500 0,1045 0,0535 0,0535 0,1108 
Dd. Jil 0,9500 0,9500 0,8955 1,0171 1,0171 1,0142 
(Nm?) N; 0,0278 0,0338 0,0298 0,0284 0,0345 0,0304 
cH, 0,0055 0,0203 0,0149 0,0059 0,0217 0,0159 
Durchsatz der Kupferwische co; 0,0200 0,0204 0,0202 0,0215 0,0219 0,0227 
(ore) 0,0485 0,0485 0,1014 0,0519 0,0519 0,1074 
Ds H, 0,9215 0,9215 0,8686 0,9866 0,9866 0,9837 
(Nm?) Ne 0,0270 0,0328 0,0289 0,0275 0,0334 0.0295 
CH, 0,0050 0,0183 0,0134 0,0053 0,0195 0,0143 
Ausbeute der Konvertierung co, 0,4722 0,7067 — 0,5405 0,7754 0,1213 
co 0,0500 0,0500 == 0,0535 0,0535 0,0063 
Oy H, 0,9500 0,9500 oo 1,0171 1,0171 0,1187 
(Nm) N, 0,0278 0,0338 — 0,0284 0,0345 0,0006 
CH, 0,0055 0,0203 — 0,0059 0,0217 0,0010 
Ausbeute der Druckwasserwdsche co, 0,0200 0,0204 0,0202 0,0215 0,0219 0,0227 
co 0,0485 0,0485 0,1014 0,0519 0,0519 0,1074 
A, 15 0 0,9215 0,9215 0,8686 0,9866 0,9866 0,9837 
(Nm!°) Nis 0,0270 0,0328 0,0289 0,0275 0,0334 0,0295 
CH, 0,0050 0,0183 0,0134 0,0053 0,0195 0,0143 
CO-Rickgas 
co, 0,0196 0,0200 0,0198 0,0209 0,0213 0,0222 
co 0,0475 0,0475 0,0993 0,0509 0.0509 0,1053 
M3 H, 0,0184 0,0184 0,0174 0,0197 0,0197 0,0197 
(Nm!) N, 0,0005 0,0007 0,0006 0,0006 0,0007 0,0006 
CH 0,0003 0,0013 0,0009 0,0004 0,0014 0,0010 
Ausgewaschene Kohlensaure 
co, 0,3577 0,5449 0,3380 0.4110 0,5985 0,4326 
co 0,0015 0,0015 0,0031 0,0016 0,0016 0,0033 
fe H, 0,0285 0.0285 0,0269 0,0305 0,0305 0,0304 
(Nms) N, 0,0008 0,0010 0,0009 0,0009 0,0010 0,0009 
CH, 0,0006 0,0020 0,0015 0,0006 0,0022 0,0016 
Frischgas co, = = nae = Soy, = 
, (ore) = = = —_ =, = 
ip H, 0,8939 0,8939 0,8426 0,95 70 0,9570 0,9543 
(Nm?) N, 0,0262 0,0318 0,0280 0,0267 0,0324 0,0286 
CH, 0,0046 0,0168 0,0123 0,0049 0,0179 0,0132 
Bei den praktisch auftretenden Betriebsaufgaben lautet die Fragestellung nun meist so, daft 


etwa die (im Betriebsablauf [ enthaltenen) Fertigprodukte | vorgegeben sind und die dafir 
notwendigen Rohprodukte rt und Betriebsbelastungen D ermittelt werden sollen. 
Beantwortung dieser Frage ist der Rang der Matrix ©, der héchstens gleich der Anzahl der 
Komponenten von ¢ sein kann, entscheidend. Fiir die Vektoren [ und t, sowie fiir die Matrix 0 


Tabelle 1. Transformierte Matrizen. 


Matrix US 


Unverzweigte Betriebskette 


Betriebskette mit Rickfiihrung 


Wassergas aus 


werde nun, u. U. nach geeigneter Umordnung von Zeilen und Spalten, gemaB 


oder ausfiihrlicher geschrieben 


= ol 


tee | a 


if 


eu 
1 


= 
11 


Pes |r 


Dos — Dior Dig? Dao 


el sl 


Dor Dep | 


mz Div steeds 
- Dio T + Dot 


gilt. Daraus erhalt man durch eine einfache Zwischenrechnung : 


a Oi Duro 


i = Ou Dio 
t Dio Doe 


eine Aufspaltung derart vorgenommen, daf die Zahl der Komponenten von { und f, sowie der 
Rang von 2); gerade gleich dem Rang von © sind, und da die Beziehung 


{ —= Dh Die | 


Si —l 


ell | 


Fiir die 


(15) 
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Tabelle 2. Transformierte Matrizen. 
Matrix 8 =ITOS 


Unverzweigte Betriebskette Betriebskette mit Rickfiihrung 

Wassergas aus iWiseenens aus 
| sserg 8 

Koks Grude pesliges Koks Grude py ialiege 

i R, R, Re R, R, R, 

Le Nm? Nm? Nn? Nm? Nm? ; Nm? 

‘ cO+H, | CO+H, | CO+H, ] CO+H, | CO+H, co +H, 
»)urchsatz der Konvertierung Ds» Nm? ner Ue ia 1,3514 — 1,2036 1,4454 0,1488 
yurchsatz der Druckwasserwasche Des Nm’ 1,5055 1,7608 1,4925 1,6454 1,9022 1,7403 
‘vurehsatz der Kupferwische D;* li Nm? 1,0220 1,0415 1,0325 1,0928 11133 1,1576 

Dd 

durchsatz der Konvyertierung ‘DS Nm*CO-+H, 1,0000 1,0000 — 1,0706 1,0706 0,1250 

‘urchsatz der Druckwasserwiasche DE Nm'CO-+ H, 1,0000 1,0000 1,0000 1,0706 1,0706 1,1250 
qurchsatz der Kupferwasche Dees Nm°CO+H, 0,9700 0,9700 0,9700 1,0385 1,0385 1,0911 
dusbeute der Konvertierung Bite Nm’®CO+H, 1,0000 1,0000 — 1,0706 1,0706 0,1250 

a lusbeute der Druckwasserwische idan = Nm'CO-+H, 0,9700 0,9700 0,9700 1,0385 1,0385 1,0911 
1 O-Riickgas ast Nm®CO+ H, 0,0659 0,0659 0,1167 0,0706 0,0706 0,1250 
10- Riickgas Pe a le 10° WE 0,1932 0,2018 0,3521 0,2077 0,2163 0,3770 
Lusgewaschene Kohlensaure Oe | Nm? 0,3877 0,5779 0,3 704 0,4446 0,6338 0,4688 

a rischgas See Nm®CO+H, 0,8939 0,8939 0,8426 0,9570 0,9570 0,9543 

Matrix © Bewertungsvektoren 
R, R, R, | co; | COMME, pen, CH, 

i 0,0778 0,2973 — Fir Nm?’ 1 1 il 1 1 
0.4444 0,45 94 = » Nm?CO+H, 0 1 1 0 0 

it, 0,5556 0.5406 — » 103 WE 0 3,017 2,574 0 8,574 

0,0278 0,0338 —— 
0,0055 0,0203 — 

| — — 0.4478 

| = = 0,1045 
ty Es a 0,8955 

6h a as 0,0298 

— — . 0,0149 


| ? Zur Bezeichnung D,* usw. s. S. 163. 


o 7 existiert, sofern der Rang von Q,> 0 ist. Mit dieser Gleichung kénnen die aus dem 


" 

| T 
_|, der vorgegeben werden kann, folgenden Komponenten des Betriebsablaufes, zu- 

t ie 

4 

vammengefaBbt in dem Vektor ral berechnet werden. Fiir den praktischen Gebrauch ist es 
| 

‘un wiinschenswert, in dem auf der linken Seite stehenden Vektor samtliche Rohstoffe, 
Durchsatze und Fertigprodukte zusammenzufassen, also auch diejenigen Komponenten des 


etriebsablaufes, die auf der rechten Seite stehen und vorgegeben werden kénnen. Man 


ob 

@rreicht dies in einfacher Weise dadurch, da man fiir diese Komponenten in die Matrix { 
| Jeilen einschiebt, die an dem Schnittpunkt mit der entsprechenden Komponentenspalte eine 1, 
shonst aber als Elemente 0 haben. Die so erweiterte Matrix werde mit §§ bezeichnet. Es ist 


| ih i : = See 
serner zweckmafig, in dem rechts stehenden Vektor |_| die Fertigprodukte [, soweit sie vor- 
Y 


| 
‘segeben werden kénnen, besonders herauszustellen und die iibrigen Komponenten als betrieb- 
‘ich gegebene Nebeubedingungen zu einem Vektor 1 zusammenzufassen, so daB also 


i 


1 


ay 
I 


t 
i] 
| 
| 
| 


silt. Dann 1aBt sich der Betriebsablauf eines ganzen Betriebssystems (jetzt in abgeanderter 


form) dargestellt durch den ,, Betriebsvektor“ 


n Abhangigkeit von dem Produktionsprogramm, dargestellt durch den ,,Produktionsvektor“ 
12 
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(4) daiok 


wiedergeben. | ' 

Durch die oben erwahnte Ranguntersuchung der Matrix 0, die in der Matrix ‘$ ihren 
Niederschlag findet, wird die Struktur eines Werkes, d. h. eines Betriebssystemes, weitgehend 
analysiert. Die Matrix $$ mége daher auch als ,,Strukturmatrix“ bezeichnet werden. Wesent- 
lich ist, daB die Auswahl der Vektoren f und 11, und damit auch die Matrix ‘$ im allgemeinen 
nicht eindeutig sind. Fiir die Auswahl der verschiedenen Méglichkeiten spielen Kapazitats- 
betrachtungen und wirtschaftliche Uberlegungen eine Rolle1. Zuweilen ist es notwendig, fir 
mehrere mogliche Falle Strukturmatrizen aufzustellen, um die verschiedensten, praktisch auf- 
tretenden Fragen rasch beantworten zu kénnen. Bei der Durchfithrung der damit verbundenen 
numerischen Aufgaben hat sich bei umfangreichen Matrizen eine Aufspaltung in kleinere 
Matrizen als zweckmaBig erwiesen?. 


c) Umkehrproblem bei Zwischenproduktentnahme. Das Umkehrproblem, d. h,| 


(16)| 


der Ubergang von der Darstellung {= 0 r auf die Darstellung {= % | i | , 1aBt sich ohne > 
Tl 


weiteres lésen, sobald die Matrix ©) mit allen ihren Elementen gegeben ist. Nun treten bei 
konformen Kopplungsbedingungen in der Matrix $$ die Elemente der Verteilungsmatrix 


c (s. Abschnitt 2) auf und es gibt sehr haufig Falle, wo diese 
, Matrix nicht von vornherein bekannt ist. Die dabei auftreten- 
eg den Fragen sollen zunachst an einem einfachen Beispiel, der 


di=r+A303 | Betriebskette mit teilweiser Rickfihrung’* naher erldutert 
werden. Es kommt vor, daB die Ausbeute der letzten Stufe einer 
Betriebskette (etwa das CO-Riickgas bei der Wasserstoffherstel- 
lung?) nur teilweise zuriickgefiihrt wird, da ein bestimmter Teil 
dieser Ausbeute an anderer Stelle benétigt wird, oder aber der 
Singularitét von (© — Un) begegnet werden soll4. In dem im 
Abb. 2 gezeichneten Beispiel mit drei Betrieben werde von der 
Ausbeute a, nur der Anteil j,a, zuriickgefiihrt, wahrend der 
Anteil 


33 = (1 —As)a, 
anderweitig als ,,Zwischenproduktentnahme“ verwendet werde. 
Lal (Durch die fiir den Ansatz von 3, gewahlte Form kommt zum 
Ausdruck, daB 3, dieselbe Zusammensetzung wie a; haben muB.} 
is-(MAs)a, Die Rechnung geht genauso wie bei der Betriebskette mit Riick- 
Abb. 2. Betriebskette mit teilweiser fihrung® vor sich und man hat zunichst 


uckfihrung. 
a; = Us; Uy, Uy d, = Ug, d, . 


Fir den Durchsatz 0, gilt aber: 
Dy =f1 +A, O13 = 1 +h, VES Die 


(© —A, Ug) dD) = 1 
D, = (C—A, Ug) 1 1. 


oder 
und weiter 


Fiir f; ergibt sich damit 
fs = Bai (E — A, Uys) t mit Bs = Bz Uy, 
wahrend fiir die Zwischenproduktentnahme 3, folgt 
3 = (1 —A,) Us (© —A, U3) 2. (17) 


1 J, Seite 135. 

2 Vel. etwa H. Unger, Z. angew. Math. Mech. 32 (1952), S. 1. 

Ss Siehe 1, Ss 119: 

eT Sicherl eo aloo 

° Die folgenden Gleichungen ergeben sich aus der ersten Spalte der Matrix (10) fiir aii. — se 
1 =4,; vel. auch I, S, 125, a | 
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Fir A, = 0 bzw. 1 ergeben sich die unverzweigte Betriebskette, bzw. die Betriebskette mit 
(| Ricklauf. Praktisch ist es nun meist so, daB A, nicht vorgegeben ist, sondern etwa die Gesamt- 
}menge von 4, oder sonst ein mit %, zusammenhangender Skalar Z;. Fiithrt man durch 6;, eine 
» Bewertung ein, so ergibt sich fiir 2, die skalare Bestimmungsgleichung , 

| Zs = 0s, 33 = (1 —As) b5, Ug, (E — Ay Us) v, (18) 
fieine Gleichung, die im allgemeinen nicht mehr linear ist. Fiir die Vorgabe von Z, gilt dabei 
§die cinschrankende Bedingung, daB sich fiir A, eine zwischen 0 und 1 liegende Lésung, die 
Grenzen mit eingeschlossen, ergeben muB, wenn die Zwischenproduktentnahme Z, vealisierbar 
‘sein soll. So wird man bei Problemen dieser Art, ahnlich wie bei Kigenwertproblemen, auf 
Snicht lineare Gleichungen gefiihrt, deren Lésung vor allem dann, wenn auch noch an anderen 
iStellen Zwischenprodukte entnommen werden, sehr miihsam werden kann. 

Der Unterschied gegeniiber dem Eigenwertproblem liegt darin, daB bei diesem die charak- 
jteristische Gleichung itber die charakteristische Determinante gewonnen wird, wihrend man 
beim Zwischenproduktproblem die Skalargleichung fiir 2 durch Multiplikation mit einem 
9 Bewertungsvektor b; erhalt. Dieser Unterschied tritt bei vektorieller Darstellung noch deut- 
Wlicher zu Tage. Setzt man in eine Betriebskette mit Riickfithrung an Stelle des Rohprodukt- 
H’vektors r den Vektor (t — 33) ein, so gilt (fiir m — 3) 


7 M3 = Us; (© — Ug,)~t (t — ga) = P(t — 35) - (19) 
[Soll 3, nun ein von a, abgezweigtes Zwischenprodukt sein, so mu® 


$3 = (1 — As) ag 


sein, wobei j, durch Z, = bj, 3, festgelegt wird. Setzt man 


a, =0 und (h==2,) == 15 
so wird aus (19) 
a=A$a+ Br, 
jeine Gleichung, die fiir den hier bedeutungslosen Fall, dafi r = 0 ist, in die Eigenwertgleichung 
C= AAO 


mit dem Eigenwert A iibergeht.1 Da aus Gleichung (19) und g, = (1 — A;)a, die Gleichung (17) 
folgt, 1aBt sich auch rein formal, wie folgt, sehen. 
Zunachst 1a46t sich fiir (19) unter Beriicksichtigung der genauso wie Gleichung (9) nach- 
| zuweisenden Identitat 
(€ — Uy)? = € + Uy (E— Ug) = E + (C— Uy) Uy 
schreiben 
a, = (© — Ug) Us: (YX — (1 — Ag) Ms) : 
'Daraus folgt 
[© + (C — Wg:)4 Usy — Ag (E — Ugy)~* Usy] a3 = (C — Uys) Ug ¥ 
'oder unter Beriicksichtigung, da®B gemaB der oben gebrachten Identitat 


E + (€C — Ug) Ug, = (EC — Uy) 


(© — U,1)—4 (C — Ag Usps) Ag = (© — Uy,) 4 Ug t . 


 Daraus folgt aber 
ds = (€ — Ay Ug)? Ug t = Ugy (© — Ag Ua) Tt, (17) 
} und dies ist, da 3, = (1 —A,) a, gleichbedeutend mit Gleichung (17). 

| Zur Behandlung des allgemeineren Falles einer Betriebskette mit Rohprodukteinsatz und 
-Zwischenproduktentnahme an jeder Stelle werde von den Gleichungen (10a) ausgegangen und 


fi 
|; 
| 
4 


Ras ber, | D; = 05, 0 » A; = be, 0; und F, = by, fi 
» bewertet. 


/ 1 Es sei in diesem Zusammenhang jedoch darauf hingewiesen, da® die Behandlung von Anfahr- und 
| Abstellvorgangen bei chemischen Betrieben auf Eigenwertprobleme fiihren kann. 


12* 


168 Pichler: Anwendung der Matrizenrechnung zur Erfassung von Betriebsablaufen, _Ingenieur-Archivyy 


Mit den normierten Vektoren 


den Abkiirzungen 
dj(A)= 65, Uis2 Uar(A)t. fir ki, d= bp, WA) Tt ‘| 
a;4(A)= be, Wie Une (A) Th (i,k =1,2,...m), (20))) 
fie(A)= Bf, Bin Wan (A) Te (i,k =1,2,...m), | | 


und weiter mit 


ie a R, > D= D, > a = A, r) {= F, 
15 Dd A, F, 
| Re | D,, As F., | 
folgt aus den Gleichungen (10a) 
a= ay,(A) Am %42(A) : Amm ®m(A) | ¥ « 
Ay 4(A) (4) Aro @am(4) (21) 
| Am —a, 1 @m(A) Am—1,2 %m9(A) : Om m(A) 


Wenn man a;;,(A) durch d,,(A) bzw. f;,(A) ersetzt, ergeben sich ganz analoge Gleichungen}y 


fiir D und f. Fir J, =—2 = 0 kommen die Matrizen fiir die unverzweigte Betriebskette mit } 
den Elementen d;, (0), a;, (0) und f;, (0), wahrend aus A; =A = 1 die Betriebskette mit (voll-_ 
kommener) Riickfithrung und den Matrizenelementen d;, (1), a;, (1) und fj, (1) folgt (siehe 
Zahlenbeispiel in Tabelle 2). Wendet man fiir die Normierung von 1;, 0; und aq; dieselben 
Bewertungen 


b’= b,, = bp, = ba, 
an, so erkennt man unmittelbar aus den Definitionsgleichungen fiir d;;, (A) und a;, (A), daB dann 
d;,(A) = 4,_, (A) fir 2 == kh ist, 
wahrend fiir 1 = k unter Benutzung der Gleichung (12) 
d;:(4)= 6’ Wied) t = B (C +A U4: Wii(A)) 
= t+ Ape eae Ue Ae 
=1+/4_,,; 
=1 +4 41 G13 
folgt. Diese Beziehungen sind eine Folge der Kopplungsbedingungen (s. Abschnitt 2) 
d; = A1 Oe eG + 
die sich bei iitberall gleicher Bewertung in 
D; =4,_4 A_, + RB; 
abwandeln. Fiihrt man auch fiir die von den Ausbeuten a; abgezweigten Zwischenprodukte bi 
die Bewertungen 
bss = Dai 
ein, so folgt aus 3; = (1 —4,) a; 
Z; = (1 —A,) 4;, 
und damit erhalt man bei vorgegebenen Zwischenproduktentnahmen Z; fiir A; aus (21) 
m Bestimmungsgleichungen 
Zi * 
ap jaar es >, hawt a;4(A) Ry + a;;(A) R;, (22) 
k=1 
kei 
wobei Aj, =A; A,1...4, (zyklisch aufgefaBt) und A =/,_1;=A,_14 =Am Arn-1.0 «(jae 
Dies sind m algebraische Gleichungen fiir J; , deren unmittelbare Behandlung bei numerischen 
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3, Rechnungen im allgemeinen sehr miihsam ist. Bemerkenswert ist jedoch der Sonderfall, daB 
R, = 0 fir k> list. Dann vereinfachen sich die Gleichungen, wie folgt: 


A, = 4,,(A) R,, 

A; — aA a;,(A) R, (i = 1) A 

| a oa = Vie ae e Eee a 

oo" man A; = qty? a ay ein, so wird Z; = (1 —4A,) A; und man erhalt folgende 
} Bestimmungsgleichungen fiir /;: 

Z, = (1-4) 4, =(1— dn) Ri, 

{ Ly = (1 — dy) Ay = (Anaya — Aur) Ry (l<k<m). 


Daraus folgt durch Addition der ersten i Gleichungen 


und speziell fiir 1 — m 


DP = De ay OR, 


uw! eine Bestimmungsgleichung fiir A. Mit A sind dann die GréBen a,,(A) , sowie d,,(A) und f;,,(A) 
Sbekannt. Fiir die Bestimmung von 4,;, D; und F; erhalt man unter Beriicksichtigung, daB 


e;aus (23) . 
An te R,— >)% iv Cha ae 
k= kal Pl 
folgt 
i—l 
ee eal) A; = 4;,(2) Ry — S87 (Evie 
kh=1 
ay (a) 
D, = (2) Ry D; = dis(2) Ry — D'S Gy Zi (Cs Qe 
kt 
S fil) 
F, = fulA) R, Fy = fix(A) R, =i De a,(1) Zi, (i > 1) : 
kw 


i Ist Z; = 0, dann ist A; =A =1 und man erhilt die Betriebskette mit vollstandigem Riicklauf. 
Ist dagegen Z,, = A,,, so erhalt man die unverzweigte Betriebskette. Es wird dann mit An 
‘auch 4 = 0, und an Stelle von a;; (A) usw. treten die entsprechenden GréBen a;; (0). Sind 


m—1 


simtliche Werte Z; vorgegeben, so folgt aus An = Ln = R,— >» Ze, 
i k=1 


} 
1 
| 
} 
Nt 


‘und damit 


ay,3(0) 


A= Sl) z (24) 
va 


=F 


|, Mit den entsprechenden Gleichungen fir D,; und F, ist damit fiir diesen Sonderfall (R, = 0 
| fir k> 1) der unverzweigten Betriebskette bei vorgegebenen Zwischenprodukten Z, das Um- 
4 kehrproblem gelést und zugleich gezeigt, daf sich der Rohproduktbedarf R, linear in Z; dar- 
5) stellen 14Bt. Ist R, auch fiir k> 1 von Null verschieden, so lassen sich, wie hier jedoch nicht 
3} weiter ausgefiihrt werden soll, die Werte J; in einfacher Weise und zwar rational aus R, und 
' Z, berechnen; dagegen ist bei vorgegebenen R, (k> 1) und Z, die GréBe R, im allgemeinen 
nicht mehr rational darstellbar. 

Unter diesen Umstanden ist es naheliegend, nach geeigneten Naherungsmethoden zu 
; suchen. Dazu geht man zweckmabig nicht von den Bestimmungsgleichungen (22) fir A; aus, 
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sondern von dem Ansatz (man vergleiche auch den entsprechenden Ansatz in Gleichung (19)) 


Ni 2 Uae Ve (T. — 34-1) + Ui (t; — 3-4) 
nes 

Ce > Win War (te — dea) » (25) 
b= 


i= >. Bin Wen (oe — 3a) - 
ail 


Man iiberzeugt sich leicht, daB die Umsatzgleichungen 
qd = U; D; und fi == B; D; > 
sowie die Kopplungsgleichungen 
De =A a Gy +t = Oy + (ti — BH) z 
durch diesen Ansatzerfiillt sind. (Dabeiist wieder zu beachten, daB gemaB (12) UA ye C=; 


ist.) Dazu kommt jetzt noch die Forderung, da8B 3; und a; dieselbe Zusammensetzung haben. 
Bei Einfithrung der Bewertung bj, und der normierten Vektoren 


A 5 oh a oe ne pee 
ae 04; 8 es ape i bi; @, Wirsi 
bedeutet dies, daB 
3: = 0; . (26) 


oder 
she Br 7h tS 
bi Oi L Pe a; L 
fiir jede Bewertung bj, sein mu8. Geht man nun von einer geeigneten Annahme fir 3; aus, 
so laBt sich damit zunachst etwa D, berechnen. Daraus folgen aber zwangslaufig unter Be- 
nutzung der Umsatzgleichungen a; = ll; 0; und der Kopplungsgleichungen 
Zi—-1 
d= (1 = faa) GG 
bij Gi} 
die itbrigen D; und schlieBlich ein neuer Wert fiir },. Damit ist ein Iterationsverfahren gegeben, 
das sich an durchgerechneten Beispielen gut bewahrt hat, dessen Konvergenzbeweis aber vom 
Verfasser noch nicht erbracht ist. Wenn der Wert R, = by, 1, die Werte R, stark itberwiegt, 
was praktisch in vielen Fallen zutrifft, so erweist sich insbesondere der Ansatz 
. 5 Wer Ya 
aps he 27 
teeraigy an et “a 
als vorteilhaft. Es lat sich namlich zeigen, daB, wenn R, = 0 fiir k> 1 ist, dieser Ansatz 
nicht nur fiir Z, = 0 eine exakt richtige Lésung bedeutet, was ja trivial ist, sondern auch 


noch ftir den Fall, daB 3,, = 4,, ist, die Betriebskette also in eine unverzweigte Betriebskette 
ausartet. Mit diesem Ansatz wird fir R, = 0 bei k> 1 aus (25) 


ke me —— 4 ave 3 Zi 
a; = (Wis Un ty) Ri — Y (Uses Weyer Ue t)%, omit =. (28) 
k=1 am * e111 
Nun gilt gemaB Gleichung (11) 


We saepa Na = Na Ua : 
Ferner gilt 


Lies eg Veg fiir a is 
dagegen a | eg ems ir ee ite | 4. 
Daraus folgt aber, daB 
Ue ts cts Wn = Vi, Wea fiir i > k, 
Uses 1 Fale anar Ver == iy 11 a Vier 
a WinMa — &) 


Se Ve Wis fiir i <k 
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ird (gema® (12), fiir k = 1 und A= 1). Damit wird aus (28) 


m 


a= (Ui, Cire tT) (R, pe hah (Ui; ¥. jee 


k=1 
und insbesondere 


Be = (Uni Wy tH) we ee, 2, Ge eat) ea 
k=1 
Ist nun 3, = Q, (Ausartung in die unverzweigte Betriebskette) und damit auch 


' 
b Om 


Qin = Lin ° 
‘so folgt aus 


m 


| A Ln = ee na hike T,) (R, san 32,) ++ La ’ 


k 


wenn man den trivialen Fall 6) Un U,,T, = 0 au®er Betracht 14Bt, 
135 eae) Zi 
k=1 
jund damit 
a; == (U5 ty )& Za» 


jsowie unter Beriicksichtigung von (27) 


‘d.h. die Forderung (26) ist erfiillt. Funee man in eee ciate mit (20) die normierten 
Vektoren 


Wi Unt — Ua ty Y 
bag Uy Ded. a;,(1) 


_ Yat _ Yat 
und 4; (0) = G Usr% 410) 


Aj; (1) = 


fein, wobei also 


| Bj, di (1) = Bo; G1 (0) =1 
rilt, so 1ABt sich 


m 


a, = a1 (1)(R— 3%) 0,0) + (0) FZ - a5, 0), 


i 


ae re (29) 
A; aa), a; =— 1) (2, 34%) + 4, (0 ) YZ 


k=i 


schreiben. Die letzte Gleichung geht fiir R, = »' Z, (Ausartung in die unverzweigte Betriebs- 
| k 


=1 
‘kette) unter Beachtung von (28) in Gleichung (24) iiber. 
Fiir den Fehler der mit dem Ansatz 
oy Wii t AI 

fas oS =; 0 

i) b/ SiG: C1 i ) 
durchgefiihrten Naherungsrechnung gibt die Differenz (a; — 3:) Z;, die ja bei der strengen 
Lésung verschwinden muB8, ein geeignetes Ma. Nun ergibt sich aus (29) 

4 ix (1) (x, a E%)-a a; (1) + 454 (0 ) F250) ) 


ae Cu 


k=1 


4;,(1) (x, eo 5%) oF a;1(0) SZ 


Somit wird 


(G; — 3:) Z; = (G; — Oi (0) Z; = (ai1(1) — 1 (0)) 


172 Pichler: Anwendung der Matrizenrechnung zur Erfassung von Betriebsablaufen. _ Ingenieur-Are 


Beriicksichtigt man, daB gemaB® der Bedeutung der Umsatzmatrizen 4;, (1) und a;,(0)> 
sind, so ergibt sich die Abschatzung 
(a; — 3) ZS (Gi (1) — 21 (0)) Z- 

Hat man fiir die unverzweigte Betriebskette und fiir die Betriebskette mit (vollstandiger 
Riickfiihiung die Matrizen mit normiertem Rohprodukteinsatz 1, berechnet, so lassen sic 
daraus die Vektoren q;, (1) und a,; (0) ohne weiteres entnehmen. Fiir das in Tabelle 1 durch 
gerechnete Beispiel ergibt sich beispielsweise bei einem Einsatz von Kokswassergas (cO + H, 
Bewertung) 


Gy, (1) = [ 0,0215 | , Gp, (0) = | 0,0200] . 
0,0519 0,0485 
0,9866 0,9215 
0,0275 0,0270 
0,0053 | 0,0050 | 
Bei einer (CO + H,)-Bewertung fiir a, ergibt sich daraus 
Se Ee = Tne Oar) seer 02 06 
Q (1) =T 9385 =| 0.0500] “24 %1(°) =9,5700 = | 0.0500 
0,9500 0,9500 
0,0265 0,0278 
[ 9.0051 | 0,0052 
und damit 
Gy (1) — C= | ELOOCOL 
0,0000 
0,0000 
— 0,0013 
— 0,0001 | 


Werden aufer R, auch noch an anderen Stellen Rohprodukte R, eingesetzt, verschwindet 
also R, nicht fiir alle k> 1, so treten in der obigen Gleichung fiir (a; — 3 ;) Z; noch weitere 
Glieder mit R, auf. Praktisch erweist sich aber der gewahlte Ansatz 3; = @;, (0) auch in solchen 
Fallen als sehr brauchbar und zwar einmal, weil R, meist iiberwiegt, dann aber auch, weil | 
sich die Zusammensetzung von t, wenn ein Kinsatz an der Stelle k tiberhaupt méglich sein 
soll, nicht zu sehr von der Zusammensetzung von (,_, unterscheiden darf. Die betriebs- 
mafigen Bedingungen, die an %, fiir k> 1 gestellt werden miissen, tragen daher haufig dazu 
bei, daB durch den obigen Ansatz brauchbare, in R, und Z, lineare Naherungslésungen méglich 
sind. Genau wie fiir A; in Gleichung (29) lassen sich auch fiir D; und F; in solchen Fallen 
mit guter Annaherung in R, und Z, lineare Darstellungen geben, die sich etwa in der Ma-. 
trizengleichung | 
oe 5G) 

a 6 

Gi 

zusammenfassen lassen und die dann weiteren Umformungen, wie sie im Abschnitt 3b be- 
handelt wurden, leicht zuganglich ist. In diesem Zusammenhang midge auch noch kurz auf 
den Fall’ eingegangen werden, bei dem die Ausbeuten durch 


a; = N;0;— a; t 
dargestellt sind. Setzt man a* = 11; D;., so wird 


A; = AF — Agi t 
und die Kopplungsbedingungen 


dD; = G4 + (tr— 44) 


gehen iiber in 
Dp = O74 + (t; — eeu ee) 
= Oh_it (t7— 34), 


a Srehenl oy Loe 
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wo 


| 


gesetzt ist. Damit ist das Problem, wenn der Vektor Qj; vorgegeben ist, auf das bereits be- 
handelte zuriickgefiihrt. Ist der Vektor a); dagegen in der Form 

Ag; = Agi 0; 
-vorgegeben, wo wieder 


ist, so setzt man 


ai ai a ON a (Z; + Ap;) a; 


und hat dann das Problem ebenfalls zuriickgefithrt. Ahnlich verfahit man, wenn fitr f; die 


: 


erweiterte Darstellung 


fi = BD; — fot 


angewandt werden muf. 


| 


Was schlieBlich die Behandlung des Problems bei allgemeinen Kopplungsbedingungen 
anlangt, so ist der Ansatz dazu ohne weiteres aufzustellen, wenn man die Gleichung (3) in der 
Form 


j= 20 
schreibt, wo & eine Diagonalmatrix mit den Elementen ©, — SS YS; ,, bedeutet. Unter Be- 
nutzung der Gleichung (6) wird daraus: J 
$= LU(C— WU) BW tr. 
Fiihrt man Bewertungen ein, so erhalt man daraus skalare Gleichungen fiir die Elemente der 
Matrix YW, die den Gleichungen (22) entsprechen. Wa&ahrend dort fiir praktisch brauchbare 
Lésungen fiir A; die Schranken 
OFA 

bestehen, gilt jetzt neben 

OF S057, = 


auch noch 


i=l 
Dazu treten praktisch noch eine Reihe weitere Bedingungen w;, = 0 fiir bestimmte 1 und k, 
da es praktisch im alJgemeinen nicht méglich ist, jede beliebige Ausbeute a in jeden be- 
liebigen Betrieb i einzusetzen. Gerade diese Einschrankung erlaubt es, auch in anderen als 
den oben behandelten Sonderfallen ohne iibermaBbigen Zeitaufwand Lésungen zu finden. 


d) Umkehrprobleme bei durchsatzbedingten Kopplungen. Bei dem Umkehr- 


problem bei Zwischenproduktentnahme sind gewisse mathematische Schwierigkeiten durch 
die Forderungen entstanden, da8 samtliche Komponenten der Ausbeuten a in gleicher Weise 


verteilt werden miissen. Nun gibt es aber viele Falle, wo diese Forderung nicht besteht und 


sich dann das Umkehiproblem auch bei allgemeiner Kopplung sehr einfach lésen 1aBt; es 


kommt auch praktisch vor, daB etwa bei Gasen oder Fliissigkeitsgemischen diese Forderung 
durch technische MaBnahmen, z. B. durch Einschalten einer Gastrennungsanlage oder Destil- 
lation, beseitigt wird. Zur Darstellung des Umkehrproblems in solchen Fallen werden nun — 
in teilweiser Abanderung der fritheren Bezeichnungen — samtliche im System vorkommenden 
Ausbeuten und Rohstoffe durch die skalaren GréBen A; und Rj, bezeichnet, wo der erste 
Zeiger i bzw. j die Art der Ausbeuten bzw. Rohstoffe und k den Betrieb, in dem die Ausbeuten 
anfallen oder in dem die Rohstoffe eingesetzt werden, bedeuten. Dabei sollen fiir A;, auch 
negative Werte zugelassen werden, und ein solcher negativer Wert soll bedeuten, daf in dem 
Betrieb k eine bestimmte Menge der Ausbeutesorte i eingesetzt wird. Es werde nun — vor- 
behaltlich einer nachstehend noch zu bringenden Erweiterung — angenommen, da der Ab- 
lauf des Betriebes k durch einen skalaren Parameter, den Durchsatz D,, eindeutig beschrieben 
und insbesondere 
Ai, = win D, » Rix =i D, 
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gesetzt werden kann, wo unbeschadet ihrer geanderten Bedeutung die Verteilungskoeffizienten 
aus Abschnitt 2 tibernommen sind. Der Durchsatz D, ist bei dieser Auffassung nicht gleich- 
bedeutend mit dem Einsatz; er kann es gelegentlich sein, er kann aber ebenso gut mit einem 
bestimmten A;, oder einer linearen Kombination aus solchen zusammenfallen. Man kann 
ihn in gewissem Sinn als ,,verallgemeinerten Durchsatz‘‘ bezeichnen, ahnlich wie man in der 
Mechanik von einem verallgemeinerten Impuls usw. spricht. Dagegen ist es notwendig, die | 
obigen Ansatze zur Beriicksichtigung von Verlusten und sonstigen durchsatzunabhangigen 
Einfliissen, wie folgt, zu erweitern: 


A;, = w;, D, + wirt, 
R;, = wy, D, + wit - 


Dabei bedeutet t einen vektoriellen Parameter (Spaltenvektor), dessen Komponenten die 
verschiedenen den Betriebsablauf beeinflussenden GréBen, z. B. Lufttemperatur, darstellen. 
Seine Dimension sei so gewahlt, daB sie die KinfluBgréBen samtlicher Betriebe umfaBt; tv, 
und tyj% sind Zeilenvektoren derselben Dimension. Meist wird t die Zeit t enthalten (s. Ab- 
schnitt 3c den Ansatz a; = Ul; D0; — a9; t); dann bedeuten D;,, A;, und R;;, die Durchsatze, 
Ausbeuten und Rohstoffeinsatze in der Zeit t. 

Die itiber k genommene Summe »' A;;, stellt nun den das aus m Betrieben bestehende 
System verlassenden Ausbeutetiberschu8 der Sorte i dar; er soll im Gegensatz zu den 
fritheren Abschnitten, wo er mit Z; bezeichnet wurde, jetzt als Fertigprodukt F; bezeichnet 
werden. Desgleichen bedeutet R; = 2' R;, den gesamten in das Betriebssystem eingesetzten 
Rohstoff der Sorte j7. Damit 1a8t sich aber unter Benutzung von (30) schreiben 


(30) 


ES 64 De +( 5 wha) t 
k=1 k=1 


Ry LeooRY ohne wD, +( 3 wit) t 
k=1 k=1 
Mit 


XS = (w;,) (n Zeilen, m Spalten) Wt = (n Zeilen, s Spalten) , 


YW* = (wF;,) (1 Zeilen, m Spalten) WF = | (1 Zeilen, s Spalten) , 


wo n die Anzahl der Ausbeutesorten, 1 die Anzahl der Rohstoffsorten, m die Anzahl der Be- 
triebe und s die Anzahl simtlicher Kinflu8gré8en (Dimension des Vektors t) bedeuten, sowie 
mit 


f=|f a Sree a dD=|D, 
F, R, D, 
ie R, Dry 


1aBt sich nun 

| MB* WE] \t t 
schreiben und daraus erhalt man, wie in Abschnitt 3b gezeigt wurde, durch entsprechende 
Umkehrung wieder die Gleichung (16), 


CH) 
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wo nun die Komponenten von t unter denen von tt enthaltem sind. Es sei noch bemerkt, daB 
sich auch bei den in den vorigen Abschnitten behandelten, durch 3 =a gekennzeichneten 
Fallen eine in D, lineare Darstellung wenigstens niherungsweise erreichen lie®; dort war 
dies aber erst nach einem langeren Umweg méglich, wahrend hier von vornherein von einem 
solchen Ansatz ausgegangen werden konnte. 

Es kann vorkommen, da ein skalarer Durchsatzparameter D, (auch neben dem Parameter f) 
zur Beschreibung des Betriebsablaufes im Betrieb k nicht ausreicht. Dann ersetzt man D, 
durch den vektoriellen Durchsatzparameter }, und ersetzt (30) durch 


Ai,= Vin De + Wie t >| 
Riu= biid, + wit t.J 


Es haben dann YW = (bj,) und Q* = (p/{) mehr Spalten und ) mehr Komponenten. Sonst 
-andert sich nichts an der Ableitung. 

| Als Beispiel fiir ein Betriebssystem, das sich in der beschriebenen Weise behandeln 1a8t, 
sei die Energieabteilung eines Chemiewerkes genannt. Diese hat im Gegensatz zu einer Uber- 
landzentrale nicht nur Strom, sondern auch Dampf verschiedener Druckstufen, Rohwasser, 
Reinwasser, Druckluft und sonstige ,,Energiearten‘‘ zu liefern, wahrend umgekehrt jeder 
einzelne der Energiebetriebe fast jede dieser Energiearten verbraucht oder liefert. 


(30a) 


4,.Zusammenfassung. Die Anwendung der Matrizenrechnung auf betriebswirtschaftliche 
Aufgaben, die bereits in einer friitheren Mitteilung an einfacheren Beispielen gezeigt wurde, 
fithrt auch bei schwierigeren Aufgaben zum Ziel. Der schmiegsame Kalkiil gestattet eine 
Darstellung der Lésung in sehr allgemeinen Fallen. Dabei zeigt es sich, daB es betriebswirt- 
schaftliche Fragestellungen gibt, die zu ahnlichen rechnerischen Schwierigkeiten fiihren wie 
die Eigenwertprobleme. Es erscheint daher zweckmaBig, auch bei der Bearbeitung betriebs- 
wirtschaftlicher Aufgaben ausreichende mathematische Hilfsmittel anzuwenden, um zu Lé- 
-sungen zu kommen, die praktisch befriedigen. 


| (Eingegangen am 28. September 1952.) 
- Anschrift des Verfassers: Dr.-Ing. O. Pichler, Leuna (Krs. Merseburg), Wickenweg 9. 
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Die Berechnung des Profilwiderstandes aus der vorgegebenen Profilform. 


Von E. Truckenbrodt. 


1. Einleitung. Die Berechnung des Profilwiderstandes umstrémter Kérper ist mit den 
Methoden der Grenzschichttheorie méglich, solange die Reibungsschicht am Kérper noch 
nicht abgelést ist. Erstmalig wurden fiir den ebenen Fall (zweidimensionale Strémung) Ver- 
fahren von J. Pretsch! und H. B. Squire und A. D. Young? angegeben. Ausgehend von der 
potentialtheoretischen Geschwindigkeitsverteilung um den Kérper ist in beiden Fallen eine 
Grenzschichtrechnung auszufithren, deren Ergebnis (Impulsverlustdicke an der Profilhinter- 
kante) dann in geeigneter Form in die Widerstandsformeln eingeht. Mit der Vervollkommnung 
und Verbesserung der Verfahren zur Berechnung der Reibungsschicht lieBen sich auch die 
Widerstandsformeln verbessern und vereinfachen. H. B. Helmbold® gelangt so zu einer ge- 
schlossenen Formel zur Berechnung des Profilwiderstandes, in der nur die Widerstands- 
beiwerte der ebenen Platte und die potentialtheoretische Geschwindigkeitsverteilung um den 
Kérper vorkommen. Durch geeignete Abstimmung der vorkommenden Konstanten beseitigt 
er die an der Hinterkante zuweilen auftretenden Schwierigkeiten. Erwahnt sei auch ein Ver- 
fahren von N. Scholz itber die Berechnung des Strémungswiderstandes schlanker Kérper mit |) 
beliebig rauher Oberflache. 

Zur Berechnung der an der Hinterkante benétigten Impulsverlustdicke wurde inzwischen 
von E. Truckenbrodt® unter Beriicksichtigung der neuesten Forschungsergebnisse insbesonders 
bei turbulenten Strémungen eine einfache Quadraturformel hergeleitet. Die potential- 
theoretische Geschwindigkeitsverteilung um den Kérper 148t sich fiir schlanke Profilkérper 
nach der Singularitaétenmethode von F. Riegels® in einfacher Weise berechnen. 

Wir wollen uns die Aufgabe stellen, die Widerstandsformel so umzuformen, das darin 
nicht mehr die Geschwindigkeitsverteilung sondern unmittelbar die Profilform vorkommt. 
Eine solche Umformung ist méglich fiir die zirkulationsfreie Umstrémung eines schlanken 
Profilkérpers. Wir geben also eine Formel zur Berechnung des Widerstandsanteils an, der 
durch die Verdrangungswirkung des Profiltropfens entsteht. Den Widerstandsanteil einer 
Sklettlinienwélbung schatzen wir ab. 


2. Die Ausgangsformeln. Den Widerstand eines mit der Geschwindigkeit U., angestrémten 
Kérpers der Breite 6 erhalt man in bekannter Weise aus der Geschwindigkeitsverteilung im 
Nachlauf in groBem Abstand hinter dem Kérper zu 


2) co 
W=c, F->U;, = ob [ U. (U,, — ux) dy. (ij 
y=— 2 
Dabei ist F = 61 die Projektionsflache des Korpers, J seine Lange und u,, (y) die Geschwindig- 
keitsverteilung im Nachlauf weit hinter dem Fliigel. 
Unter EKinfithrung der Impulsverlustdicke in unendlichem Abstand hinter dem Kérper 


Pee Oe @) 


]aBt sich (1) in der Form 
Ce (la) 
schreiben, wo c, den Widerstandsbeiwert bedeutet. 


1 J. Pretsch, Jahrb. 1938 d. Deutsch. Luftfahrtforsch., S. I, 61. 
oe B. Squire und A. D. Young, The calculation of the profile drag of aerofoils, ARC-Rep. 1838 

° HH. B. Helmbold, Ing.-Arch. 17 (1949), S. 273. Ein ahnliches Verfahren wird angegeben von 
E. M. Dowlen, Journal Roy. Aeron; Soc., 56 (1952), S. 109. 

* N. Scholz, Jahrb. d. Schiffbautechn. Gesellsch. 45 (1951), S. 244. 

5 E. Truckenbrodt, Ing.-Arch. 20 (1952), S. 211. 

° F. Riegels, Ing.-Arch. 16 (1948), S. 373; 17 (1949), S. 94 und 18 (1950), S. 329. Vel. auch E, Trucken- 
brodt, Ing.-Arch. 18 (1950), S. 324 und 19 (1951), S. 365. 
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Der wesentliche Inhalt der Verfahren von Pretsch und Squire- Young besteht nun darin 
eine Beziehung zwischen der Impulsverlustdicke an der Profilhinterkante und der imgpale: 
verlustdicke in unendlichem Abstand hinter dem Profil herzustellen. 

Fir die Impulsverlustdicke an der Hinterkante x = x, hervorgerufen durch die Reibungs- 
)wirkung einer Profilseite (Ober- bzw. Unterseite) schreiben wir 


cea Ba | Uh 
dh | = (1 a iy (3) 
=U) 


. Nach Squire- Young ergibt sich auf halbempirischem Wege fiir das Impulsverlustdicken- 
-verhaltnis die Beziehung 
Hy+5 

Oe Un\ 2 
' w=? (a) (4) 
H;, bedeutet den Grenzschichtformparameter (Verhaltnis von Verdrangungsdicke zu Impuls- 
-verlustdicke) an der Hinterkante. Da das Geschwindigkeitsverhaltnis U,/U,, an der Hinter- 
-kante nur wenig vom Wert eins abweicht!, AA 
ispielt die GréBe (H,+5)/2 keine sehr DH, 
sgroBe Rolle. Wahrend 4H. Schlichting? y 
eden Wert 3,2 vorschlagt, rechnet H. B. 
| Helmbold mit 3,5. 
In Abb. 1 haben wir das Impulsverlust- 
dickenverhaltnis }_/@, fiir die Formpara- 
jmeter H,—1.4 und 2,0 in Abhangigkeit 
-vom Geschwindigkeitsverhaltnis U,/U,, 
-aufgetragen. 
| Nach Pretsch gilt fiir das Impulsverlust- 
-dickenverhaltnis die Beziehung 


Bo 


'wobei sich die GroéBe € ein fiir allemal in 
| Abhangigkeit vom Geschwindigkeitsver- 
|haltnis U,/U., und vom Gruschwitzschen 
|Formparameter y darstellen lat. Wir 
fentnehmen diesen Zusammenhang einer 
| Arbeit von A. Walz? und habenihninAbb. 1 
wiedergegeben. Die Ubereinstimmung mit 
‘den Kurven von Squire- Young kann als 
gut bezeichnet werden. Bei unseren Rech- 
nungen wollen wir fiir das Impulsverlust- 
dickenverhaltnis setzen 


nach Fretsch 
——— nach Squire-Young 


a U, \3 

| Pu 9( Ua)? on 

| dD, tie We 

5 ee: mage | | 10 
Auch dieser Zusammenhang ist in Abb. 1 ~g9 088 0.50 0% 094 0.96 0.98 700 
dargestellt. Abb. 1. Das Impulsverlustdickenverh4ltnis Joo /#, in Abhingigkeit vom 


Die Impulsverlustdicke an der Hinter- Geschwindigkeitsverhiltnis Uj;/Uco nach J. Pretsch und H. B. Squire 
. . . . 

Kante einer Profilseite berechnen wir nach und A. D. Young 

‘der Quadraturformel von E. Truckenbrodt*. Diese lautet: 


%u/1 ae xp|! 1 
On 1 (Gls 3 1+n U \5 x 2 ltn U \3+2n x l-+n 
satay (Payer) sa [ao 
0 x,[l 


1 Diese Aussage gilt fiir die wirklichen Geschwindigkeiten, wahrend sich fiir die Potentialgeschwindig- 
keiten bei endlichem Hinterkantenwinkel an der Hinterkante der Wert null ergibt. 

2 H. Schlichting, Grenzschicht-Theorie, 8. 467. Karlsruhe (1951). ; 

8 4, Walz, Einige Beispiele zur theoretischen Berechnung der Polare eines Tragfligelprofils. 
Forsch. Ber. d. Deutsch. Luftfahrtforsch. Nr. 1848 (1943). 

4 Siehe FufSnote 5 auf Seite 176. 
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Hierin ist x der vom Staupunkt aus gemessene Abstand auf der Profilkontur. Er ist bei 
schlanken Profilen, mit denen wir uns im folgenden befassen wollen, in guter Naherung gleich 
dem Abstand gemessen auf der Profilsehne. Es gilt dann 4%, =I. Die GréBe x, gibt die Lage 
des Umschlagspunktes vom Jaminaren in den turbulenten Strémungszustand an und wird im 
folgenden als bekannt angenommen. U/U,, ist die potentialtheoretische Geschwindigkeits- 
verteilung um den Profilkérper. Es bedeuten cy, und cy, die Widerstandsbeiwerte einer mit der 
Geschwindigkeit U,, langsangestrémten einseitig benetzten ebenen Platte von der Lange I. 
Fiir glatte Profiloberflachen sind in Abb.2 die Widerstandsgesetze der ebenen Platte dar- 
gestellt. Ist die Oberflache rauh, so setzt man fiir cy die entsprechenden rauhen Wider- 
standsbeiwerte der ebenen Platte ein. Die Gré&e n wird von E. Truckenbrodt mit 
1 
Te . (Fj 


angegeben. 


or: 10° 


795 70 10 ee 10° 
Abb. 2. Die Widerstandsgesetze der liangsangestrémten ebenen glatten Platte. 
1,328 
(1) Blasius = laminar 
VRe 
7 
(2) Prandtl ee a8 2 
VRe 
(3) Prandtl- 0,455 oes 
aed Ope SSS t 
Schlichting f (le Re) 2,58 urbulent 
0,03 
(4) Falkner ON oe 
VRe 


Die Formel fiir die Impulsverlustdicke #, kénnen wir angenahert auch so schreiben, dab 
wir jeden Term fiir sich mit 1/(1--n) potenzieren’. Physikalisch bedeutet das, daf hinter dem 


Umschlagspunkt die turbulente Reibungsschicht sich so ausbildet, als wenn keine laminare 


Vorstrecke vorhanden ware. Diese Annahme ist bei gréBeren Reynoldsschen Zahlen gut 
erfiillt. Man vergleiche hierzu die noch spater gemachten Angaben fiir den Widerstand der 


langsangestrémten ebenen Platte bei verschiedenen Lagen des Umschlagspunktes. 


' 
| 


Wir setzen jetzt samtliche Beziehungen in die Widerstandsformel (la) ein und finden 


schlieBlich 
; ay 1/2 1 6/7 
ere. U \5 U \10/3 
7 | (Be) dx| +e, | ta dx||_. (9) 
0 %y 
Der vereinfachten Schreibweise wegen wollen wir fortan schreiben 
au 
1 — xX ° 


1 Mathematisch gesehen setzen wir also (x-+y)¢~xa-+ ya, Der maximale Fehler, den man bei 
einer solchen Entwicklung macht, ergibt sich bei x=y. Er betragt o TASS oRa he SOMA i0 
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Wie schon in der Einleitung erwahnt wurde, wollen wir die Geschwindigkeitsverteilung 
| U/U., nach dem Singularitatenverfahren von F. Riegels! berechnen. Fiir den hier betrachteten 
Fall des in Sehnenrichtung angestrémten symmetrischen Profilkérpers gilt 

il 
U 1 1 dy dx’ 

cy ee | 4 

P= 0 . 
‘Hierin stellen y = y/I die Ordinaten des Profiltropfens dar. Der Faktor 1//1-+-y? bildet eine 
Korrektur fiir die Stellen starker Kriimmung der Profiloberflaiche (Nase). Da wir uns mit 
schlanken Kérpern befassen wollen, wollen wir ihn gleich eins setzen. Fiir unsere weiteren 
Rechnungen ist es zweckmaBig, die Formel (10) in etwas anderer Form zu verwenden. Aus 
der Theorie der Quell-Senkenbelegungen erhalt man 


1 
U Af . , 2 ( —x’)? a 
| ga ()=1 +4 tim Coe een dx’. (11) 


20 


‘Durch Bildung des Grenzwertes und entsprechende partielle Integration kann man die Iden- 
titat mit (10) zeigen. 

3. Die Widerstandsformel. In (9) kommt die Geschwindigkeitsverteilung mit der Potenz 
+3-+2n (n=1 fiir das erste Glied und n=1/, fiir das zweite Glied) vor. Unter der Annahme 
schlanker Profile, d.h. kleiner Verdrangungsgeschwindigkeiten im Vergleich zur Anstrémungs- 
geschwindigkeit, U~U., kénnen wir in linearisierter Form schreiben 

1 
y \3t2" Fares bag, 
(Zia) =1 4842 Mim | y(a) eo eae x (12) 


wy oa 


“Mit diesem Ansatz gehen wir in (9) hinein und finden nach Ausfiithrung der Integration iiber x 
und nachfolgender Grenzwertbildung sowie weiterer Linearisierung 


1 


1 
Co 5 (x) 1 20 y (x) 6/7 
as cp (1 + 55 53 eer dx) x? + ey (1 = ae Cae ax) (L—-~*,,)a". (13) 
0 0 


Xu 


Die Integrale sind im Sinne des Cauchyschen Hauptwertes zu nehmen. Fiir die Integrations- 
-veranderliche haben wir x statt x’ geschrieben, da Verwechslungen ausgeschlossen sind. 
Wie wir aus (13) sofort ersehen, andert sich der Widerstand linear mit der Profildicke. 

Fir das Profil verschwindender Dicke (y =0, Gp 
ebene Platte) erhalt man aus (13) sofort 


40 


Cw 1/2 6/7 
mirieh x)! “i e7(1—-*%,) ae (14) 


Zum Vergleich geben wir noch die Widerstandsformel 


aes 


08 


der ebenen Platte an, die man gewinnt, wenn man 
‘die oben gemachte Linearisierung nicht vornimmt. g¢ | 
N 

ae 7/6 7/12 7/6 va ae 

Tala? tef%a—a)] - (15) ) 
In Abb.3 haben wir die Kurven Cyp/2C,, nach den For- ey 
meln (14) und (15) fiir die beidon ReynoldsschenZahlen , , ee 
Re =108 (cy, = 0,0045) und Re — 10®(cz = 0,0021) p 
dargestellt. Wenn man beachtet, daB man die Lage Ly 
des Umschlagspunktes nicht sehr genau kennt, so 02 04 06 08 10 


kommt den Unterschieden der beiden Formeln bei Abb.3. Die ‘Waderstandsbeiwerte der langsangestrémten 
den technisch vorkommenden Reynoldsschen Zahlen — cbenen Platte in oe von der Umschlags- 
(Flugzeuge Re ~ 2.10’) keine groBe Bedeutung zu. 

Die Widerstandsdifferenzen nach beiden Formeln sind nicht gréfer als diejenigen infolge einer 
geringen Verschiebung des Umschlagspunktes, und zwar Ax, =0,08 bei Re =10° und 4x, =0,02 
bei Re =10%. 


1 Siehe Fufnote 6 auf Seite 176. 
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Bei voll-turbulenter Strémung ergibt sich aus (13) mit x, =0 


1 
Cw as 20 y (x) 
ue en [1 4 EA da. (169 
6 
Wir wollen jetzt (13), ahnlich wie es E. Truckenbrodt schon bei der Berechnung der poten- 
tialtheoretischen Profileigenschaften getan hat, in eine Summenformel umformen, in der auBer 
ein fiir allemal berechneten Konstanten nur noch die Profilordinaten an bestimmten Stellen 


der Sehne voikommen. 


Nach dem Vorschlag von F. Riegels? beschreiben wir den Profiltropfen durch den Reihen- 


ansatz 


=> b, sin vy. (17) 
T 
Zwischen ~ und x besteht dabei der Zusammenhang 
x= (1+c08 Q) . (18) 


Fiihren wir den Ausdruck 
1 


BSN Fen) nt dy (19) 
0 


A Xu 
0 


ein, dann kénnen wir fiir (13) auch schreiben 


Cw » dt Jers /2 20 J(xu=1) — J (xu 6/7 
Z= H|1+5 = ae | wt ! ab sep = a hz)". 20) 
Setzen wir (17) in (19) ein, dann folgt nach einer bekannten Formel von H.Glauert? 


1 
Pepe ae COSY Du - (21) 


Nehmen wir jetzt zur Beschreibung des Profiltropfens nach (17) eine endliche Reihe y< N an, 
dann sind die Koeffizienten durch die Besselschen Formeln (vgl. Hiitte 27. Auf]. 5.208) folgen- 
dermaSen bestimmt : 

4 N= 


b=5 > 


=1 


sin aie Wonks (22) 


Es bedeuten yp, = (y(/L)(%n) die Ordinatenwerte des Profiltropfens an den Stellen 


=F (I | cos ae (23) | 


N 
Die Lage des Umschlagspunktes ist gegeben durch 
ox il mun 
Cy — m=z (1-008 NY: (24) 


wobei n= N fir x, =0 und n=0 fiir x, =1 gilt. 
Fiir die Summe (21) wird dann 


2 See ae (25) 
m=1 Ss NV eg Seen 


Bei der Summenbildung ist das Glied m=n auszulassen. 
Den Ausdruck (25) setzen wir in (20) ein und kénnen schreiben 


Cw = — 
Fae Sy Poe 239) +e( Qa + > eon2ym) (26) 
m= m=1 


1 E. Truckenbrodt, Ing.-Arch. 18 (1950), S. 324. 
2 Siehe FuBnote 6 auf Seite 176. 
° H. Glauert, Tragfliigel- und Luftschraubentheorie, S. 82. Berlin (1929). 
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fierin gelten fiir die Konstanten die folgenden Formeln: 
1 nm n\]1/2 1 6/7 
= ie (1 + cos al 4 Qu= E ( ces al : | 
/ 5 9 | Ie am [a= eee n 
= (1) orig ( (c a (m-+-n) ba inn) 
a a reer 2N 3 Pare Wie ak eh «oes? 
N 
10 2 1/7 ( ae tm t= (yr mt (m-+ ) a 
= a) Jetg ctg Mg OR "| 
| “or a] OS oN 2 aes Brew Ge: 


0,, —0. 


. Konstanten sind fiir N=12 fiir verschiedene Lagen des Umschlagspunktes berechnet 
vorden und sind in Tabelle 1 mitgeteilt. 


Tabelle 1. Die Konstanten in der Widerstandsformel (26) fiir N = 12 


xy 0 0.0670 0,25 0,50 0,75 1,00 
A 12 10 8 6 4 0 
ee 
P, Qn P, Qn P, Qn P, Qn P, Qn Ee O- 
0 1,0000 0,2588 | 0,9423 0.5000 | 0,7814 | 0,7071 | 0,5521 0,8660 | 0,3047 1,0000 0 
a RE i eee 
4g *m Pinn Qmn Pinn Qn Pinn Qmn Pin Omn Pinn Qmn Pinn Gan 
1 0,9830 0 1,8399 | —0,0078 | 1,8606 | —0,0187 | 1,9280 | —0,0402 | 2,0675 |—0,1020 | 2,3655 1,6099 0 
3 0,8536 0 0,6734 | —0,0284 | 0,6887 | —0,0715 | 0,7442 | —0,1726 | 0,8977 | —0,7217 | 1,6914 0,5893 0 
5 0,6294 0 0,4930 | —0,0736 | 0,5199 | —0,2107 | 0,6390 | —0,8735 | 1,3239 1,1480 |—0,7840 0,4314 0 
¢4 0,3706 0 0,4930 | —0,2314 | 0,5671 | —1,1257 | 1,1838 1,4835 |—0,7795 0.6197 |—0,1467 0,4314 0 
9 0,1465 0 0,6734 | —1.6383 | 1,1696 2,4285 |—0,7442 1,0059 |—0,1540 0,7217 |—0,0498 0,5893 0 
ial 0,0170 0 1,8399 8,1995 |—0,5914 3.3963 |—0,1051 2,3168 |—0,0358 1,8697 |—0,0130 1,6099 0 


Die Berechnung des Profilwiderstandes geht nun folgenderma8en vor sich: Gegeben sind 
lie Reynoldssche Zah] Re = U,l/y, die Lage des Umschlagspunktes x, = x,/l und die Profil- 
contur y =y(/1. Aus der Reynoldsschen Zah] und gegebenenfalls der Rauhigkeit der Profil- 
berflache bestimmt man die Widerstandsbeiwerte der mit der Geschwindigkeit U,, langs- 
ngestrémten ebenen Platte von der Lange / bei voll-laminarer und voll-turbulenter Strémung, 
f, und cy. Ist die Profilkontur durch einen einfachen analytischen Ausdruck gegeben, so 


‘ann man versuchen, den Widerstandsbeiwert nach (13) analytisch zu berechnen. Fiir andere 
-rofilformen wendet man die Summenformel (26) an. Man berechnet durch Bildung der 
-roduktsummen die Widerstandsbeiwerte fiir einige Umschlagspunktlagen. Das Ergebnis 
ragt man iiber x, auf und ermittelt den gesuchten Widerstandsbeiwert bei der gegebenen 
Jmschlagspunktlage. 

4. Beispiele. a) Joukowski-Profil. Die linearisierte Form des symmetrischen Jou- 


cowski-Profils lautet: 


y= 2e(1— x) /x(1— x) = = sin y (l—cos 9) . (28) 


fwischen der gréften Dicke t =t/l, die sich an der Stelle x,=0,25 befindet, und der GréBe « 
vesteht der Zusammenhang 


4, 
be == j= VUMODBE - 29 
3/2 (29) 
Jie Fourierschen Koeffizienten in (17) bestimmen wir sofort zu 
=e vadh (30) 
Weiter ergibt sich nach (21) unter Beachtung von (18) 
(a= ia 3 x, 4 2 a1). (31) 
Wir setzen sdmtliche Beziehungen in (20) ein und erhalten 
w 1/2 6/7 
= eg [1+ 1,924 (3 —2,) t] mu + &y, esate TO ONi == xa) Uae) (32) 


13 
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Fiir ein 15% dickes Profil (t =0,15) haben wir in Abb. 4 die laminaren und turbulenten Wide 
standsanteile in der Form (Ac,/2cf); und (Ac, /2cs), tiber der Lage des Umschlagspunktes 
dargestellt. Das Ergebnis entspricht der Lésung der Gleichungen (13) und (26). Zum Vergleic 
haben wir auch die Kurven eingetragen, die sich nach (9) durch graphische Integrationen d 
potenzierten Geschwindigkeitsverteilung iiber der Bogenlange ergeben. Die Ubereinstimmun 
mit den aus den linearisierten Formeln gefundenen Kurven ist gut. In Abb.5 haben wir au 
gehend von diesen Ergebnissen die Kurven ¢,/2c,, ther der Lage des Umschlagspunktes fii 
die beiden Reynoldsschen Zahlen Re =10° und Re = 108 aufgetragen. 

b) Ellipsen-Profil. Die Form des Ellipsen-Profils ist gegeben durch 


y=tyx(l—x) = + sing. (3 

Mit b,=t und J=~ (2x,—1) erhalt man nach (20) 
w : 6 

= = ef (1+ 2,500 t) xy!” + ef (1+2,857 t) (1—,)""" . (3 


I! Interessant an dieser Formel ist, dai sie i 


Yoo bezug auf die Abhangigkeit von der Lage de 
Umschlagspunktes derjenigen der ebenen Platt 
b-1 


(14) entspricht. Die Anwendbarkeit der For 


mel (34), die ja nur fiir den Fall anliegende 
el Strémung gilt, ist dadurch stark eingeschrank 
daB bei dem vorgegebenen elliptischen Prof 
Ln 
2C 
40 St 
BA 
\ 
LOSS 
. ib 
05 es NN \ 
\ 6 
\fe=105, 6, = 00045 
40 
: \ N 
0 az. o4 06 08 40 \ 
tz (Asm 
acr J ?\ 20 58 
45 Fb ?\2F/t 
uf he=108 : 
j ; e=10 ; \ 
a laminarer Anteil oy = 0.0021 ~ 
\ 
O4 
nach (13) 
We turbulenter Anteil ——— nach (9) 
G2 
Ly Ly 
0 G2 G4 06 G8 40 0 G2 G4 G6 48 40 
Abb. 4. Joukowski-Profil, t = 0,15, Geschwindigkeits- Abb. 5. Joukowski-Profil, t = 0,15, Widerstandsbeiwerte 
verteilung und Widerstandsanteile in Abhingigkeit von in Abhiangigkeit von der Umschlagspunktlage bei ver- 
der Umschlagspunktlage. schiedenen Reynoldsschen Zahlen. 


im hinteren Teil des Profils gewéhnlich bereits abgeléste Strémung herrscht. Die Widerstand 
werden also bei schon relativ kleinen Profildicken erheblich gréGer sein als die Werte, di 
man nach (34) errechnet. 


c) Abgeandertes Ellipsen-Profil. Durch den Ansatz 


noe yx) 
eee rie: (3 
mit = / —*% und b= esti set (36 


t 
kénnen wir das eben behandelte Ellipsen-Profil so abandern, da® wir Profile mit verdnder 


licher Dickenriicklage x, erzeugen kénnen. In Abb.6 haben wir einige Profile dargestellt 
Diese Profile sind im Bereich 0,191 <x,<0,809 wendepunktfrei. 
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Die Berechnung des Widerstandsbeiwertes wollen wir fiir den Fall vollturbulenter Stré- 
mung durchfiihren, indem wir (35) in (16) einsetzen. Zweckmafig fithrt man dabei zur ein- 
fachen Ausfiithrung der Integration Polarkoordinaten nach (18) ein. Man erhalt nach kleiner 
Zwischenrechnung 


+ = cy (1+ 2,857 t) . (37) 

Es ergibt sich das bemerkenswerte Ergebnis, da® bei den abgednderten Ellipsenprofilen 
der Widerstandsbeiwert bei voll-turbulenter Strémung nur von der Profildicke t, dagegen 
aber nicht von der Dickenriicklage x, abhangt. 
Auch bei diesen Profilen wird schon bei kleineren 
Profildicken im Profilhinterteil Ablésung eintreten. 

d) Abgeandertes Parabel-Profil. Von 
E. Truckenbrodt! wurde bereits an anderer Stelle 


7 


, 1 4% 
t 26, 

30|-_2 J | 

Gia tk Se es 
. S , 

N 

! aN 5 

6,07 a | al 

ee Ore ee oe) eee 
10 


abgeanadertes Ellipsenproti! 
— — — abgeinderies Parabefprotil 
° Joukowsk/-Protil 
OS a a | WACA=00., 

a) NACA -63-0.. 

° NACA -65-0.. 


o = NAUA-66-0.. 
GS 
ie ut 
0 A OY 06 Ge 40 0 GI G2 G3 O4 05 Oo G7 
Abb. 6. Ordinaten verschiedener Profilformen. Abb. 7. Die Widerstandsinderungen infolge der Verdrangungs- 
a) Abgeanderte Ellipsen-Profile, b) Abgeanderte Parabel- wirkung des Profiitropfens bei voll-turbulenter Strémung. 


Profile, c) NACA-Profile 


ein symmetrischer Profiltropfen mit verdnderlicher Dickenriicklage aus dem bekannten Parabel- 
zweieck in folgender Form angegeben: 


a, BC) 
Sarre ee) 


E 1—2 
rte es tl d= ae 
2 x} se 


(39) 

t 

Die Profile sind im ganzen Bereich 0<x,<1 wendepunktfrei. An der Vorder- und Hinter- 

kante laufen die Profile jeweils in eine spitze Kante aus. In Abb. 6 sind einige Profile dargestellt. 
Setzen wie (38) in (16) ein, so 1aBt sich die Integration sehr einfach ausfiihren. Fiir den 

Widerstandsbeiwert bei vollturbulenter Strémung erhalt man dann 


Cw 0,909 1—x, 
ee et ; 40 
‘ au f | 1—2n, ne 2 iP iN 
Fir das Parabelzweieck, x,=0,5, folgt hieraus: 
> = cy (1 +1,819 12). (41) 


1 Siehe FuBnote 1 auf Seite 180. 
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In Abb.7 haben wir die Widerstandsinderung durch den eben beschiriebenen Profiltropfe 
gegeniiber der langsangestrémten ebenen Platte in der Form 1/t- Ac,/2c, tiber der Lage de 
eréBten Dicke x, aufgetragen. 

e) NACA-Profile. Fiir die in Tabelle 2 angegebenen NACA-Profile wurden die Wider 


standsbeiwerte mittels der Summenforme] (26) ermittelt. 


Die Profilordinaten wurden dem NACA-Report Nr.824 entnommen und sind in Abb, 
dargestellt. 


bh Cp 10? Tabelle2. Die Widerstandsénderungen 
10 i bet NACA- Profilen. 
1 4%w 
Profil ore Dep, 
NACA—00. . 2,437 
G3—0R ae 2,051 
65—0 . . 2,063 
66—0.. Zales 


O4 


ve Fiir die Profile NACA 00, . ist eine vollstandige 

Rechnung fiir alle Umschlagspunktlagen durchgefihrt 

Us 7 worden. Das Ergebnis ist fiir die Reynoldssche Zahl 
a Re =6.10%in Abb. 8 iiber der Profildicke aufgetragen. 
[ees t Zum Vergleich ist auch die dem NACA-Report 

C Le G10 a G20 Nr. 824 entnommene Kurve der Messungen mit ein- 
song Dis Waacremarhivene de Poe NACS HD. getragen. Da die Umschlagepunktlagen fir de 
schlagspunktlagen. — — -- — Messung nach NACA- Messungen nicht bekannt sind, ist ein genauer Ver- 


Report Nr. 824. : . 
: gleich zwischen der Rechnung und den Messungen 


nicht moéglich. Da8® der Anstieg der MeBkurve stirker ist als derjenige der theoretischen 
Kurven bei konstant gehaltener Umschlagspunktlage diirfte damit zusammenhangen, daft 
der Umschlagspunkt mit wachsender Profildicke nach vorn wandert. 

Fir die anderen Profile wurde nur die Widerstandsanderung bei voll-turbulenter Strémung 
berechnet. Die Ergebnisse sind in Tabelle 2 und Abb.7 wiedergegeben. In dieser Abbildung 
ist auch der konstante Wert fiir die abgeinderten Ellipsen-Profile eingezeichnet. Unsere 
Werte fiir die NACA-Profile und das Joukowski-Profil liegen zwischen dieser Kurve und 
derjenigen fiir die abgedinderten Parabel-Profile. Die mittlere prozentuale Widerstands- 
erhéhung durch die Verdrangungswirkung des Profiltropfens betragt fiir die neueren ameri- 
kanischen Profile etwa den doppelten Betrag des prozentualen Profildickenverh4ltnisses!. 


5. Der Einflu8 einer Profilwélbung. Im folgenden wollen wir noch den Einflu® einer 
Profilwélbung auf den Widerstandsbeiwert abschatzen. Wir machen unsere Uberlegungen 
am Beispiel der voll-turbulenten Stiémung, x,—0. Wir gehen von der Widerstandsformel (9) 
aus, die wir naherungsweise wie folgt schreiben?: 


=o {ete “ 


Hierin bedeuten Uy und U, die Geschwindigkeitsverteilungen auf der Ober- und Unterseite 
des Profils. Fiir die zirkulationsbehaftete Umstrémung des Profils kénnen wir schreiben 


U U U 
a, > (os) (ae), ( 
Ks sind U, die bereits behandelten Verdi dngungsgeschwindigkeiten des Profiltropfens und U, 


die Geschwindigkeiten, die durch die Skelettlinienwélbung entstehen. Das obere Vorzeichen 
gilt fiir die Profiloberseite und das untere entsprechend fiir die Profilunterseite. Wir setzen 


1 N. Scholz (vgl. FuBnote 4 auf Seite 176) findet auf Grund seiner Beispielrechnungen etwa den 2,4- 
fachen Betrag des prozentualen Dickenverhiltnisses fiir die prozentuale Widerstandserhéhung infolge 
des Profiltropfens: 


* Die Vereinfachung besteht darin, daB® wir in der Formel fiir die Impulsverlustdicke (7) n=@ 
gesetzt haben. 
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1(43) in (42) ein und erhalten die Widerstandsinderung infolge der Profilwélbung zu 


I 
: 1 ow U \2 
Ac, Tat 6 Ch | Gal fis dx 5 (44) 
F s 


) Fir nicht allzu dicke Profile kénnen wir den EinfluB der Verdrangungsgeschwindigkeit des 
| Profiltropfens vernachlassigen. Wir setzen also (U/U,),~1. Dann wird aus (44) 


1 
Acw a fh OES ; 
es (5 dx. | (45) 
é 


\Die Widerstandsainderung wichst, da die Geschwindigkeiten U, proportional der Walbungs- 
jhéhe sind, mit dem Quadrat der Wélbungshohe. 

| Ehe wir jetzt fiir einige spezielle Skelettlinienformen (45) auswerten, miissen wir zundchst 
noch den Anstellwinkel festlegen, bei dem wir die Widerstandsinderung berechnen wollen. 
Die zusdtzlichen Geschwindigkeiten, die durch eine Anstellung des Skelett-Profils (Profil 
pverschwindender Dicke) entstehen, ergeben sich aus der ersten Birnbaumschen Normal- 
verteilung zu 


J are 
ey aes : 
)Wiirde man hiermit in (45) hereingehen und die Integration ausfiihren, so erhalt man eine 
funendlich groBe Widerstandserhéhung. Es gibt aber einen Anstellwinkel, bei dem der Anteil 
nach (46) verschwindet, ~n—a,;. Diesen Anstellwinkel bezeichnet man als den des sto freien 
Hintritts. Fiir ihn ergibt sich an der Profilvorderkante (x =0) keine unendlich groBe Ge- 
schwindigkeit, wie sie sich sonst nach (46) errechnet. Der Zustand des stoBfreien Eintritts 
stellt also das Widerstandsminimum dar. Die Berechnung des Anstellwinkels des stoBfreien 
Kintritts bzw. des zugeordneten Auftriebsbeiwertes kann man in bequemer Weise mit den 
on E. Truckenbrodt! angegebenen Formeln durchfiihren. Fiir den Auftriebsbeiwert gilt z. B. 


1 
y(s) 
Ca SS 47 
: | Vx (1—x)° St 
0 


Fiir drei Beispiele wollen wir die Widerstandserhéhung infolge der SkelettlinienwéJbung 
bestimmen. 

' a) Die neueren NACA-Profile (6-Serie) (vgl. NACA-Report Nr.824) besitzen Skelett- 
linien, die aus einer vorgegebenen Geschwindigkeitsverteilung errechnet worden sind. Dabei 
wurde angenommen, da die Geschwindigkeiten im vorderen Teil des Profils (0< x< a) 
konstant verlaufen und im hinteren Teil (a<(x<1) linear bis auf den Wert null abnehmen. 


U (Gee 
0<e<e: (G).=aatg (48) 
) Oe te = 7s (l1— x). 


Einige hieraus sich ergebende Skelettlinienformen sind in Abb.9 dargestellt. Fiir die gréBte 
Walbungshéhe gilt etwa f= (0,055 —0,077)c,,. Wir setzen (48) in (45) ein und finden fir die 
Widerstandsanderung 


Acw * 1 12a, 2 (49) 
2c 


St 
In Abb.9 haben wir dieses Ergebnis iiber a aufgetragen. Bei konstanter Geschwindigkeits- 
verteilung lings der ganzen Profiltiefe (a =1) wird Ac,/2 cz, = 0,188 cq, Nimmt man den Auf- 
triebsbeiwert des stoBfreien Eintritts mit c,,—=0,4 an, dann ergibt sich eine Widerstands- 
erhéhung infolge der Wélbung von etwa 3%. Allgemein kann man sagen, daB die Profil- 
wolbung nur unwesentliche Widerstandserhéhungen mit sich bringt, ein Ergebnis, das auch 
durch die im NACA-Report Nr.824 mitgeteilten MeBergebnisse gut bestatigt wird, 


7ig (IEP VS Gr 


1 Siehe Fu@Bnote 1 auf Seite 180, 
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b) Das bekannte Parabel-Skelett 
y= a (1—2) (50)}) 


Tt 


hat entsprechend der zweiten Birnbaumschen Normalverteilung eine elliptische Geschwindig-| 
keitsverteilung 


(vei ae a) 


aie Hiermit ergibt sich die Widerstandserhéhung zu 
“ag 
910 - : Acw _ 2 92 (52)| 
ee | 
74, oF 4 Vf 


Dieses ist etwa derselbe Wert wie fiir die Skelettlinie | 
mit konstanter Geschwindigkeitsverteilung. Die 
Wolbungshohe ist jedoch beim Parabel-Skelett um 
45% gréBer als bei dem Skelett mit konstanter Ge- 
schwindigkeitsverteilung. 

c) Eine S-Schlag-Skelettlinie ergibt sich in 
bekannter Weise aus 


16 
JY = Cry 35 


x(1— x)(1—2x). (53) 


Da fiir diese Skelettform c,,—0 ist, haben wir die 


Skelettlinienordinaten auf den Nullmomentenbeiwert | 
Cm, bezogen. Die Geschwindigkeitsverteilung ergibt 
sich aus der dritten Birnbaumschen Normalvertei- 
lung zu 


U 16 a E | 
ee Cm (1—2x) Yx(I— 4) « (54) 
Die Widerstandsanderung findet man aus (45) zu 
Accueil oe eee : : 
ic, Sar aed (6% 


Nehmen wir einen Nullmomentenbeiwert vonc,,, =0,1 
an, dann betragt die Widerstandserhéhung durch 
den S-Schlag der Skelettlinie auch wieder nur 3%. 


6. Zusammenfassung. In die bekannten Verfahren 
zur Berechnung des Profilwiderstandes schlanker 
Kérper von Pretsch und Squire- Young werden die 
neuesten Ergebnisse der Reibungsschichtforschung — 
006 | | eingefiihrt. Die neue Formel wird so umgeformt, daB 

die darin vorkommende Geschwindigkeitsverteilung | 
, um den Kérper fiir den Fall der zirkulationsfreien — 
, QP ~Ss«U mstrémung durch die Ordinaten des Profiltropfens _ 
Abb. 9, Die Widerstandsinderungen infolge Profil ausgedriickt wird. Die Reibungsschichteinfliisse gehen — 
bore SE ate cto Stremung. Skelett-Linie inform der laminaren und turbulenten Widerstands- — 

c eport Nr, ° 

beiwerte der mit der Anstrémungsgeschwindigkeit 
langsangestrémten ebenen Platte von der Lange des Profils ein. Die Profiloberflache kann 
dabei glatt oder rauh sein. Es werden einfache Quadratur- und Summenformeln ge- 
wonnen, mit denen die Widerstandsbeiwerte bei verdnderlicher Umschlagspunktlage fiir eine 
beliebig vorgegebene Profiltropfenform leicht und schnell ermittelt werden kénnen. An Bei- 
spielen werden die Formeln erprobt. Die Widerstandserhéhung infolge der Verdrangungs- 
wirkung des Profiltropfens ist proportional der Profildicke. Sie betragt bei voll-turbulenter 
Strémung prozentual etwa den doppelten Betrag des prozentualen Profildickenverhiltnisses. 
Der Einflu8 einer Skelettlinienwélbung wird abgeschatzt. Er ist bei den ttblichen Wolbungs- 
héhen nur gering. 


° Parabelskelett 


010 + 


(Eingegangen am 15.Oktober 1952.) 
Anschrift des Verfassers: Priv.-Doz. Dr.-Ing. E. Truckenbrodt, Braunschweig, Kérnerstr. 12. 
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) Zur Systematik des plastischen Momenten- und Spannungsausgleichs 
in Stabwerken unter Biegung. 


Von H. Craemer. 


1. Grundlagen. Wir legen einen ideal-plastischen Baustoff mit dem bekannten Span- 
tungs-Dehnungs- Gesetz 


o=—e fir o<e=e,, | 
Es Cie) 
| Ol— One hur COR OO) 
tei Belastung und 
o = = (e— &) (1.2) 


/ei Entlastung zu Grunde; hierin ist o, und ¢, Spannung und Dehnung an der Streckgrenze 
tnd ¢, die Restdehnung bei vélliger Entlastung nach Erreichen der Streckgrenze. Die Gleich- 
jewichtsbedingungen werden wir am unverformten System ansetzen. Nachdem unter einem 
Biegemoment M am Querschnittsrand die Plastizierung eingesetzt hat, breitet sich bekannt- 
fich der Vorgang tiber den ganzen Querschnitt aus, bis das vollplastische Moment 


M, = 0,5 (1.3) 


reicht ist; dabei ist S die Summe der statischen Momente der Querschnittsteile beiderseits 
der Flachenhalbierenden um diese. Es tritt alsoecin Spannungsausgleichim Querschnitt 
in. 

Mehr Beachtung als dieser hat der Momentenausgleich gefunden; er tritt in statisch 
inbestimmten Balkentragwerken auf, falls S$ konstant ist und kommt dadurch zustande, daB 
‘lurch FlieBen von Querschnittsteilen der Verformungswiderstand dort abnimmt, so da eine 
euverteilung der Momente gegeniiber den elastostatischen Werten eintritt. Die Héchstlast 
? = P, des Tragwerks ist erreicht, wenn durch volle Plastizierung sich FlieBgelenke mit 
/bwechselndem Drehsinn in solcher Zahl gebildet haben, daB das Tragwerk oder Teile des- 
elben beweglich werden. 

In vielen Anwendungen des Stahlbaus und fast immer im Stahlbeton ist aber S keine Kon- 
itante. Der Vorgang ist dann an sich derselbe wie oben fiir den Momentenausgleich beschrie- 
pen, jedoch gleichen sich dabei nicht die Momente, sondern nur die Biegespannungen aus; 
man wird also auf den Oberbegriff des Spannungsausgleichsim Tragwerk gefihrt, der 
owohl auf konstantes wie variables S anwendbar ist. Bei statisch bestimmten Systemen tritt 
benfalls dieser Spannungsausgleich ein; er ist aber nicht von einer Umlagerung der Momente 
vegleitet, sondern gleichbedeutend mit dem Spannungsausgleich innerhalb der meistbean- 
spruchten Schnitte. 

_ Zur Darstellung ziehen wir spaterhin die folgenden dimensionslosen, als Ausnutzung be- 
Be cineten GréBen 


| o M ib E 
1 — me > i Shera M. ° QO = P; * QO = a (1.4) 
a Mit », und w, bezeichnen wir die‘ Ausnutzung des betreffenden Systems bei FlieS- 


beginn?. 

2. Besondere Arten des Spannungs- und Momentenausgleichs. Bei den oben bezeichneten 
lrei Arten des Ausgleichs (Spannung im Querschnitt, Momente und Spannungen im Tragwerk) 
<gnnen verschiedene Sonderfille auftreten. Man kann zunachst die Falle betrachten, wo der 
m Endzustand vorhandene Ausgleich bereits im elastischen Zustand besteht und wahrend des 
FlieBens erhalten bleibt. Dies sei als elastischer Ausgleich und der entgegengesetzte Fall 
us plastischer Ausgleich bezeichnet. 

Man kann ferner danach fragen, ob der Ausgleich sich auf samtliche Elemente des betrach- 
eten Systems oder nur auf Teile desselben erstreckt; im ersten Falle werden wir von kon- 


1 Hf. Craemer, Ing.-Arch. 20 (1952), S, 129. 
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tinuierlichem Ausgleich, andernfalls von értlichem Ausgleich sprechen. Unter Ele- 
menten sind dabei zu verstehen: beim Spannungsausgleich eines Querschnitts seine Flachen 
differentiale, beim Momentenausgleich die Querschnitte des Tragwerks, beim Spannungs- 
ausgleich im Tragwerk die Flachendifferentiale seiner simtlichen Querschnitte. 

Bei durchlaufenden Balken auf unbeschrankt tragfahigen Stiitzen vollzieht sich im ail 
gemeinen der Momenten- oder Spannungsausgleich nur in einem einzigen Felde, indem sich 
iiber oder in der Nahe von zwei Stiitzen (bei Endfeldern: einer Stiitze) und im Felde FlieB- 
gelenke bilden. Die Nachbarfelder bleiben dabei elastisch und erfahren nur Formanderungen, 
die gegeniiber den plastischen verschwinden; schlieBen also rechts oder links zwei oder mehr} 
Felder an, so bleiben diese statisch unbestimmt. Diese Art des Ausgleichs sei unvollstan- 
diger Ausgleich genannt. | 

Nur bei ganz spezieller Bemessung des Tragwerks und nur fiir einen bestimmten aus im} 
allgemeinen mehreren méglichen Belastungsfallen erfolgt der Ausgleich in samtlichen Feldern 
gleichzeitig und heiBe dann vollstandiger Ausgleich. Nur in diesem Sonderfalle ist diet 
wiederholt aufgestellte Behauptung! richtig, da fiir den Zusammenbruch des n-fach un-| 
bestimmten Systems n +1 FlieBgelenke nétig seien, so daB es vor Bildung des letzten Ge- 
lenkes statisch bestimmt wird. | 

Der unvollstandige Ausgleich ist auch bei durchlaufenden Rahmen, die aber hier wegen des 
Auftretens von Normalkraften nicht behandelt werden, und bei Tragerrosten? die Regel. In 
den spateren Beispielen werden die Begriffe des elastischen, kontinuierlichen usw. Ausgleichs} 
jeweils nur auf ein Feld eines durchlaufenden Balkens bezogen. 


3. Spannungsausgleich im Querschnitt. a) Elastischer Spannungsausgleich. Siehti 
man beim idealen I-Trager davon ab, daB die Spannungen im Flansch mit wachsendem Ab- 
stand vom Steg etwas abnehmen, was aber nur fiir sehr breite Flanschen Bedeutung hat, soi 
erkennt man, daB hier die Spannungsverteilung statisch bestimmt ist und sich infolge- 
dessen wahrend des ganzen Belastungsvorgangs nicht a4ndern kann. Es liegt also elastischer 
Ausgleich vor (Abb. 1), und es ist bei Be- oder Entlastung 


te Wihtiri 00S ph (3.1) 
Oe fr Oe aL 
lings/orgcos stir Gevi==1ita| 
Bei Entlastung nach Erreichen der Héchstlast ist 
C=O +, (3.3) 
wenn 0, die zu é) gehérige Ausnutzung bezeichnet. 
b) Plastischer Spannungsausgleich. Bei jeder anderen Querschnittsform hingeger 
kann das Moment nur durch Heranziehung der plastischen Eigenschaften voll entwickelt wer- 
den: plastischer Ausgleich. Aus dem 


Verformungsgesetz (1.1) bei Belastung 
folgt 


und 


bzw. 


i= on ons Ve (3.4) 


die Ausnutzungslinien fiir Spannung 


= 
7 


Abb. 1. Ausnutzung bei elastischem Spannungsausgleich 


cewoustonie. Abb. 2. Ausnutzung bei kontinuierlichem‘plastischem Spannungsausgleich 


im Querschnitt. 


und Dehnung stimmen also insoweit iiberein, f(u, v) =f (0, v). Weiter ist fiir simtliche Ele- 
mente #=1, o=o fir yl; (3.5) 
1 F. Bleich, Stahlhochbauten, 1. Bd., S. 406, Berlin 1932; J. F. Baker, Design of steel frames, Struc. 


Engr. Okt. 1952; J. A. van den Broek, Theory of limit design, S. 32, New York. 
® H. Craemer, Ing.-Arch, 14 (1943), S. 233. 
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d.h. alle Spannungsausnutzungslinien laufen im Punkt 4 =v =I1 zusammen und diejenigen 
fiir Dehnung haben eine Asymptote in »y = 1. Die Entlastungslinien stimmen mit denen fiir 
Belastung iiberein, solange » <v, ist; daritber hinaus sind sie den Asten mit 0 AW aE 
parallel (siehe Abb. 2). 

_ c)kontinuierlicher Spannungsausgleich. Mit (3.5) ist auch die Frage nach dem Auf- 
treten des kontinuierlichen Ausgleichs beantwortet: er tritt in allen Fallen auf. 


4, Momentenausgleich. a) Ortlicher plastischer Momentenausgleich. Dieser 
Fall kommt am haufigsten vor und ist bekannt. 

b) Ortlicher elastischer Momentenausgleich. Das bekannteste Beispiel hierfiir 
ist der voll eingespannte Balken mit konstantem Tragheitsmoment J und einer Einzellast P 
‘in der Mitte, also elastostatischen Gré8tmomenten M = +- PI/8. Die Ausnutzungsfunktion 
-lautet hier fiir die voll plastizierten Schnitte a 
| y=+to, (4.1) 

fiir alle anderen (b, c) hingegen 

Pe hoe. tot © = Vokes tT (4.2) 

(siehe Abb. 3). Bei Entlastung gelten die gleichen Bezichungen, es bleiben keine Restmomente 
im Tragwerk zuriick. Es gibt eine Reihe anderer Systeme mit elastischem Ausgleich, z. B. den 
| Zweifeldtrager, wo er kiinstlich durch Stiitzensenkung herbeigefiihrt wird. 
Beim Grtlichen, elastischen oder plastischen, Momentenausgleich 
_besteht die aus den plastischen Verformungen allein gewonnene Biege - 
linie aus geraden Stiicken mit unvermitteltem Richtungswechsel in 
den FlieBgelenken. Im Falle der Abb. 3 ist sie also ein Dreieck mit 
der Spitze in der Mitte. 

ce) Kontinuierlicher, plastischer Momentenausgleich. 
Dieser Fall ist sehr selten. Aus kinematischen Griinden miissen 
namlich Gebiete mit positiver und negativer Kriimmung abwechseln, [vI 
beim kontinuierlichen Ausgleich also unmittelbar aneinander an- RU ecu aR a8 
stoBen. Dies ist aber nur méglich bei Belastung durch auBere Mo- ae i 
mente, z. B. bei dem in Abb. 4 dargestellten beiderseits durchlaufen- 
den Balken, wo vollplastisch M,= — M, wird. Ist elastisuh M, <<|M,|, 
was vom Verhaltnis c/l und der Steifigkeit und Belastung der an- 
schlieBenden Felder abhangt, so wird die Ausnutzungslinie wie an- 
gegeben; dabei wurde idealer I-Trager angenommen, da nur dann 
die Ausnutzungslinien stiickweise gerade sind. 

d) Kontinuierlicher elastischer Momentenausgleich. 
Bei entsprechender Abstimmung von e/! auf das Verhalten der 
| Nachbarfelder, z. B. c/] = 0,5 und J = oo in den Anschluffeldern, 
_ wird schon elastisch M, = — M,, d.h. der Ausgleich ist kontinuier- 
lich elastisch. Die Ausnutzungslinien a und 6 gehen dann in eine 
einzige Gerade durch den Ursprung iiber, und zwar unabhangig von Fea hk EET 
der Querschnittsform. scher Momentenausgleich. 

Die Biegelinie besteht beim kontinuierlichen, plastischen oder 


elastischen, Ausgleich nur aus gekriimmten Stiicken mit wechseln- eye 
der Kriimmung, d. h. aus einer Gelenkkette. 7 
5. Spannungsausgleich in statisch bestimmten Tragwerken. 
Wie bekannt, gleichen sich in den meist beanspruchten Schnitten a fee 
des Tragwerks die Spannungen aus: értlicher plastischer Aus- bay 
gleich. 


Wahlt man die Widerstandsmomente W so, daf in allen Schnit- be 
-ten M/W =—konst., so wird das System traditionell, aber irre- ayy, 5. qusnutzung bei Byiliebers 
 fiihrend als Tr ager gleichen Widerstands bezeichnet, wabrend — clastischem dao dele ol im 
. ragwerk. 

-es richtiger Trager gleicher Spannung, genauer: elastostatisch ¢ sd 
bestimmter Randspannung, hei®en mite. Statt dessen kann man auch die Bedingung 

M 

Sg konst. (5,1) 
einfiihren, die bei wechselnder Querschnittsform nicht mit der vorgenannten identisch ist. 
- Unter der Hichstlast werden dann alle Schnitte gleichzeitig mit o, ausgenutzt. Wir haben dann 


190 Craemer: Zur Systematik des plastischen Momenten- u. Spannungsausgleichs usw. _Ingenieur-Archiy 


mit dem kontinuierlichen plastischen Spannungsausgleich im statisch be- 
stimmten Tragwerk zutun. Da P/M zufolge der statischen Bestimmtheit unveranderlich 
ist, beeinflussen sich die Querschnitte nicht gegenseitig und die Ausnutzungsfunktion f(y, @) 
stimmt mit derjenigen f(u, 7), z. B. in Abb. 2, ttberein. 

Wahlt man bei zunachst beliebigem Verlauf von S einen idealen I-Trager, so ist der Quer- 
schnitt selbst ebenfalls statisch bestimmt. Alle Linien f(u, @) sind dann Gerade durch den 
Ursprung. Falls in zwei oder mehr von ihnen volle Ausnutzung erreicht wird, haben wir mit 
dem értlichen elastischen Ausgleich zu tun. Die zugehérigen Ausnutzungslinien sind 
dann durch Abb. 5 gegeben; in den voll ausgenutzten Schnitten kénnen Restdehnungen, aber 
keine Restspannungen auftreten. 

Der Spannungsausgleich wird ferner kontinuierlich elastisch, wenn auferdem die 
Bedingung (5.1) erfiillt ist. In Abb. 5 fallen dabei alle Geraden mit /w = g/m <1 fort. 

6. Spannungsausgleich im statisch unbestimmten Tragwerk. a) Ortlicher plastischer 
Spannungsausgleich. Ist S verdnderlich, so wird unter der Héchstlast in den FlieB- 
gelenken M/S = o,, wobei M/S ein Extremum gegeniiber den Nachbarschnitten bildet ; die Un- 
tersuchung kann nach einer friitheren Arbeit! geschehen. Da in den tibrigen Querschnitten 
o <0, ist, haben wir nur drtlichen Ausgleich. 

b) Kontinuierlicher plastischer Spannungsausgleich. Man kann nun, ohne obige 
Extremalbedingung zu verletzen, die GréBen S in den Schnitten zwischen den FlieBgelenken — 
solange verringern, bis iberall M/S =o, ist, so da8B kontinuierlicher plastischer Ausgleich 
vorliegt. 

Aus kinematischen Griinden wechseln aber positive und negative Momente ab, so da8 bei 
einseitiger (beidseitiger) Kontinuitat eines Feldes ein (zwei) Momentennullpunkte auftreten. 
Beim kontinuierlichen Ausgleich mu® dann an diesen Stellen S = 0 sein. Dies ist aber nur 
moglich, wenn gleichzeitig das Tragheitsmoment J = 0 ist. Hierdurch nun wird das betr. 
Feld statisch bestimmt, d.h. ein Ausgleich vollzieht sich nur innerhalb der Querschnitte. 
Es bedarf wohl kaum des Hinweises, daB S = 0 mit Riicksicht auf die Schubspannungen und 
auf wechselnde Lastfalle praktisch nicht verwirklicht werden kann, sondern nur als lehrreicher 
theoretischer Grenzfall hier herangezogen wurde. 

Ist S unbekannt, d. h. soll das Tragwerk be messen werden, so kann bekanntlich in einem 
durchlaufenden Balken die SchluBlinie im Rahmen der kinematischen Bedingungen beliebig 
angenommen werden. Bemi®t man dann in der kinematisch notwendigen Zahl von Quer- 
schnitten nach S = M/o,, so erreicht man fiir den betreffenden Lastfall den vollstandigen 
plastischen Spannungsausgleich; tut man dasselbe in allen Schnitten, so ist der Ausgleich 
auBerdem kontinuierlich. 

ce) Kontinuierlicher elastischer Spannungsausgleich. Dieser besteht, wenn man 
einen idealen I-Trager als Querschnitt wahlt, da dann auch innerhalb der Querschnitte keine 
Neuverteilung méglich ist. 

d) Ortlicher elastischer Spannungsausgleich. Derselbe kann auftreten, wenn bei 
Verwendung idealer I-Trager die GréBen S auf Steifigkeit und Belastung des ganzen Systems 
entsprechend abgestimmt werden, vgl. auch das unter 4b Gesagte. 

7. Zusammenfassung der Ergebnisse. Der Spannungsausgleich im Querschnitt ist stets 
kontinuierlich und kann sowohl plastisch wie elastisch sein. Der Momentenausgleich im Trag- 
werk ist im Normalfalle drtlich plastisch, in Sonderfallen értlich elastisch; die an sich denk- 
baren Falle des kontinuierlichen, plastischen oder elastischen, Ausgleichs kommen praktisch 
kaum vor. Beim Spannungsausgleich im Tragwerk ist wieder der haufigste Fall der értliche 
plastische Ausgleich, wahrend der értliche elastische Ausgleichin Sonderfallen vorkommen kann. 
Die weiteren Falle des kontinuierlichen, plastischen oder elastischen, Ausgleichs im statisch 
unbestimmten Tragwerk lassen sich praktisch nur angenahert verwirklichen, sind aber als 
Grenzfalle wichtig. 

Die Ausnutzungslinien f(u, v), f(y, w) und f(u, w) bestehen beim elastischen Ausgleich aus 
Geraden durch den Ursprung und es ist », baw. w, = 1. Beim kontinuierlichen Ausgleich gehen 
unter der Héchstlast alle Ausnutzungslinien durch einen Punkt. 


(Eingegangen am 27.Oktober 1952.) 
Anschrift des Verfassers: Professor Dr. H. Craemer, Alexandria (Agypten), Bulkeley, Immeuble Binorr. 


1 H. Craemer, Applications following from the maximum conditi , Abh, Int. Ver. Brii : : 
bau 12. Bd., 1952. ao nt. Ver. Briicken- u. Hoch 
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Zur Berechnung von Wirbelverteilung und Auftrieb eines diinnen 
Unterschallprofils in zwei hintereinander angeordneten Fliigelgittern 
bei kompressiblen Strémungen. 


Von Paul F. Byrd und Mary T. Huggins. 


_ 1. Einleitung. Vermittels einer Methode, die auf einem einfachen Singularitétenverfahren 
mur Erzeugung von Filiigelgittern und auf der dabei erhaltenen Integralgleichung beruht, 
berechneten J.Dérr1 und K. Nickel? in ihren kiirzlich hier erschienenen Aufsitzen Wirbel- 
verteilung und Auftrieb fiir ein Profil im senkrechten Gitterverband. Man interessiert sich 
heute jedoch auch fiir die Potentialstrémung durch mehr als ein Schaufelgitter®. Im folgenden 
wird daher gezeigt, da} die Arbeiten von Dérr und Nickel sich ohne Schwierigkeit auf den 
Fall von zwei hintereinander angeordneten Gittern ausdehnen lassen. (Das Verfahren labt 
sich auch auf beliebig viele (m) Fligelgitter erweitern. ) 


2. Die Integralgleichung. Betrachtet man zwei hintereinander liegende Fligelgitter, so 
lautet die (singulare) Integralgleichung* zwischen den Wirbelbelegungen J\(x), [(x) und der 
induzierten Abwartsgeschwindigkeit W(x) der Skelettlinie der Fliigel 


W(x) = — 55 { Hla) Gta ee —m)dx, — 54 [ T(x) Ctag(s —a)du (0 <a<b<e) (1) 


b 
it 
‘a pie) 
g=%. p=li—M. (2) 
obei M die Machzahl der Abstrémungsgeschwindigkeit U,, und t die Gitterteilung bedeuten. 
ie Punkte 0, b liegen auf den Vorderkanten der zwei ungestaffelten Gitter und die Punkte 
, c auf den hinteren Kanten. (Ist c=b, dann geht der Ausdruck (1) selbstverstandlich in 
ie Integralgleichung fiir das Einzelgitter iiber; und im Grenziibergang t—> 00 hat die Integral- 
gleichung genau die Form der Tragfliigelgleichung fiir zwei hintereinander liegende Einzel- 
fliigel.) 
Da die Tangente an die Kontur des Fliigelgitters parallel zur resultierenden Strémungs- 
tichtung ist, gilt in erster Naherung die kinematische St1émungsbedingung 


‘dz 
W(x) = Un (5 bi (3) 
() fale ; 
wobei 2(x) die Skelettkontur bedeutet. 


3. Strenge Lésung der Integralgleichung. Wir fiihren zunachst die Substitutionen 
i 78%, ny = €78%1, N= e28%s (4) 


| . 
.. Damit geht die Beziehung (1) in eine Integralgleichung iiber (abnlich der von H. Lomax? 
und P. Byrd® betrachteten), die sich in geschlossener Form lésen !48t. Unter der Abstrém- 
bedingung, daB 

F(a=0 und I,(c) =0 (5) 


| 1 J. Dérr, Ing.-Arch. 19 (1951), 5. 66. 
2 K. Nickel, Ing.-Arch. 20 (1952), S. 6. 
3 W. Spraglin, Flow Through Cascades in Tandem, NACA TN 2393, 1951. : 
4 Diese Integralgleichung ist nach der Singularitatenmethode und aus der Birnbaum-Glauerischen 
Theorie hergeleitet. Vorausgesetzt ist, daB in jedem der zwei Gitter die Profile aller Fligel gleich seien. 
Ferner werden die Wirbelbelegungen auf der x-Achse anstatt auf der Profilkontur angeordnet, da man 
sich hier auf sehr schlanke Formen beschrankt (d.h., man halt die Ordinaten der Fliigel fir erheblich 
klei ls die Fliigeltiefe und die Schaufelteilung). / oe 
Be 7 Lomax u. °P. Byrd, Theoretical Aerodynamic Characteristics of a Family of Wing-Tail-Body 
Combinations, p. 8, Eq. (19), NACA TN 2554, 1951. 
| 


192 Byrd u. Huggins: Zur Berechnung von Wirbelverteilung usw. Ingenieur-Archi 


ist, findet man schlieBlich fiir (1) die Umkehrformeln: 
vighhe (<< oe < aie 


_ __ 2g 1 /(Fe ga — Ts gx) (Tase—Iggx) |) 1 op i py 
BENG) emer Tren wre ie 


(Iq gb — Tq gx.) Vg gx, a | 
-|- [a—ctos>—an 1 ep W(x) won (6a) 


0 


~~ (Sq gx — Tg gb) Tg gxe W Peis | 
+ [oer *2)] Vee a (3) dp, | 


und fiir b< x <c: 


_ __ 2g 7/(Sg ex — Tq ga) (Tage — Tg gx) | 1 PAT wal | 
ae np (Ig gx — Tg gb) Tq gx ey cae 2 | 
ee = (Tg gb — Fg gx) Tg gm (x,) dx aa | 
+{l Etgal *)] Va Sie te ee a (6b) 
0 
(Xg gx. — Tq gb) Vg Exe )d 
+ [0—erg=— ‘lie a ee = (%2) | I 


Die Konstante A ist 


Whe yo—*aes (1 — &g ga) 
1—Zgqgb ° 


die GréBe (L,+L,) hat die einfache Bedeutung der Gesamtzirkulation um ein waagerechtes 
Fliigelpaar und wird nach Gl. (9) berechnet. (Die Lésung (6) fiihrt fiir c—b und 6 =1 un- 
mittelbar auf das von Dérr erhaltene Ergebnis zuriick.) 


4, Beispiel. Im Sonderfall W(x) = — Ua findet man aus (6) fiir das Vorder- und Hinter- | 


gitter die zwei Beziehungen 


Taye (Tg gb — Tq gx) Ta gx 
7 BD JV (gga — Tg gx) (Tg ge — Tg gx) ’ 


P(x) _ 2m0U. |e 
BD (Tg gx — Tg gb) Tg gx 


mit der Konstanten 


_ 7 Vb= S460) (L— aga) | [yaa ees eee 
- =a) 1— %q gb < (2g gb — Ig gm) Tg gx ae 


(8) 


-«[y/ Tg exe — Tg ga) (Tage— Tg gm) 4. 
(Xg gx_— Tg gb) Tg sxe 2 


Die Formelausdriicke (7) geben die Wirbelverteilung um diinne ebene Platten mit dem An- 
stellwinkel a. (Im Fall des Einzelfligelgitters, d.h. fiir c=b, erhalt man sofort die Ergebnisse 
von Grammel! und Nickel?.) 


5. Gesamtzirkulation um Gitterprofile von beliebiger Kontur, Bei Verwendung der 
Theoreme?® iiber die ,,Umkehrbarkeit“« der Anstrémung kénnen wir aus (7) fiir beliebige 
Skelettprofile eine allgemeine Formel fiir die Gesamtzirkulation 


L411, = [ LD (%) a +f L(x) dx, 


1 R. Grammel, Die hydrodynamischen ereniaeen des Fluges, S. 93—97, Braunschweig 1917. 
2 Siehe FuBnote BS) Whe 


% M, Heaslet u. J. Spreiter, Reciprocity Relations in Aerodynamics, NACA TN 2700, 1952. 
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sder den Auftrieb (Kutta- Jukowskischer Satz) 
H= 0... U. (Ly ++ Es) 
um ein Fliigelpaar gewinnen. Dabei kommen wir zu der Gleichung 


2m (Xgq gb — Bq gx.) Tq gx 
7 L hi = ly 98 § §%1) 29 8%, 
| ae Bo (Ug ga — Tq gx) (Tg ge — Tq gx) WA CE ar A 


(Ug gxg — Tq gb) Tq gx, | 
W (x5).0X, 
aa (yp PEE TICL LEERY TAMA Gee | 


Dies ist die GréBe, die in (6a) und (6b) eingeht. Die Gesamtzirkulation L, und L, um ein 
Einzelglied jedes Gitters lassen sich dann durch eine Integration aus (6a) und ‘(6b) Beseacnees 
Setzt man c=b, 6 =1, so ergibt sich zugleich die von iN che hergeleitete Gleichung.) Zur 
Ausfiihrung der Integrationen in allen der obigen Formeln ist es zweckmafig, die Substitu- 
tionen (4) zu verwenden. 


_ 6. Symmetrische Gitterprofile. Zum Schlu8 méchten wir hinzufiigen, daB das Verfahren 
auch fiir zwei hintereinander angeordnete Schaufelgitter gilt, deren schlanke symmetrische 
Profile Anstellwinkel Null haben. Hier bringt die Integralgleichung die waagerechte indu- 
tierte Strémungsgeschwindigkeit U(x) in Beziehung zu den Quell-Senken-Belegungen Q,(x), 
(x) der Fliigel. Diese Gleichung, die beinahe dieselbe Form wie Ausdruck (1) fiir die Skelett- 
inie hat, lautet 


U(x) = 58, i Ou(m) Ctg ge —a) dey +985 [ Qala) tg g(x — a) dey. (10) 
0 b 


obei Q(x) =2 U,, dz,/dx und z, die Ordinaten der Oberseite des symmetrischen Kérpers be- 
euten. Unter der Bedingung, da8 die Gesamtquellstarken verschwinden (weil die Profile 
eschlossen sein sollen), kann man wie fiir Gleichung (1) eine Umkehrformel der Integral- 
leichung (10) erhalten, die fiir gegebene Geschwindigkeitsverteilungen die symmetrischen 
rofilkonturen der zwei Fliigelgitter zu berechnen erlaubt. Voéllig ausgeschrieben sind die 
wei Gleichungen dieser Umkehrformel etwas langer als (6a) und (6b) und wurden deshalb 
uier nicht gegeben. Im besonderen Fall des ig atereieris fiihrt die Formel auf 


oa ak j Vis go —F5 em) Boum Ula W 
0) ie 


(Eingegangen am 31.Oktohber 1952.) 


Anschrift der Verfasser: Paul F. Byrd und Mary T. Huggins, 
Ames Laboratory, NACA, Moffett *Pield (California) (USA): 
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Ein Beitrag zur nichtlinearen Elastizitatstheorie. 
Von H. Jung. 


1. Einleitung. An der Oberflache eines isotropen Kérpers sind die Verschiebungen bz h 
ein Kriftesystem gegeben. Das elastische Verhalten des Kérpers werde durch eine Schubf 
funktion und eine Kompressionsfunktion gekennzeichnet, die beide aus Versuchen zu ermittel : 
sind. In Weiterfiuhrung der bisherigen Ansadtze von C. Weber!, A. Philippidis? und H. Ka i 
derer? wird die Voraussetzung, daB bei den Verzerrungen die Glieder zweiter und héherer Ord} 
nung gestrichen werden diirfen, nicht gemacht, sondern es werden erst die Glieder dritter un 
héherer Ordnung vernachlassigt. In Ziff. 2 werden die Verzerrungen und der Verzerrungss 
tensor hergeleitet, wobei vom deformierten Element ausgegangen wird. Die auf das defor} 
mierte Element bezogenen Spannungen sind in Ziff. 3 dargestellt. Mit den aus Ziff. 2 und é 
erhaltenen Ergebnissen wird in Ziff. 4 ein allgemeines Elastizitatsgesetz aufgestellt, und in 
Ziff. 5 werden die Zusammenhange mit dem Plastizitatsgesetz und dem Hookeschen Gesety 
aufgezeigt. Elementare Lésungen werden dann in Ziff. 6 bis 8 gegeben. Der ebene Verzerrunge} 
zustand wird in Ziff. 9 behandelt und daraus-eine technische Biegelehre (Ziff. 10) hergeleitet 


f 


2. Verzerrungen, Wird ein isotroper Kérper durch ein Kraftesystem deformiert, so geh 


der Punkt P des Kérpers mit den Koordinaten x, y, Zin den Punkt P mit den Koordinate 
x, y,ztiber. Bezeichnet man mit u(x, y, z), v(x, y, 2) und w(x, y, z) die Verschiebungen in dep 
x-, y- und z-Richtung, so ist der Zusammenhang der Koordinaten gegeben durch 


=x — u(x, ¥, 2), ¥ = ¥ — v(x, ¥, 2), Z = 2— w(x, y, 2). (i) 


Aus (1) lassen sich die Verzerrungen bestimmen‘. Ks ergibt sich 


Ou du\2 dv\2 , /dw\2 
ro2 = 25, (65) + (Ge) + (Ga) | 
eee 9 du\2 Ov \2 
v9 = 254 — [lay] +) 

— 9 Ow du\2 _ (dv\? 
Yes = Oz (3:) ! (] 


a), Cue ee du du __ ov Ov Ow Ow 
‘*Y “Oy ' 0% Ox Oy Oxdy dx dy’ 


__ du Ow du du Ov Ov Ow Ow j 
Vas-— gyi QR ADE as UR Dat ORR A (2b 


_ a | dw  dudu dv dv Bw Ow 
95 Ox Nay S Oyds 9 Oy den mayOne 


| 1 1 1 
ee Vx x 9 Vay are | 
Lae 1 1 | 
o= 9 Vey Ges 9 Vyz (3 
1 1 il 
9 Ves 9 Vyz 9 Vee 


C. Weber, Z. angew. Math. Mech. 28 (1948), S. 189. 

A, Philippidis, Z. angew. Math. Mech, 25/27 (1947), S. 31. 
HI. Kauderer, Ing.-Arch, 17 (1949), S. 450. 

G. Hamel, Theoretische Mechanik, S.190, Berlin 1949, 


ponprH 
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Die Invarianten des Tensors (3) sind 
1! 


Ji a. ye (eer Yyytb P| ’ 
—— 1 9 2 
me aa Als (Vex Yyy Yaw Vex Tey tea Px le Vys) ? 
oa y ) , (4) 
1 xX é xy } NZ 
Or ae es oe 
Vxx Vyz Vsz 
Fiibrt man die Tensoren 
1 1 1 1 i ) 
Ear 0 0 | 9 Vex recs 9 Vey g Vas | 
| 1 1 
1 1 1 | 
So= | 0 mls 0 und: 5O.= ! 9 Vey a Yxy— ye Ss Q Vy# | (5) 
1 i 1 1 1 
0 0 aa | g Vee =e Vyz “9 pee ey 
ein, so laBt sich der Verzerrungstensor in der Form 
6=6+6 (6) 
darstellen. Die durch die Verzerrungen hervorgerufene Volumenanderung ist 
Ap =1—JI—2),+45,—8 Jy. (7) 


Bezeichnet man mit J; (i = 1, 2, 3) die Invarianten des Tensors G, so ist die durch die Kom- 


ponenten von S bedingte Volumenanderung 


sis Mee (8) 
Werden in (2a) und (2b) die quadratischen Glieder vernachlassigt, so wird 
eM 


da J, und J, nur quadratische Glieder enthalten. 


3. Die Spannungen. Schneidet man aus dem deformierten Kérper ein Volumenelement 
heraus, so lauten die Gleichgewichtsbedingungen fiir dieses Element 


00x OTx y | Ox ad | : 


raat 


dy | @ 
OTx ¥ doy Oty: l 
| PUT eg en | (9) 
Ox Oty x 00: vl 
Ox ain dy aH Oz is | 


Bei der Aufstellung der Gleichungen (9) ist angenommen, da die Massenkrafte und andere Vo- 
lumenkrafte vernachlassigt werden diirfen. Die am Volumenelement auftretenden Spannungen 
i sind die Komponenten des Tensors 


[ o, Bice Tae 
ealme & Gy \. | Tys (10) 
Tx 2 Tyz Ox 


-mit den Invarianten 


Pi 6,0, 4-6, = 3 0, 
P, = KO Oy ai Oy Oz + Oy Oz ee Pe yz 

On Tx Tz (Pt) 
Py =|Txy oo Vys 

Txz Ty x O; 
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Spaltet man den Tensor (10) in einen Kugeltensor ‘f), der fiir das Element den hydrostatischen 
Spannungszustand gibt, und einen Deviator S$ auf, so 1aBt sich (10) in der Form schreiben 


% 0 0 Ox — 99 Txy Tx x = 
aS == || () Co 0 -+- Tx y Oy ma 0 Vy x = Bo ‘gra 8 e (12) 
0 C9 | | Tx x Ty 0, — 03 J 


4, Das Elastizititsgesetz. Bezeichnet man mit y; (t = 1, 2,3) die Hauptdehnungen, so 
lassen sich diese aus (4) bestimmen. Ks gilt 


Yan + V¥yy + Yee = V1 + V2 + Ys > \ 


(13) 
Yxx Vyy Mee Vex 59 Vex Vee Mee ole == NA) + Vos +7715 J 
You Yay Yur yee 7 0 0 
Yay Yyy Vyz| = 0 v2 0 (14) 
Vxz Vyz Vez 0 0 Y3|° 


In entspirechender Weise ermittelt man aus (11) die Hauptspannungen (0;, 62, 05). 
Nach Ros und Lichinger? gilt 


2 Lee (Rene 
Vo — 6)? ++ (6, 0)? + (63 — 0%)” = All (n 3 J,) Al (v2 3 J,) t bys oy) ). 
oder in den Tensorinvarianten aus (3) und (10) geschrieben 
y2 yz PP? =P» = (2 + Je a (15) | 
Eine Schubfunktion wird durch 
aps 
el? ag hee 
Vitis, 


definiert. Fordert man, daB das Elastizitatsgesetz beim Ubergang zu kleinen Formanderungen 


das Hookesche Gesetz ergibt, so wird der Zusammenhang zwischen §§ und © am zweckmifig- 
sten in der Form gegeben 


(16) 


Z =6,(/2/2 n—4,). (17) 


Fiihrt man noch eine Kompressionsfunktion fly i) ein, so 1aBt sich der hydrostatische Span- 
nungszustand geben durch 


By = Sho(z i). (18) 


In (17) und (18) wurde die Voraussetzung gemacht, daB® die Hauptspannungsrichtungen mit 
den Hauptdehnungsrichtungen zusammenfallen. 
Das Elastizitatsgesetz lautet damit in Tensorform 


B= Es(2/tn— 4) +64 (45) (19a) 


und in Komponenten 


t 1 Sy lee 1 1 

Oy =(4 yy — 4 i) {V2 \3 Ji J) 3 Shs a (19b) 
1 1 a pecmaer = ft i 

0; =( Vez 3 5) «(V2 We ws 3,) | 3 omar ip 


aon a ae u. A. Eichinger, Versuche zur Klarung der Frage der Bruchgefahr, III. Metalle, Ziirich 
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Tx, = & vee (12 es dt = 3) ° (19c) 


1 — 1 
G ty =H Yet (2 PEF Je). | 


Definiert man eine Querdehnungsfunktion durch 
1 Sai omt 
fr e J,) =~ (12 te R-J,) 
= = _ (20) 


(Vins) 


so lassen sich die Gleichungen (19b) in der Form schreiben 


Vy 
, 


; | ce ee 
o. = (40x +2) (12 VE It J). 
= a a i ne 
Oy = GB Yyy : 3) a (V2 y+ Jf ee 1,) rs Hf (21) 
1 =n Vay ee 
Ora ( 9 er = 3 5) a(2 Ve Ji ee 3.) ‘ 
Setzt man 
DE OS Naeem aie 
Ce et Pale (22) 
& =F (}P), ey) 


wobei G die Umkehifunktion von g und F, die Umkehrfunktion von f, bedeutet, so 14Bt sich 
das Elastizitatsgesetz in der Form schreiben 


6 = Hellz /t Pe P) +8 F,(4 7). (24) 


In Komponenten erhalt man 


| ) 
LecloERlo(AVEA—R)tan(he).| es 
ee 


(26) 


(fs ee (5) (27) 


14 
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ein, so geht (25) tiber in 


Sry = (0, 7 iH a antes “ | 


5. Spezialisierung des Gesetzes. Das durch (19a) bzw. is gegebene Elastizitatsgesetz muB | 
beim Ubergang zu kleinen Formanderungen das Hookesche Gesetz ergeben. Entwickelt man 


die Schubfunktion 
= 1 
6/2/34 2-4) 


in eine Potenzreihe nach der Veranderlichen 


25 R— h. 


so ist 
2 2 5 ces Ea? 2g.) ae 29 
(3/4 J J. Jz) = & +8172 z Ji— + B° 2 as 2) 39 ° (29) 
Die Kompressionsfunktion wird ebenfalls in eine Potenzreihe nach = J, entwickelt; es wird 
u i 1 
f(z h)=~t+aghtagdit (30) 


Geht man mit (29) und (30) in (19a) und (19b) ein, so kommt 
oa ‘ 2 

Ox g Vex 3 J) (8. + a V2 V4 3 J, Jo + Vedas S (e+ az a ee a 

1 1 

a, =(3 Yyy aAlt +g, /2 ee he 

1 

3 


) 
/ 
(% +a V2V2Jt— t+ ipa TAC) vac ), 


“ease \ 


(31) 


a 
Tee = ra (co +a: y2)/4 1 er ), 
y= 5 Veal Bo gy) See . 


Streicht man in (31) similiche Glieder, die Quadrate und héhere Potenzen der Dehnungen ent- 
halten, so erhalt man das Hookesche Gesetz. 

Soll aus (24) das Hookesche Gesetz abgeleitet werden, so geht man mit den Reihenentwick- 
lungen 


= 1 = 1 1 
(/3/F Pr P,)= 6,4 G2) 4 Pr P, +6,-2(} Pp ) eee (am 
1 1 1 
(= P| 0, ee erg eee (33) 
in (25) und (26) ein und erhalt 
1 P mye Ls 1 
te rs (6 G22 P?— P, 4 | | ene toe +), 


(34) 
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L © Wise eae 


Soll (34) auf das Hookesche Gesetz fithren, so miissen in (32) und (33) die Glieder 


x 1 
G, (v2 Le p? — P,)< Ceeind = E i | < Ey (ES) 


. damit in den Reihenentwicklungen (34) die Glieder héherer Or dnung gestrichen werden 
onnen. 


Um die klassische Plastizitatstheorie zu erhalten, setzt man 


wes 
oe 


pale 
(2 Ve R— 4] = 


09 toe J = 08 


BIT, 
Die Spannungen werden damit 
1 §-1 
Ox — 09 ==> Vex = 
0 9 y /2 = ve 
noe ee 
Ap 0 9 ro 2 V7 
7 2 
mies 1 §-1 
Ox 09 2 Van 2 ine a ° 
1 oy 
ty Pp 
2 
1 8-1 
Tyz = Vane ae r) 
See ays, 
1 8-1 
Tez — Vex = 
Ape Y 2d 


6. Gleichférmige Spannungszustinde. a) Der hydrostatische Spannungszustand. 
Setzt man in (25) und (26) 


Ox, = Oy G60, =s koust:, 


Txy= Tx, =Ty, = 9, 


so werden die Gleichgewichtsbedingungen erfiillt, und es ist 


Vax =Vyy = Yer = 2K, (35) 
Vay = Var =Vy2 = 9. (36) 
Geht man mit (36) in (2b) ein, dann erhalt man 
du | av du du dv av dw dw 9 | 
dy | dx axdy Oxdy Ox dy ? 
ee eu 
dv , Ow du du dev dw Ow 9 | 
dz | dy dy Qz Oy Oz Oy Oz j 
Die Gleichungen (37) werden befriedigt durch den Ansatz 
b= 0%, VS Cy,  W = Cy%. (38) 


14* 
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Aus (38), (2) und (35) erhalt man 


Ci S=iCy a= Cy == Co, (39 
2 Kg sete 0 


Aus (39) ergibt sich die physikalisch mégliche Lésung 
ec =1—yl—2 K,=1—Jl—20,F, (%). (40 


b) Der einachsige Spannungszustand. Setzt man in (25) und (26) 


Oy = 0, =Txy Tuz =Tys Om 


Cy ons 
dann ist 
1 
Oo are 
und 
et eirie Fa ee 
75 Nae — a Me (Xo) Ta % 2 (9) - (41 
il 1 1 1 
2 Vy y 2 Vuxz a a 3 wo G (29) +z Ho Fs (xo) © (42) 
Vege Ve ee (43) 
Mit (43) und (2) wird 
Ul Ce el nee (44) 


Aus, (44), (41) und (42) erhalt man mit (2) 


4 2 A 
q=1 -/i — 3 #0 © (%o) — = Ho Fo (%) ’ 


2 2 
Coe oe 1-V1 by Ho © (%) — 3 Xo Fs (xo) : 


c) Der zweiachsige Spannungszustand. Setzt man in (25) und (26) 


Ox = Try = Te, =Ty, = 0, 


Oy = Xo> O;, = *4> 
dann ist 
1 


= 3 (% +4) > P, = x %- 


eri 


Die Verzerrungen sind 


1 1 ae ) 
g vee = Pee + 4) o(/2 y= Mo Hy +A) as = (% + %) F, Be (Ho +1)] » (49) 


3 
zi =(5%— 3) CVF Pm FA) + hate) Fa(L be +m))> (6) 
xs) G (V2 x2 Xo 4) ; : (% + %) F, e (% + “)) ‘ (47) 


U=C%, UV=CY, W= Cz 


=1—[1 +5 bo +m) 6(VZ (aim FA) — 2 by +10) Fa(L Gm +] 
Gea 1 (5 ite : m) GC We x2 — Hy Hy, + 72) — 2 (29 ++ 2) ae (%) + “))) > 
Co f e ey : 4) G (V2 V2 — Hy %, + “| — & (%) +%,) F, . (%) + ~))] 


ist damit die gesuchte Lésung. Der dreiachsige Spannungszustand 1aBt sich auf dieselbe Weise 
behandeln. 
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d) Der einachsige Schub. Setzt man in (25) und (26) 
Ox Oy O; Tx y Tx x 0 ’ Tyz = > 
so werden die Gleichgewichtshedingungen erfiillt. Die Verzerrungen sind damit 
Vx x Yyy Vex Yxy Vuz 0, Vyn = %y GC (4) (48) 


Setzt man (48) in (2) ein, so erhalt man das Gleichungssystem zur Bestimmung der Verschie- 
bungen 

eS eiaeice 

BUTT ey 

Slee oan 


Ou ov du du Ov _ Ov Ow Ow — (|) 
dy ' dx ax dy ax dy dx dy a 


Ou Ow ou Ou Ov _ Ov Ow Ow 0 
Oz Ox Ox dz Ox Oz ONGNOR! Pia mas 


Ov Ow du du Ov dv Ow dw 
ors Oy oy cosety hk on™ bye gre A” 4) 


Um das Gleichungssystem (49) zu befriedigen, macht man den Ansatz 

Si= (() VU =C,2, W= C52. (50) 
Geht man mit (50) in (49) ein, so erhalt man 
; Co y2 cs — e2, 

Cy = 2,6 (mm), ¢, =1—Y1 +x,G6(m). 

_e) Allgemeiner Ansatz fiir die elementaren Liésungen. Samtliche hier durch- 
ec chneten Beispiele und noch weitere lassen sich aus dem Ansatz fiir die Verschiebungen 
) u=ax-+ay +432, 
) : fee en| (51) 


W = Cy, X—- Cy YY, -|='€3% . 


herleiten. 


7. Torsion. Bezeichnet man mit P(r, p>» z) einen Punkt vor der Deformation und mit 
P(r, vy, 2) diesen Punkt nach der Deformation, so ist der Zusammenhang der Zylinderkoordi- 
mnaten gegeben durch 


(52) 


wenn mit 
u die Radialverschiebung, 


rv die auf dem Kreisbogen r = konst. durchlaufene Strecke, 
w die Axialverschiebung 


bezeichnet wird. Aus (52) erhalt man die Verzerrungen? 
2 2 (dv\? Ow\2 
Ugeud) & or}? 
ae! u \2 dv\? (5 u \2 1 /dw\2 1 /du\2 
J +250 —F) — Ge) (17) — oe) — a): 


me) () (raw (5) (53) 


= 
x 

x 

| 

D| WY 
TIS 
D| 
718 
SSS se 


1 FuBnote 4, S. 1. 
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1 fou 2 Ov du Ou Ow Ow ff 2 Ov Ov 
Yr@ => le pene we) Or Or dp Or Op ar) y 


Ow Ou Ou Ow Ow Gi 1)? Ov Ov 
Yrs 5, VT Or Or dz Or Oz Or Oz” 


1 [dw 20v  dudu dw Ow 9 Ov Ov 
Vox => Ee sai ira eee Ce oe ree ae Beil ; 
In Zylinderkoordinaten lauten die Gleichgewichtsbedingungen 
00r i a Otre Ono oe 
~ oi t r q Oz r 1 
ee 1 Pe | ee ao Zi 
Fir 7 Op | = > ( (54) 
Otrz | 1] OTp« | 00: | ale Saye 
On. lire Oph osu cr 
Setzt man in (54) C;, = 09 10) = a7 = Ue 
so gehen die Gleichgewichtsbedingungen iiber in 
Otr OT pz 55 
Op ay Op sae ( ) 
OTr » OT pz | Otro 56 
or 1 az | r r , (56) 
Mit Tro =P, (7), Tox = P,(r) 
ist (55) erfiillt und (56) geht tiber in 
_ +272 ‘ (56a) 
Macht man fiir die Verschiebungen me Ansatz 
=P =, v=, % 
und vernachlassigt (5) (Theorie erster Ordnung), so erhalt man aus (53) 
a 
Yrr a Yoo Vez Yro Yrz 0 > | e 
= { (57) 
ern 


Mit (57) wird 
> gl TC, 
7 ae(T aR (58) 
Durch (58) ist die Spannungsverteilung in einem kreiszylindrischen Kérper gegeben, der an 
seinen Enden durch zwei Momente verdreht wird. Bezeichnet man mit r; den Innenhalbmesser — 
und mit r, den AuBenhalbmesser eines Rohres, so wird das auf das Rohr ausgeiibte Drehmoment 


pa er wae (59) | 


Der Torsionsversuch ]48t sich damit zur berg enne der Schubfunktion nur im Rahmen einer 
Theorie erster Ordnung heranziehen. 


8. Der Zugversuch. Fiir den rotationssymmetrischen Verschiebungszustand erhalt man 


aus (53) mit v= 0, p= 
Ou du\2 Ow\2 
DT oe ) 2 (Fr) % 


(60) 


ei Ow Ou Ou Ow Ow 
ie 5, Vay ar oy mere ora 
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Die Gleichgewichtshedingungen (54) lauten dann 


Oor | Or — 09 Otr z 


or r Oz 


2.0. 
(61) 


Otrz ! 00 | Trz 


or az | gee bx 


 Setzt man 


O,=%}> G, =O, = 7%;, =—9- 


so weiden die Gisicleewichtshedingungen (61) erfiillt, und man erhdlt fiir die Verzerrungen 


1 1 1 I 

5) Yrr 3 Yow = — 3 %4 & (%) 3% Fy (4) » 

(62) 
; ae 3 x, G (*,) +z 4 (%) | 


u=c Tr, 1 = Ce 
so wird ? 

} 

Ver BG = ey © (4) 5) 

ee ay 6 pas 

Yoo = Ea 0 Se Rar ag (%;) + 3 2 (%) » (63) 

4 

Vax = 2 €, — 3 = 4G (%q) Pa Fe (m) | 

J 


Mit (63) ist 


(64) 


Der Zugversuch eignet sich nur im Rahmen der Gleichmafdehnung zur Bestimmung der Ver- 
festigungsfunktion, da beim Auftreten einer Einschnitrung die Bedingung 


0, =9,=0 
nicht mehr fiir alle r und 2z erfillt ist. 


9. Der ebene Verzerrungszustand. Macht man fiir die Verschiebungen den Ansatz 


u=u(x,y), v= v(x, y), w=0, (65) 


ree = 255 |(5s) +(e) | 
maak (Bh) 


so geht (2) iiber in 


(66) 
du , Ov du Ou Ov Cv 
Yey— Gy | Ox  Ondy  Oxdy’ 
Yur = Vys = Yar = 9. 
Die Gleichgewichtsbedingungen lauten damit 
CLS CE ee 
0x ' dy | 
(67) 


wax 4 Sy 0, | 
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Mit (66) werden die Spannungen 


1 = 
Oy = 6 (2 Vie = Yyy) «(V2 


iii 
pe yy es) 82/3 fase. 
ones —FOetyde(BY3 J i,) | = ene (68) 


1 Ne hae ae 
he =F reya(2Vt Jr i), 


ah 


— 


a 
re 


i 2 i 
1 = Yau Yyy) > J = GZ (vex Vr — Vey) : 


AuBer den schon bekannten elementaren Lésungen lassen sich auch fiir den ebenen Ver- 
zerrungszustand keine weiteren exakten Lésungen auf einfache Weise finden. 
An dem durch Abb. 1 gegebenen Belastungsfall soll nun gezeigt werden, wie technische 
Beispiele naherungsweise durchgerechnet weiden kénnen. Fiir diesen Belastungsfall lauten 
die Randbedingungen 


zh - aA fy ae (69) 
AL Zz : 
eles: Soa Macht man die eer da t,, identisch verschwindet, so 
7 ania fithrt (67) auf 
Lr 
y dx = D,(y), oy = P(x). (70) 


Abb.1. Setzt man 
Gyace Ole 


so sind die Randbedingungen (69) erfiillt. Aus der zweiten Gleichung (68) folgt 
ik a See + + 
$2 ry — ves) 8(V2/ 4 BL + AS) =0. (71) 


Um die Gleichung (71) nach y,, auflésen zu kénnen, macht man fiir die Verschiebung den 
Ansatz 


Wie (72) 


Nice ees dv\2 /du\2 /dv\2 d Ou du dv Ov m : 
ernachlassigt man az! Vay) ay}. neo de Oy ox oy (o Cobale aman mit (72) 
Ving <= Oy ey an 


Mit der Annahme, dafi die Schubfunktion und die Kompressionsfunktion nur von der Haupt- 
verzerrung 7, abhangt, laBt sich aus (71) y,, naiherungsweise bestimmen: 


—__slEons) alban) 
| der (ye ssl ae ys] 
Aus (73) ergibt sich dann 
[ater tl) ae 
; meres il) +4,(4 vss) : 


(73) 


(74) 
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Die Spannungen werden damit 


7 g pile We x 1/0) i ) x6 x 
| at (Ve ae ai (75) 
. ea 


| | 5 (\/% tre) teres a ee ee 
Bonne apo er 


Beriicksichtigt man bei den in (75) auftretenden Tensorinvarianten i und of noch das durch 
(73) gegebene y,,, so ist 


dk ee — ; 
yA 4 2 
2e(I/2 ives! +f(d res] Abb. 2. 


10. Der unendlich lange Plattenstreifen. Der Belastungsfall sei durch Abb. 2 gegeben. Die 
Abmessungen der Platte sind 
h h 


<=) Sa Oe SS he ——s 7 <=-+ 


Auf den Randein x = + akann die Platte als starr eingespannt oder frei aufliegend angesehen 
werden. 


Aus zahlreichen Versuchen ist bekannt, daB auch bei groBen Ausbiegungen die Dehnungen 


Al 
a == ay (76) 
sind. Mit (76) erhalt man 
J il 
| OP et Vox ieee Ne (77) 


Geht man mit (76) in (77) ein, so wird, wenn in den Reihenentwicklungen nur die drei ersten 
Glieder beriicksichtigt werden, 


Yax = 2G ¥— cpy?. (78) 
Der Ansatz (76) fiihrt damit auf dieselben Verzerrungen wie der Ansatz (72). Mit (78) werden 


die Tensorinvarianten, wenn y,,y vernachlassigt wird, 


j.=aqy—tay, f,=0. (79) 


Die Kompressionsfunktion und die Schubfunktion werden mit (79) 
/é = 1 
fly 4) =o + a) lar—zar) to 
ie = o (80) 
91/1 fz 2 | 1 242 | 2 ae iar | 
B\V2z J2— 2) =2lV 3 jay—g7ie|)=a tas 


] 
CRN iy 1 
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Geht man mit (80) und (78) in (68) ein, so erhalt man, wenn bei den Reihenentwicklungen nur 
die beiden ersten Glieder beriicksichtigt werden, 


1 1 oF 
o.=1(2q47— a y*—t 995) |e tale | 
zx 1 
& tal/z (a9 — det y')], 

1 Nee pe 

Oy =| (2 Vyy — 24 ¥ + ETN") [80 FAY = 
Tei 1 es 1 

alee E Nyy Fay — 5 cb 9) |e +¢,|/4 (ax—fet st}. 


Vernachlassigt man den Einflu8 der Querkraft, so ist 


1 
Gy te 


(81) 


{l ] 1 
3 (ey 9 cy? 9 V9) 


1 
1) 


Ona 


Aus (82) kann dann y,,y bestimmt werden: 
Ee 
Ei ey Via les 
Bot > : eo + 2 ae 81 


Setzt man (83) in (81) ein, so kommt 


1 
yy = (ay 5) iy’) 


1 
Bi Jem yg 


asta] 

ay— yey —a(ay—tety)] (84) 
if 

( ‘ 

talay—yety'] 


Aus (84) erhalt man durch Reihenentwicklung nach c", wenn c” (n > 2) vernachlassigt wird, 


1 1 j 2 
Ox =tlax—fay)p & + +2 2 81 
1 1 cm 
3 (ay 5) ct9*) 60— & +/2 B 
2 
&0 — &o +) E 


ae 
289+ & + le 2 $1 


1 
CR oaks hr 


x 


1 
ay— sey? 


2 
a 0 81 il 


1 (2 8 + &y)” Sh ey)” 


ee | ¢2 IAPs 4 5 
Ogee, Sy Y. 2 2) + e% tery |y, ype ee ]. (Sam 
Mit (85) wird das auf einen Querschnitt x = konst. ausgetibte Moment 
+ h/2 
M(x) = [ yo,dy ='de, + Be, (86) 
—hj2 


wobei 


V5 | 
ye TE 8 do dee)? alee ee | 


36 2g + ey 2 32 2g, +e 
gesetzt ist. — | 
Aus (86) ergibt sich 
ar AO) ae Vatee = (oq 
Setzt man (87) in (76) ein, so kommt 
2 M(x) 
Ex i a5 a ema Ay. (88) 


Mit der in der technischen Balkenbiegungslehie itblichen Annahme ist 


dx exdx A Az M(x) 
Q y =| 2B Viet B Jas. (89) 
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Bezeichnet man mit & die Durchbiegung der Platine, so erhilt man aus (89) 
dE 


dx A At M(x) 
eee el (90) 
d&\2\3 /2 2B 4 B2 B 
(+ (Fy) 


‘Da in (88) ¢ nicht mehr klein von erster Ordnung ist, kann in (90) die Kriimmung 9 nicht durch 
die zweite Ableitung ersetzt werden. Die Annahme 


1 _ @é 
o.. dx 
fiihrt auf eine Theorie erster Ordnung. 
Setzt man 
LE, 38 (a 
oo. dx? (1 2 (=) 
in (90) ein und entwickelt den Wurzelausdruck in eine Potenzreihe, so erhalt man 
as 3 (a\? l _ B yp wee 91 
(1 mate ! -)=aM oe aie (91) 


Beschrankt man sich in (91) auf Glieder erster Ordnung, so geht (91) in die bekannte Diffe- 
rentialgleichung der technischen Biegelehre iiber. Die Lésung der Differentialgleichung (91) 
wird am zweckmafigsten mit Hilfe der Stérungsrechnung durchgefiihrt. 


(EKingegangen am 3. November 1952.) 
Anschrift des Verfassers: Dr.-Ing. H. Jung, Stuttgart-N, Robert-Mayer-Str. 66. 
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Zur Berechnung von Pendel-Schwingungstilgern. 
Von K. Desoyer und A. Slibar. 


1. Einleitung. Als Fliechkraftpendel ausgefiihrte Schwingungstilger werden nach den ein- 
fiihrenden Arbeiten von R. Sarazin!, E. S. Taylor? und O. Kraemer® sowie weiteren Unter- 
suchungen. verschiedener Autoren* im modernen Motorenbau immer mehr verwendet. W. 
Schick* rechnet unter der Voraussetzung kleiner Pendelausschlage die Einwirkung von meh- 
reren Fliehkraftpendeln auf eine einfach zusammenhangende elastische Kette als Ersatz- 
system fiir die gegebene Kurbelwellenanordnung eines Reihenmotors. A. Kimmel und I. 
Lutzweiler* fiihren diese Naherungsrechnung fiir eine als mehrfach zusammenhangende 
elastische Kette aufgefaBte Kurbelwelle® durch. Die Rechnung fihrt in beiden Fallen zu der 
als ,,Abstimmbedingung“* bekannten Beziehung 

(6 == 0 (1) 
zwischen der Kreisfrequenz 2 der erregenden Momente, der mittleren Winkel geschwindigkeit 
W) der Welle und den Abmessungen L und I (siehe Abb. 1). Da die Kreisfrequenz Q der er- 
regenden Momente ein v-faches der Winkelgeschwindigkeit w, der Maschine ist, vereinfacht 
sich die Gleichung (1) zu: 


~|h 


(1a) 


vy durchlauft dabei fiir Zweitaktmaschinen die Werte 1, 2, 3, .-., fiir Viertaktmaschinen die 
Werte 4/5, 1245 222 eee 

Bis jetzt wurde angegeben, zur Erreichung einer befriedigenden Wirkung derartiger Flieh- 
kraftpendel die Gleichung (la) je nach der zu tilgenden Ordnung einzuhalten und im iibrigen 
die Masse des Fliehkrafttilgers méglichst groB auszufiithren®. Eine genaue Angabe fiir die giin- 
stigste Masse des Tilgers zu liefern, waren die bisherigen Naherungen nicht imstande, da bei 
Annahme kleiner Pendelausschlage in den damit linearisierten Bewegungsgleichungen der von 
der Corioliskraft herrithrende Term nicht mehr auftritt. Gerade dieser Term aber ist es, der 
die Bedingung fiir die Bestimmung der optimalen Pendelmasse liefert. Dies soll im folgenden, 
zunachst fiir den einfachen Fall des als mathematisches Pendel aufgefaBten Tilgers, gezeigt 
werden’. 


2. Die Bewegungsgleichungen, Der Rechnung liegt eine Anordnung gema4®B Abb. 1 zu 
Grunde. An einer Schwungmasse vom Trégheitsmoment J sei im Abstand L von der Dreh- 
achse ein Pendel der Masse m und der Lange / reibungslos angelenkt. Auf die Schwungmasse 
wirke ein auBeres Drehmoment 


M=M,+AMsinvg, (2) 


wovon der konstante Betrag M, zur Aufrechterhaltung der konstanten mittleren Winkel- 
geschwindigkeit w, der Schwungmasse, also zur Uberwindung der Arbeitswiderstande, ab- 
gespaltet werde. Die Bedeutung der weiteren Bezeichnungen geht aus Abb. 1 hervor. 

Die Bewegungsgleichungen fiir die Bewegung eines Massenpunktes m in einem rotierenden 
System lauten in vektorieller Schreibweise bekanntlich: 


mb, = 2) B + m(t xX to) + m (tv Xt) X tw +2 m(b,4 Xt) . (3) 


1 R. Sarazin, D.R.P. 597091. 

2 E.S. Taylor, J. Soc. automot. Engr. 38 (1936), S. 81. 

° O. Kraemer, Z. VDI82 (1938), S. 1297, Z. VDI 83 (1939), S. 901, Motortechn. Z. 1 (1939), S. 3. 

* W. Schick, Ing.-Arch, 10 (1939), S. 303; A. Kimmel und I. Lutsweiler, Ing.-Arch. 12 (1941), S. 100; 
R. Lambrich, Das Fliehkraftresonanzpendel. Jb. dtsch. Luftfahrtforschung 1940, S$. 1170; H. Kam- 
merer, Die Eigenfrequenzkurven bei Drehschwingungssystemen mit Fliehkraftpendeln, Jb. dtsch. Luft- 
fahrtforschung 1940, S. 1180. 

5 R. Grammel, K. Klotter und K. v. Sanden, Ing.-Arch, 7 (1936), S. 439. 

° H. Schrén, Die Dynamik der Verbrennungskraftmaschine, Wien 1947; K. Haug, Die Drehschwin- 
gungen von Kolbenmaschinen, Berlin 1952. 

* Die Behandlung des Tilgers als physikalisches Pendel soll spiter erfolgen. 
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Darin bedeuten »b,,, bzw. b,,; die Relativgeschwindigkeit bzw. die Relativbeschleunigung 
des Massenpunktes m gegeniiber dem rotierenden System, ferner ty bzw. vy den Vektor der 
Winkel geschwindigkeit bzw. der Winkelbeschleunigung des rotierenden Fiihrungssystems und 
t den vom Drehpunkt des Fiihrungssystems zum Massenpunkt weisenden Radiusvektor, 
x das vektorielle Produkt. Die Terme auf der rechten Seite von (3) bedeuten 

> $ die vektorielle Summe aller an der Masse m wirkenden physikalischen Krafte, 


m(tX ty) die zufolge der Winkelbeschleunigung des Fithrungssystems auftretende Schein- 
kraft, 


m(to Xt) xXto die Fliehkraft, 

2 MV, X10) die Corioliskraft. 

Die Orientierungen, Richtungen und Betrage der drei letztgenannten Scheinkrafte sind in 
Abb. 1 eingetragen. Sieht man wie iiblich vom Gewicht des Pendels ab und spaltet man die 
Gleichung (3) in Komponenten in Richtung BA bzw. senkrecht dazu auf, so folgen die Glei- 


chungen 


m1 = S$ + mr sin y— mro*cosy—2mlwé, | 3) 
ml § =— mra@cosp— mrw? sin y, | 

worin S die Stabkraft im Pendel bedeutet. Aus Abb. 1 liest man die Beziehungen ab: 
-siny LL 
sind” or 
Die Bewegungsgleichung fiir die 
Schwungmasse mit angekoppeltem 
~ Pendel lautet 

Ja =AMsinvgy + SLsin#. (5) 
Aus (3), (4) und (5) folgt 
(J+ m L?sin?3) @ = AMsinyp 

+ LImsin #. (10? +w? (l-+ Leos?) 

+2lw#], 
(19 +0 (l+Leos 9) +o*Lsin d= 0. 
- Trennt man die Winkelgeschwindigkeit 
) w in den konstanten Anteil Wo und die Abb. 1. Schema eines Fliehkraftpendels an einer Drehmasse. 


| Stérung dw: w =a,+ 4o, (7) 


,» reosy=1-+ Leos}. (4) 


(6) 


dann ist wegen @,—0 die Gréke @ durch Aq zu ersetzen. Fiihrt man (7) in (6) ein und streicht 
wie itiblich Aw gegen w, so folgt 


(J + mL? sin29) Ad = AMsiny wt + Lmsin 9 - [102+ a2(1-+L cos 9) + 21099), | 
19 + Aw@(I+L cos9) +02 Lsind =0. 


3. Aussagen der linearisierten Gleichungen: Wird in den Gleichungen (6’) der Ausschlag- 
winkel # als klein angenommen, so gewinnt man die linearisierten Gleichungen in der Form: 
JAo = AMsiny ot +mL(L +1) a8, i (8) 
) 16) +-(L +1) do + Lak od = 0 f 


mit den fiir den stationdren Vorgang giiltigen partikularen Lisungen 


(6’) 


AM: (L-+l) ) 
Sa we (J(L— vl) + mL(L + 1)*) oe | ie 
( 
AM - (L—7?1l) ; 
ee ome ee. | 


Die zweite Gleichung (9) liefert die bekannte Abstimmbedingung (la), jedoch bei erfiillter Ab- 
stimmbedingung keine Aussage iiber die Gréfe der optimalen Pendelmasse ; fir nicht erfillte 
Abstimmbedingung folgt daraus die bekannte Richtlinie, die Masse m méglichst groB zu 
-machen. 
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4, Lésungen fiir gréfere Pendelausschlige. Fiihrt man in die nicht linearisierten Glei- 
chungen (6’) die dimensionslosen Abkiirzungen ein: 


I AM _ 
i =A ? 7 a 0 ? J w ad ? 
und kennzeichnet man die Ableitung d/dr durch einen Strich, so nimmt das System (6’) die 
Form an: 
(1 + Osin? 9) - dw’ =a, K siny t +a) O sin -[AW?+A + cos +218], 
LO’ + sin & (11) 
9 A+ cos 3 


Wyt=t (10) 


Aw’ = 


Beriicksichtigt man, da® 3, im BogenmaB gemessen, praktisch stets kleiner als 1 ist und ver- 
wendet man demzufolge die Entwicklungen 
5 he 7? 
sin § = o ——, cos 0 1, 
so geht die zweite Gleichung (11) mit 


1 1 1 @ 
pene ee Ty ay 


iiber in 
AOS ip sue ee paca oe) (12) 


Wo 


wobei zur Abkiirzung gesetzt wurde 
A ] 2—A A 


VSTigg?/ MT SS Teaco. Sie 
Die erste Gleichung (11) geht mit (12) tiber in 
a, 8’ +a,h + 4,8 +a, 0'R + 4,973 +a,00 = Ksinyt, (13) | 
wenn die Abkiirzungen 
tA eay a, = ¢,-c0% 
a, =c,—(1+A)0, a, = —AO, (13a) 
a, = 0+2 (60, +4 +4), of H= 270 


verwendet werden. 


Die Gleichung (13) stellt eine nichtlineare Differentialgleichung fiir #(t) dar, deren fiir den — 
stationdren Vorgang giiltige partikulare Lésung im folgenden durch Iteration gefunden wird. | 
Eine partikulare Liésung des linearen Anteiles von (13) 


a, 0, + 4,8, = Ksinyt 
ist ' 
Oe 2 Ksinyt. (14) 


Bringt man in (13) die nichtlinearen Terme auf die rechte Seite und setzt hierin fiir } und dessen 
Ableitungen die erste Naherung 3, gema8 (14) ein, so ergibt sich fiir die zweite Naherung #0, die 
Differential gleichung 


a ye” rs] PEN Ke A a,— v? a, K3 cag) 
1 Gs — oe eee eas sin" VT 
2 1 5 
ee EM Neve (15). 
Bic ean Go cos PC BIDET — 
eae 1 


v Ag 


(ag tage KP sin yr cos. | 


Diese 1a8t sich mit Hilfe der Identitaten 


2 sin « sin B = cos (a — B) — cos (a =| 


2 sin acos B = sin (xn +B) + sin (a — B) J (16) 
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auf die Form 


uv 3 = 
a, 7, sae a0, =|K — aa “oe % K?| sin y t 
2 1 
ee (1) 
2(a,— 2 a,)2 K? sin 2y 7 + ee AG a3 sin 39 T 
bringen. Der Ansatz 
8, = Cysinvt + Cysin2y7 + Cysin 3 yt 
liefert die partikulare Lésung von (17): 
= K _ 3(a3— v? ay) + ¥? a5 x73]. Y Ae . 
° a, — 7? ay 4. (a, — v? a,)* os Pa 2 [a, — (2 v)? a;] (ag — v? a)? K? sin 2yc| 
18) 
| a,— »? (a, +45) Ques 
fire Nera K sin 397. 


Durch weitere Iterationen lieBe sich auf diesem Wege eine beliebig lange Fourierceihe fiir 
#(t) gewinnen. Wir wollen uns zur Bestimmung der optimalen Pendelmasse im folgenden mit 
der Naherung @, begniigen. 

Setzt man den so erhaltenen Ausdruck (18) in die Gleichung (12) ein, so geht diese mit den 
Identitaten (16), wenn man nur die Glieder bis einschlieBlich dritter Ordnung in K beriick- 
sichtigt, nach Integration nach t itiber in 
Aw Uifes— se; 


ee y 2 
jz i =r (cy v? cy) 


Bise ae elie Sat aorain Kel cos ve 


4 (a, — v? a,)* 4 (a, — v? a,)8 


[eg — (2)? c,] a 
| ie ae STG a,)? K? cos 2t (19) 


1_ {lez — (3)? ey] [ag — v? (ay + a5)] Cz — V" Cy 
12 3 [a, — (3 v)? ay] (a, — v7 a,)3 (a,— = Be 
Die Rechtfertigung fiir die Streichung der héheren Potenzen von K folgt sofort aus einer 
GréBenabschatzung fiir diesen Parameter an Hand von ausgefiihrten Motoren. K ergibt sich 
danach in der GréBenordnung von 10 bis 10~% je nach Tangentialdruckdiagramm, Kurbel- 
wellenkonstruktion und Drehzahl der Maschine. 

Die Gleichung (19) zeigt, daB sich das Glied erster GréBenordnung in K durch Erfiillen der 
bekannten Abstimmbedingung 


; 1—4=0, (4=7) (1a) 
tilgen 14Bt, da nach (12a) und (13a) 


2 ear} 
Co — ¥? cy 1—r?1 


a,—va,  (1—vrA)+ (1+ O 


ist. Weiter aber erweist sich, daB das Glied der zweiten GréBenordnung in K, welches von den 
Corioliskraften herrithrt, durch entsprechende Wahl der Pendelmasse m (d. h. des Parameters 
@) zu einem Minimum gemacht werden kann. Nach Erfiillung der Bedingung (la) bleibt 
namlich nach Einsetzen von a, bis a, bzw. c, bis c, nach (13a) bzw. (12a) in (19) fiir die Stérung 
Aw der Winkelgeschwindigkeit der Maschine: 
Aw y? 3 v8 
@  8@+1) 6 K'cos 1t +3734 1) Ol? +1)? O— 394] 
| 2 Moe? aT] 

8 (v2 +1)? 6? [(v? +1)? © — 8 v4] 


K2 cos 27 


(20) 
K° cos 3yT. 


Differentiiert man somit den Koeffizienten von K? cos 2 v7 in (20) nach dem Parameter O 
ae setzt diesen Differential quotienten gleich Null, so ergibt sich 


3 v8 (v2 +1)?- 20—304 (21) 


2 (v2? +1)? 67[(r? +1)? O—3r4?? 


und daraus fiir die optimale Pendelmasse wegen (10) die Beziehung 


__ Mopt L? 3 v4 99 
Oop = ai ~ 2(v? -+1)2° F ( ) 
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3 

Die Gleichung (22) liefert z.B. fiir »=1: m, = oS J , 
24 

fiir V2: Mont = 55 5 


Durch Einsetzen von (22) in (20) ergibt sich bis auf die Glieder mit vierten und| 
héheren Potenzen von K das so erreichbare Minimum der Winkelgeschwindigkeitsschwankung 
der Anordnung zu 


ee (+1)? 


2 y2+t] 
Do 27 ye K? cos » Wot — s 5 K* cos 2vapot T Tt [(3 v)2+ 1]- K3 cos 3y Wot. (23) 


In Abb.2 ist ohne Beriick-} 


‘Amplitude von ae 


F - ae 3 sichtigung der Glieder mit K® fiir 
| erfillte Abstimmbedingung (la) die 
a = I i Amplitude der dimensionslosen 
| Gx 6Be ( he ) als Funktion des Para-| 
Y alll 4 @, K 
2 
<ul 24 ed Lad % V=3 meters 90 = em fiir die Werte v=1, 
5 a 1-+}}-—-— at} 
. Rg | S| 2,3 aufgetragen. Es zeigt sich, daB| 
2 S41 al 2 N entgegen der bisher ttblichen Auf: 
Oppp | fr S Oppp fir Ss 1 S| fassung, modglichst groBe Pendel- 
p= 


massen zu verwenden, die jeweils) 

Pe ee ae doppelten Werte der optimalen 
7 Nee aid aa : : 5 Pende!l massen unerwiinschte Rese- 
; ae nanzen hervorrufen. Fir héhere 
Werte von v konvergiert O,,, gegen 
3/2, Onis gegen 3. 

5. Zusammenfassung. In der 
vorliegenden Arbeit wird die Kin- 
wirkung eines an einer Drehmasse 
angelenkten Fliehkrafttilgers, der 
als mathematisches Pendel auf- 
gefaBt wird, fiir groBe Pendelaus- 
schlage untersucht. Das besondere 
Interesse gilt dabei der Auffindung 


no einer Bestimmungsgleichung fiir die 
Abb. 2. Die Amplitude von a als Funktion des Parameters anzubringende ginstigste Pendel- 
Se re masse. Die Rechnung zeigt, daB sich 
@=- ; ‘ 3 : 
J Beziehungen fiir deren Bestimmung 


erst ergeben, wenn man den von) 
den Corioliskraften herrithrenden Term in den Bewegungsgleichungen des Systems beriick- 
sichtigt. Entgegen der bisherigen Auffassung, giinstigste Wirkung sei durch méglichst groBel 
Pendelmasse bei erfillter Abstimmbedingung L/l = y? zu erzielen, erweist die Rechnung die: 
Existenz von optimalen Werten fiir die anzubringende Pendelmasse und zeigt weiter, daft 
die doppelten dieser Werte zu unerwiinschten Resonanzen fiihren. | 
Die Erweiterung der vorstehenden Betrachtungen auf die Wirkung bei Anlenkung von meh- 
reren Tilgern an elastische zusammenhingende Ketten (als Reduktionsergebnis einer Kurbel- 
welle) soll in einer nachfolgenden Arbeit gezeigt werden. 


(Eingegangen am 10. November 1952.) 


Anschrift der Verfasser: Dr. K. Desoyer und Dozent Dr. A. Slibar, Wien IV, Karlsplatz 13, 
Technische Hochschule. 
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Zur Theorie der Photonenraketen. 
Von E. Sanger. 


1, Einleitung. Wenn man die technische Beherrschung stationarer Kernreaktionen vor- 

-aussetzt, begegnet deren direkte Ausniitzung in Raketen bei den Teilchengeschwindigkeiten 

“von etwa 4 bis 12% der Lichtgeschwindigkeit, entsprechend spezifischen Impulsen zwischen 

1220 und 3660 Sekundentonnen Antrieb je Kilogramm Treibstoff, der bekannten Schwierig- 

- keit, daB® sich die schnellen Reaktionsprodukte mit technischen Mitteln nicht in eine ge- 
wiinschte Richtung lenken lassen. 

Es wird daher oft erwogen, die Energie der Kerntreibstoffe auf weitere, inerte Massen zu 
iibertragen, die sich mechanisch oder elektrisch in die gewiinschte Richtung lenken lassen, und 
diese mit geringerer Geschwindigkeit auszustoBen, z.B. in thermischen Atomraketen mit spe- 
-zifischen Impulsen bis etwa 2 secton/kg oder in elektrischen Ionenraketen mit etwa 10 sec- 

ton/kg, gegeniiber den bekannten Werten chemischer Raketen von nur etwa 0,3 secton/kg. 

Daneben besteht jedoch noch der hier naéher untersuchte grundsatzliche, umgekehrte Weg, 
die Energie der Kernreaktion nur auf einen Teil der Kerntreibstoff-Massen zu iibertragen, und 
diesen Teil mit héherer, z. B. Lichtgeschwindigkeit auszustoBen, entsprechend spezifischen 

Impulsen von etwa 24 bis 216 secton/kg, wahrend der nichtbeschleunigte, verbrauchte Treib- 
| stoffteil geschwindigkeitslos von Bord gegeben, oder an Bord behalten wird. 

Fiir derartige partielle Photonenraketen erscheinen technische Lésungsméglichkeiten in ab- 
sehbarer Zukunft nicht vollig ausgeschlossen, z. B. wenn es gelingt, die Kernreaktionsenergie 
‘in ,,Kernlampen* nach Art elektrischer Gasentladungslampen weitgehend in Photonen zu 

-verwandeln und diese mit optischen Mitteln, etwa ,,Reflektoren“ héchsten Reflektionsvermé6- 
eens in die gewiinschte Richtung zu lenken. 

Sie versprechen gegeniiber thermischen Atomraketen geringeren spezifischen Treibstoff- 
_verbrauch und gegeniiber elektrischen Ionenraketen besonders geringere spezifische Trieb- 

werksgewichte, so daB sich mit ihnen Aussichten auf Fluggeschwindigkeiten von mehreren 
hundert Kilometern in der Sekunde eréffnen. 

Vom Vorgang der Beschleunigung einer gegebenen Ruhemasse m, unter Energiezufuhr, 
z. B. der Beschleunigung eines Elektrons in der Elektronenschleuder, ist bekannt, daB mit 
_Annaherung der Bewegungsgeschwindigkeit w an die Lichtgeschwindigkeit ¢ bei gleichblei- 
bender Ruhemasse m, sowohl die trage Masse m, = my/)/1—w?/c? als auch die kinetische Energie 
und der Impuls nach Unendlich gehen, also die Energiezufuhr schlieBlich unendlich gro wer- 
den miBte, so daB bei diesem Vorgang die Lichtgeschwindigkeit nicht erreicht werden kann. 

Wesentlich verschieden davon stellt sich der hier zu behandelnde Fall einer adiabatischen 

Beschleunigung der gegebenen Ruhemasse m, ohne duBere Energiezufuhr, also auf Kosten des 
Energieinhaltes der Ruhemasse selbst dar, wie er z. B. dem Raketenprinzip inharent ist. 
) Dabei nahert sich mit abnehmender Ruhemasse my die Bewegungsgeschwindigkeit w eben- 
falls der Lichtgeschwindigkeit c, um sie bei verschwindender Ruhemasse m  tatsachlich zu 
-erreichen, wahrend die trage Masse m, wahrend des ganzen Beschleunigungsvorganges kon- 
stant und gleich der Anfangsruhemasse my bleibt. 

Infolge des Aquivalenzprinzipes von Masse und Energie E = myc? gilt jederzeit die Adia- 

batenbedingung 
E= m, ce = m, c= mjc? + (my — m))c? = konst (1) 
Wir sprechen im Grenzfall mj — 0 von Zerstrahlung der Materie und nennen das entstan- 
dene lichtschnelle Partikel ein Photon, dem also wohl trige Masse, aber keine Ruhemasse mehr 
zukommt. 


2. Totalraketen. Man kann einer im Zustand der Ruhe gegebenen Masse my einen latenten 
Energieinhalt m)c? und einen latenten Impuls myc zuschreiben, ahnlich, wie man einem Photon 
der trigen Masse hy/c? die Energie hy und den Impuls h»/c beilegt. 

In der folgenden Raketentheorie spielen die Bruchteile E’ bzw. J’ des latenten Energie- 
bzw. Impulsinhaltes, die man technisch nutzbar machen kann, eine grundlegende Rolle. 


15 
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Wir wollen diese beiden Ausbeuten mit den Sigeln 


é= E'/m,c? fir die Energieausbeute (2) 
und 
t= J'/m,c fir die Impulsausbeute (3) 
bezeichnen. 
Der in der Raketentechnik ibliche spezifische Impuls ist dann 
J a= 0 (4) 
bzw. der spezifische Treibstoffverbrauch je Impuls- oder Antriebseinheit 
gd Leis (5) 


Aktiviert man einen Bruchteil ¢ dieser latenten Energie m,c? der Ruhemasse mp, entzieht 
ihn aber nicht dem System, sondern verwandelt ihn vollstandig in kinetische Energie der ge- 
samten restlichen Ruhemasse mj) = (1 — ¢) my, so wird die erreichbare Bewegungsgeschwin- 
digkeit wim Verhaltnis zur Lichtgeschwindigkeit ¢ aus der Adiabatenbedingung und dem he- 


kannten Verhiltnis )1 — w?/c? von Ruhemasse mj zu trager Masse my gleich 
oC — J 1—m?2/m? = 2 a= a. bzw. e= 1— l= wie (6) 
Die kinetische Energie ist gleich der aktivierten Energie 
E’= em,c?= (m—m,) c= myc? (1 — yi— w?/c?) : (7) 


Der gewonnene Impuls ist 


J’ = myw+ em, w= mw = myc /2e— &, (8) | 
daher die Impulsausbeute 
Osea Mi ati j2e—#= w/e , (8a) 
damit ist der spezifische Impuls 
J'/my= t+ e=cy2Ze— 2 = w (8b) 


und der Impuls je Einheit des Energieaufwandes 


i MeL p cae freee wie 
J'/E’ = y2/e— 1]/e = iiianyareaaae (8c) | 
Mit wachsender Energieausbeute ¢ wird daher der Impulsgewinn je aufgewendeter Masse immer 
gréBer, je aufgewendeter Energie immer kleiner. 

Diese Verhaltnisse der adiabatischen Be- 
schleunigung einer Masse sind in Abb.1 graphisch 
dargestellt. | 

Wahrend mit wachsender Energieausbeute ¢ 
die erreichte Strahlgeschwindigkeit w/c in einem | 
Kreisbogen zunachst rasch, spater langsamer 
wachst, und die beschleunigte, unverwandelte 
Ruhemasse m, linear mit ¢ abnimmt, bis sie bei 
kinelische Energie —| volliger Zerstrahlung schlieBlich verschwindet, 


| 
trage Masse mt/Mp=7 | _ 


Geschwindigkelt 10/6, 
=[mpulsausbheure b 
=Vee-e? 


NG 


ae Eng o*=€ bleibt die trage Masse m, dauernd gleich m, und 
: konstant. Dementsprechend wachst der ge- 
| aS wonnene Impuls wie w/c. 
a =. Im nichtrelativistischen Bereich kleiner Energie- 
BU As | ausbeuten € wachst der Impuls sehr rasch mit der 
resthiche Ruhemasse™. aufgewendeten Energie, der Antrieb je aufgewen- 
Mgitg Pe : Ree _ deter Energieeinheit ist giinstig, je aufgewendeter 
0 05 = Masseneinheit ungiinstig. Bei hohen ¢ kehren sich 
Abb. 1. Adiabatische Massenbeschleunigung diese Verhaltnisse um. 


in. Totalraketen. Schon relativ kleine Energieausbeuten ergeben 


verhaltnismaBig hohe Impulsausbeute, z. B. chemische Reaktionen mit ¢ = > w/c? ~~ 10% 


ein 6 /2€ = 1,4- 10~ oder spezifische Impulse von J’/m, —~- cy2e= w = 4,2 - 103 m/sec oder 
0,42 secton/kg. 
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; Kernreaktionen mit ¢ = 5 - 10-3 ergeben J’/m) = w = 3-107 m/sec oder 3000 secton/kg, 
} wahrend bei hundertprozentiger Umwandlung (¢ = 1) der spezifische Impuls auf 30500 sec- 
| ton/kg als Grenzwert anwachst. 

Dem ArbeitsprozeB der meisten derzeit bearbeiteten Raketentriebwerke liegt dieser spe- 
j zielle Vorgang adiabatischer Beschleunigung einer Masse mit totaler Energieiibertragung an 
idie Treibstoffmassen zugrunde. Wir haben diese adiabatischen Raketensysteme zur Unter- 
‘ scheidung von den im nachsten Abschnitt behandelten daher kurz Totalraketen genannt. 
Fir die mit ihnen erreichbare Fluggeschwindigkeit v im widerstands- und schwerefreien 
| Raum in Abhangigkeit von der relativen Treibstoffladung m)/M, gilt bekanntlich! die rela- 
‘ tivistische, sogenannte Raketengrundgleichung 


si pa Band Oh 

| My * iceaal (°) 
) worin M, die Anfangsruhemasse der Rakete und my ihre anfangliche Treibstoff-Ruhemasse ist. 
| Nach Einfithrung der Energieausbeute ¢ kann man sie auch in der Form schreiben 

= 1 


pies rey (9 


My 1+ vic 


) worin sich ¢ wie bisher auf das triebwerksfeste- Koordinatensystem bezieht. Im nichtrela- 
) tivistischen Bereich kleiner und mittlerer v und w geht sie in die klassische Raketengrund- 


| gleichung iiber 


ie — e— v/v 

a ie (10) 
> bzw. 

Se a (10a) 


3. Partialraketen. Es ist fiir raketentechnische Zwecke interessant, den bisher behandelten 

_speziellen Vorgang noch dahingehend zu verallgemeinern, da man die freigemachte Energie 

/em,c* nicht gleichmafig auf die totale restliche Ruhemasse mj)= (1 — €) m, aufteilt, sondern 

-einem Bruchteil ~ von mj keine Energie zukommen 1a4Bt, also nur den Rest (1 — 2) mj = 

(1 — a) (1 — €)m, auf eine Geschwindigkeit w diabatisch beschleunigt, wahrend der Bruchteil 

am, = % (1 — €) m, in Ruhe zum Triebwerk bleibt. 

a) Theorie der Partialraketen. Aus der Adiabatenbedingung, daB die Summe 

der tragen Massen bei dem Vorgang wieder konstant und gleich m, bleiben muB, folgt fir die 

-erreichbare Geschwindigkeit w folgende Uberlegung. 

Die Anfangsruhemasse m, verwandelt sich zu einem Teil ¢m,in Energie, wihrend der Rest- 
betrag mo = (1— e)m, iibrigbleibt. Von diesem Resthetrag bleibt der Anteil 7m) = 2(1—e)m, 
endgiilltig in Ruhe zum Triebwerk, wahrend der andere Teil (l1—a)mo = (1 — x) (l—e)m, 

mit Hilfe des verwandelten Energieadquivalentes ¢m, auf die Geschwindigkeit w gebracht wird, 

wobei seine gesamte triage Masse also (1—z) (l—«e)m)+em)= m, [1—a (1—e)] wird, wah- 

‘rend seine Ruhemasse natiirlich (1 — 2) (l—e)m, bleibt. 

Da die Ruhemasse sich zur tragen Masse immer wie y1— w/c? verhalt, folgt 


wade V2e(1—a) — e2(1— 27) (11) 
c 1—x (I—e) , 


Die kinetische Energie ist wieder gleich der aktivierten Energie: 


(1—z) (1—/1— w/e?) (12) 


Y1— w/c? + (1—a)(1—1— wife?) 


EE em, == (mi, — me = mi, 


Der gewonnene Impuls ist 


J'= (l—2) myw+ em w= mw[l—a(1—e)J= m,cy2e(1—x)— 2 (1— 22), (13) 


daher die Impulsausheute 
t= J'm,c= 26 — a) — #1 — 2m) = [1—a (1—«)] w]e (14) 


1 J, Ackeret, Helv. Physica Acta 19 (1946), S. 103. 
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und der massenspezifische Impuls 
J'|m,= c+ y2e(1— a) —e& (1— 22) = [l—a (1— 8) Jw, (15) 
schlieBlich der energiespezifische Impuls 


pee Pel aye? (S22) 1 w (1—a(1—e)][VI—w2/e? + (1 —x) (1— V1 — w/c?) ] 16 
T= — ; lees (16) 
ec coe (1— x) (1—J1— w/e? ) 

Die Impulsausbeute ¢ und damit auch der massenspezifische Impuls sind in Abb. 2 als 
Funktion der Energieausbeute ¢ und der Partialitat ~ dargestellt. In den Grenzfallen a = 0 
und wz = 1 verhalten sie sich wie 

EA (We) ee (17) 


t(% =1) 


d.h. in allen von ¢ = 1 verschiedenen Fallen hat die Photonenrakete (7 = 1) schlechtere Im- 
pulsausbeute als die Totalrakete (7= 0). 

In Abb. 2 bedeutet namlich die Partialitat 7 — 0 wieder den schon eingangs behandelten 
Fall der Total-Raketen, wo die freigesetzte Energie zur Beschleunigung der gesamten ver- 
bleibenden Masse verwendet wird. Dieser Fall gibt 


Sale aa bei gleicher Energieausbeute immer die héchste Im- 
ie Me Photonen—Rokesen pulsausbeute. 

SS N By Partialitat a —=1 bedeutet, da® nur die frei- 

y & z gesetzte Energie in Form von Photonen ausgesendet 

\) | wird, wahrend die gesamte, nicht umgewandelte Treib- 


R Photonenraketen ist immer am kleinsten. 

05 = -Partia/-Rakelens Die Unterschiede der Impulsausbeute sind natur- 
gemaB am gréBten, wenn die Energieausbeute ¢ 
klein bleibt, da dann die gréBten Massenbetrage fiir 
den Impuls unausgewertet bleiben, die Unterschiede 


7 werden dementsprechend immer geringer, je gréBer 
 € wird. 
Fnergre- ? Dae 
| cusheuk e Bei ¢ = 1, also vélliger Massenzerstrahlung, hat 
0 05 40 die Impulsausbeute ihr absolutes Maximum Eins, 


Abb. 2. Impulsausbeute , der Partialraketen in Ab- alle Werte von a beziehen sich hier auf die gleiche, 
haingigkeit von der Energieausbeute ¢ und der Par- : ° 
tialitat z. (Der massenspezifische Impuls ergibt sich verschwundene Restmasse Null, und sind unter sich 
aus ¢ durch Repl paketion aa der Lichtgeschwindig- also gleichwertig, es liegt der Fall der totalen Pho- 
elt c. 
tonenrakete vor. 
Solange wir in der Raketentechnik nur die geringen Energieausbeuten chemischer Reak- 


tionen zur Verfiigung haben, ist das Interesse, diese geringen Energien der gesamten Treib- 


stoffmasse mitzuteilen, daher gro8. Die Impulsausbeute verhalt sich hier in den Fallen | 


a= 0 und w= 1 etwa wie y2/e. 

Mit dem ¢ = 107! chemischer Reaktionen fanden wir fiir 7 = 0 den massenspezifischen 
Impuls zu etwa 0,42 secton/kg, als Photonenrakete mit 7 = 1 betriige er in diesem Fall nur 
etwa ein Hunderttausendstel, namlich 3 - 10~* secton/kg. Die Photonenrakete auf chemischer 
Treibstoffbasis ist daher wohl als Scheinwerfer, nicht aber als Triebwerk von Interesse. 

Mit dem ¢ = 7- 107% der Wasserstoff-Helium-Reaktion wird der massenspezifische Im- 
puls fiir 7 = 0 nach (8b) etwa 3,55 - 107 m/sec oder 3550 secton/kg ; fiir 7 = 1 nach (15) etwa 
2,1 - 10° m/sec oder 210 secton/kg, d. h. die Wasserstoff—Helium-Photonenrakete hat nurmehr 
etwa siebzehnmal geringere Impulsausbeute, als die totale Wasserstoff-Helium-Rakete. Sie 
kénnte daher durchaus praktisches Interesse bieten, wenn ihre Verwirklichung geringere tech- 
nische Schwierigkeiten bereitet, als jene der entsprechenden Totalrakete. 

b) Partialraketen mit Inertmassen. In den Bereich der Partialraketen 
gehéren von vornherein auch alle Raketen-Vorschlage, bei denen inerte Massen My, durch 
Energietrager mo, mit der Ausbeute ¢, auf thermischem oder elektrischem Wege beschleunigt 
werden, wahrend die nichtumgewandelten Energietragermassen nicht beschleunigt werden, 


z. B. thermische Atomraketen oder elektrische Ionenraketen 12345 wo die — etwa im Uran-_ 


1 J. Ackeret, a.a. O. 2 H. Oberth, Wege zur Raumschiffahrt, Miinchen 1929 und Ann Arbour 1945. 
8 L.R.Shepherd u. A.V.Cleuver, J. Brit. Interpl. Soc. 8 (1949) S. 59. 
4 J. Bredt, Naturforschg. 6a, 2 (1951) S.103. 5 L.Spitzer, J. Brit. Interpl. Soc. 10 (1951) S, 249. 


stoffmasse in Ruhe bleibt. Die Impulsausbeute dieser | 
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meiler gewonnenen — Energien zur entweder thermischen Beschleunigung von Wasserstoff- 
plasma auf Geschwindigkeiten von der GréSenordnung 2-104 m/sec, oder elektrischen 
| Beschleunigung von Ionen auf Geschwindigkeiten zwischen 105 und 106 m/sec dienen sollen. 
Mit einer Energieausbeute ¢, des Energiestoffes my,im Meiler und der zusatzlichen Plasmen- 
oder Jonen-Ruhemasse m,, wird die gesamte Energieausbeute 


1 
ious 1 + mo3/Mo1 


(18) 
und die Partialitat 
1—e 
t= - ag 9 (19) 


| womit alle weiteren Gréen nach den Gleichungen (11) bis (17) berechenbar sind. Die darin 
immer wiederkehrenden Ausdriicke werden 


a= 2 e(L—m) — FB) = |) Seal) alr) Frain) (90) 


(1 + mo9/1 91)? (1 —e, + mg/m) 


| 
7] TT oY) y 3 
fe : = : : iV . 10 10 
& Cc 
on 7 BS ar i < 15 = 
z ar 
ee lsat | SS ee ie 
eal y S NI sia SS “ 
NI N Si) = ' 
x 8 N Slay ss 3 0 
L aN ~ n = 
ES : t 8 10° S 41007? 40 
FS ‘S xs 3] 8 N 
=a ie Ss ye SN z SS 
ies = <2] 9) & Ss 2) 
2 SN & SS & S 
iS t st S| 8 N g 
x S See S 
) ls ant? Scie cde 8 
x xX & S 
as { ut 0’ & 40? § sve 
ms = s 2 S 
/ S ‘ = ie) 8 
> S {—~ § x 8 s&s 
tJ /mg=EC +f | 8S 3 $ 
J ee eee lees 
+ © C | S S S 
S & 3S) 
s 
E si i =: a 0° 70’ 70°§ 
{ : {| ~2 
a i \ fee 
ab \ Q 
~N 
7 & 
| Partiohtit- | s 
sli Sa? 4 
4— SS 
=a 8 
& 
f g 
1 Se 
0 (| Til 
wo? w? wo 0% ww? wt 0° 0 0 0? = 10% §— 10% roe _ inerteMasse 


Tor — Mejlermasse 


Abb. 3. Partial-Raketen mit Inertmassen; Verinderlichkeit iiber my /m,, der Auspuffgeschwindigkeit und des massenspezifischen 
und energiespezifischen Impulses, ferner der durchschnittlichen Impulsausbeute, fiir Atom-, Ionen-, und Photonen-Raketen mit 
Meilerheizung yon ec, = 10°* Energieausbeute. 


und 


kel x (1 e)= &1 + Mo9/™Mo1 (21) 


1 + mpo./™o1 


und damit 


w/e= a/b, (22) 
foo 4. e=— Ow, (23) 
J'/E'= alec= (w/c) b/ec. (24) 


In Abb. 3 sind diese GidBen iiber my/m,, fiir ¢,= 10-8, wie es etwa einem leichten Uranmeiler 


-entsprechen kénnte, aufgetragen. 
Der Grenzfall mj./m),= 0 bedeutet die Photonenrakete mit 


&=s; t= 1; 2=e; b= a3 wle=—1; J']/my= ee; J'/E' = le. 
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Der massenspezifische Impuls und der Wert a haben ein Maximum, wenn die beschleunigte 
Masse myp+ €,m, gleich ist der Meilermasse mp), also bei 


Mpo/™Mo, = 1 — &,, (25) 


und zwar von der Héhe 


ee Ve (26) 


2— & 


Dieser Fall hat offenbar gleiche Antriebseigenschaften wie jener der Totalrakete, nur daf hier | 
nicht die unumgewandelte Meilermasse (1 — ¢,) m,, sondern die gleich groBe Inertmasse mp, be- 
schleunigt wird. Dieses Maximum ist daher 1//2e,—¢?-mal gréBer, als die Impulsaus- 
beute der entsprechenden Photonenrakete gleicher Energieausbeute ¢,. Der Verwirklichung 
dieser giinstigsten Partialrakete stehen jedoch dieselben praktischen Schwierigkeiten ent- 
gegen, wie der entsprechenden totalen Atomrakete. Eine der Ursachen dieser Schwierigkeiten 
ist der in Abb. 3 ebenfalls als Energieaufwand je Antriebseinheit aufgetragene Kehrwert des 
energiespezifischen Impulses J’/E’. Mit der Gréfe dieses Energieaufwandes wachsen bei vielen 
Bauarten die unvermeidlichen Warmeiiberginge an das Triebwerk selbst, wodurch die 
thermischen Atomraketen gegenwartig auf den eingezeichneten Bereich der Auspuff- 
geschwindigkeiten um w= 20000 m/sec beschrankt erscheinen, wahrend man bei elektrischen 
Ionenraketen durch Vermeidung hoher Arbeitstemperaturen glaubt bis an den in Abb. 3 
gleichfalls eingezeichneten Bereich der Auspuffgeschwindigkeiten von 10° bis 10° m/sec gehen 
zu kénnen. 

Die Photonenrakete gleicher Energieausbeute ¢, = 10-3 hat etwa adhnlichen massenspezi- | 
fischen Impuls, wie die genannten Jonenraketen, jedoch bei mehrtausendfach héherem Ener- 
gieaufwand je Impulseinheit. Sie ist also gegen die Ionenrakete nur dann wettbewerbsfahig, 
wenn der zu befiirchtende Warmeiibergang aus Ursachen, die in ihrer spezifischen Arbeits- 
weise liegen, wesentlich kleiner wird, als bei Ionenraketen. 

In Abb. 3 ist noch bemerkenswert, daB im Bereich groBer mo,/m,, die Auspuffgeschwindig- 
keit w, der massenspezifische Impuls w und der Energieaufwand je Impulseinheit w/2 parallel 
verlaufen, wie wir dies von nichtrelativistischen Totalraketen gewéhnt sind. In der Tat 
sind die typischen Eigenschaften der Partialraketen auf den Bereich der mo/m), zwischen 
etwa 10-4 und 10*? beschrankt, unterhalb ist die Wirkungsweise schon praktisch die der 
partiellen Photonenraketen, oberhalb die von Totalraketen, wie auch die a-Werte iiber 
Mp./Mp, in Abb. 3 zeigen. 

Es ist fiir die hier betrachtete zweite, diabatische Art von Raketen, die Partialraketen, 
wichtig, festzuhalten, daB die Auspuffgeschwindigkeit w selbst kein Ma8B mehr fiir den massen- 
spezifischen Impuls ist, wie auch der Unterschied der Gleichungen (8b) und (15) zeigt. Aus der 
letzten folgt beispielsweise, daB eine Photonenrakete mit 7 = 1 und w/c = 1 den spezifischen 
Impuls 


Jim EC (27) | 


hat, der also bei konstanter Auspuffgeschwindigkeit w= c nur von der Energieausbeute ¢ _ 
der Treibstoffe abhangt. 

c) Flugleistungen mit Partialraketen. Bei diesen Raketen zweiter Art, 
die wir kurz Partialraketen nannten, da sie einen Teil ihrer verbrauchten Treibstoffe nicht ab- 
stoBen, erhebt sich das Problem, was mit dem nichtausgestoBenen, verbrauchten Treibstoff- 
rest geschieht. Dies ist zwar fiir den bisher untersuchten spezifischen Impuls und dessen Kehr- 
wert, den spezifischen Treibstoffverbrauch, gleichgiiltig; wenn sie jedoch an Bord belassen 
werden, belasten sie das spezifische Triebwerksgewicht und beeinflussen damit die bei ge- 
gebenem Treibstoffvorrat erzielbaren Flugleistungen, bzw. die mit dem Raketenfahrzeug er- 
reichbare Endgeschwindigkeit v im ungiinstigen Sinn. 

Werden die verbrauchten, aber nicht beschleunigten Treibstoffmassen x (1 — ¢)m, dauernd 
an Bord behalten, dann kann man fiir die erreichbare Fluggeschwindigkeit v die Raketen- 
grundgleichung (9) fiir Totalraketen beniitzen, wenn man z (1 — e)m, iiberhaupt nicht zu den 
Treibstoffmassen m, zahlt und fiir w den Wert aus (11) einfiihrt. 

Werden jedoch die verbrauchten, aber nicht beschleunigten Treibstoffmassen ohne Relativ- 
geschwindigkeit zum Fahrzeug und kontinuierlich von Bord gegeben, so andert sich die rela- 
tivistische Raketengrundgleichung, da die Auspuffgeschwindigkeit nicht mehr mit der vollen 
Treibstoffmasse, sondern mit der um den geschwindigkeitslos abgetrennten Teil zt (1—e)m, 
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verringerten Masse [1 — a (1—e) my] wirkt, in 


My __ 1—v/c ae aaa 
M, : ( 1 =e) ‘ (28) 


‘wie man auch durch eine der Ackeretschen Herleitung der Gleichung (9) analoge Ableitung 
finden kann. 

| Im nichtrelativistischen Bereich kleiner und mittlerer v und w degeneriert diese Beziehung 
zu der klassischen Grundgleichung der Partialraketen 


Le — vo/w (l1— ax(l1—~e 
“poise = Ge) (29) 
Die Anfangsruhemasse M, des Fahrzeuges umfaft: 

1. die Treibstoffmasse m,, die ihrerseits bei Partialraketen zerfallt in 

a) den zu beschleunigenden Anteil [1 — 2 (1 — e)]m), 

b) den nicht zu beschleunigenden Anteil 2 (1 — ¢)m); 


2. die Triebwerksmasse mp, (Eigentliches Triebwerk, Tankanlagen, Triebwerksnebeneinrich- 
tungen usw.); 


3. die restliche Fahrzeugmasse my, (Zelle, Nutzlast, Ausriistung usw.); 
Im Falle des an Bord behaltenen Anteiles x (1 — ¢€) m, gilt daher Gleichung (9) in der Form 


My) 1 San Lk —le 2w a _ 

M, — 1-+ mg/m, + mof/m, r cena fi m (1 e)). (30) 
Fir w/c kann man je nach dem vorliegenden Fall die Gleichungen (6), (11) oder (22) einfithren. 
Wiinscht man v/c als Funktion des massenspezifischen Impulses oder des spezifischen Treib- 
stoffverbrauches auszudriicken, so kénnen dazu die Gleichungen (8b) baw. (15) oder (23) 


dienen. 


' 
; 


——=_ = SS Srl tC hc Ll 


Das spezifische Triebwerksgewicht ( eee ist in (30) ebenfalls enthalten. Da dieses 
'Triebwerksgewicht je Impulseinheit bei den hier behandelten Raketen auferordentlich von 
der angewendeten Fahrzeugbeschleunigung abhangt und daher weniger wie eine technische 
Triebwerkskennzah] beniitzt werden kann, als das Triebwerksgewicht je Schubeinheit, werden 
wir davon hier keinen weiteren Garewach machen. 


| Im Falle des geschwindigkeitslos von Bord gegebenen Anteiles 7(1 — ¢)m, gilt dagegen (28) 
-unmittelbar: 

| 1 

Tigh - 1 irre) 2/26 (1—a) — e (1 — 22) (31) 


M, 1-+mgi/m,+ mim, I +o/e 


: 


| wofiir dieselben Zusatzbemerkungen wie fiir (30) gelten. 

| In beiden Fallen 1aBt sich die im schwerefreien und widerstandsfreien Raum erreichbare 
'Fluggeschwindigkeit v/c als Funktion von m,/M,, ¢ und x darstellen. Dies ist in Abb. 4 ge- 
‘schehen fiir den Gesamtbereich der Totalraketen, Partialraketen und Partialraketen mit 
‘Inertmassen, so da dieses Diagramm also chemische Raketen, thermische Atomraketen, elek- 
trische Ionenraketen und Photonenraketen umfaft. Dabei wurden folgende Beziehungen be- 
| niitzt : 


ohne Massenauswurf, im relativistischen Bereich: 


2 V2e(1— =) — e*(1— 272) 


if ™o \ 1—a(1—e) 

Pals fe M [1—z(1 er (32) 
a ; 2 2e(1— a) —e (1— 22) 
EO — Lae : 1—a(1—e) 
1441 M, [1—x(1 ely 
ohne Massenauswurf, im nichtrelativistischen Bereich 

» 22 (1—a)—e(1-20) J) oy 1 33 
= Pe In | ™M, [ Tu ( e)]) > (33) 
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mit Massenauswurf, im relativistischen Bereich 


2/2 e(1 —x) — 2 (1—2 2) 
10-3) 
0 
cS m, \2 /2eG—™)—# (1-22) (34) 
0 
1+( ar 


mit Massenauswurf im nichtrelativistischen Bereich 


2=— y2e1—m— A 2aa)Im (1— 5p) - (35)| 
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Abb, 4. Erreichbare charakteristische Fluggeschwindigkeiten o/c mit einstufigen Total- und Partial-Raketen der Treibstoffladung 
m,/M,, tiber der Energieausbeute e. 


Fir 7=0 (Totalraketen) geben die Gleichungen (32) und (34) bzw. (33) und (35) natiirlich 


jeweils gleiche Ergebnisse : 


Totalraketen, relativistisch (entsprechend (9a)) 


SE elec (36) 
¢ eae” 
1+ (1-72) 
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Totalraketen, klassisch (entsprechend (10a)) 

v ey ae ss m 

+= 2e— e* In iLe-wre) : (37) 
Fiir den hier besonders interessanten Fall x = 1, ¢ + 1 (partielle Photonenraketen) vereinfachen 
sich die vier Gleichungen folgendermaBen: 
partielle Photonenraketen, ohne Massenauswurf, relativistiseh 


v 1—(]1—em,/M,)? 


e601 +(1—e mM)? ” (37) 


partielle Photonenraketen, ohne Massenauswurf, klassisch 

v 

~=— In (l1—em,/M,). (38) 
partielle Photonenraketen, mit Massenauswurf, relativistisch 


vw 1—(1—m,/M,)** (39) 
Oe Te Wie aa Ae ae 


partielle Photonenraketen, mit Massenauswurf, klassisch 
v 
or ae éeln (l—m,/M,), (40) 


Fir = 1, ¢ =1 (Totale Photonenraketen) werden die Resultate beider Gleichungsgruppen 
wieder gleich und hier unabhangig von z: | 


totale Photonenrakete, relativistisch 


vt (1—m/ My 
= 1+ (1—m,/M,)? ’ ep) 
totale Photonenrakete, klassisch 
= —In (l—m)/M,) . (42) 


Einige dieser Beziehungen sind in Abb.4 zahlenmafig dargestellt, so daB sich die mit 
einstufigen Total- und Partialraketen erreichbaren charakteristischen Fluggeschwindigkeiten 
v in Bruchteilen der Lichtgeschwindigkeit c direkt ablesen lassen. 

Die charakteristische Fluggeschwindigkeit v ist die bei beliebigem Beschleunigungsvorgang 
im widerstands- und schwerefreien Raum unmittelbar erreichbare Fahrzeuggeschwindigkeit. 

Jeder praktischen Flugaufgabe kann man bekanntlich* eine solche charakteristische Ge- 
schwindigkeit zuordnen, die nicht tatsdchlich erreicht zu werden braucht, die aber die not- 
wendige Treibstoffmasse zur Erfiillung des Fahrtauftrages kennzeichnet. Beispielsweise be- 
tragt sie fiir einen Uberschallflug um den halben Erdumfang etwa v/e~2-10-°, fiir 
den Aufstieg zu einer ErdauBenstation v/c~3-10~, fiir das Verlassen des Erdschwerefeldes 
vje~4,5-10-%, fiir eine Fahrt von der Erdoberflache zur Marsoberflache und zuriick etwa 
v/e =1,7-:10-*, zum Verlassen des Schwerefeldes der Sonne rund v/¢e ~1,3-10%. 

Als Abszissen sind die Energieausbeutene aufgetragen, iiber den Bereichen der chemischen 
Reaktionen, der Gemische von Kernreaktionsprodukten mit Inertmassen, der reinen Kern- 
reaktionen, bis zur vollstandigen Massenzerstrahlung, also von ¢=—10~1° bis «=; 

Von den Partialitaten z sind, um das Bild nicht zu iiberladen, nur die beiden Grenz- 
falle 7 =0 (Totalraketen) und 7=1 (Photonenraketen), letztere mit und ohne Massenaus- 
wurf, in Form von drei Kurvenscharen eingezeichnet, deren einzelne Kurven jeweils die prak- 
tisch interessierenden Verhaltnisse m)/M, von Treibstoffmasse zu Anfangs-Fahrzeugmasse 
enthalten. 

Man erkennt sofort wieder, daB im ¢-Bereich chemischer Reaktionen die Photonenraketen 
véllig uninteressant sind und praktisch nur Totalraketen cine Rolle spielen kénnen, wie dies 
auch wirklich der Fall ist. 

Die Gemische von Kernreaktionsenergie mit Inertmassen, wie sie beispielsweise von 
thermischen Atomraketen und elektrischen Ionenraketen verwendet werden sollen, haben 


1 A, C. Clarke, Interplanetary Flight, London 1950. 
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nach (18) durchschnittliche Energieausbeuten zwischen denen der chemischen und Kern- 
Reaktionen, wahrend diese Raketensysteme praktisch immer (siehe Abb.3) als Totalraketen 
betrachtet werden kénnen (7 =0), so da ihre Flugleistungen aus deren Kurvenschar abzu- 
lesen sind. 

Wenn man bedenkt, da8 bei einstufigen, thermischen Atomraketen Treibstoffladungen 
von m,|M, =0.8 vorstellbar sind, bei elektrischen Ionenraketen wegen der auBerordentlichen 
Triebwerksgewichte wahrscheinlich nur wesentlich kleinere Treibstoffladungen von viel- 
leicht m,/M,=0,2, dann erscheint selbst bei sehr schwach beschleunigten Fahrzeugen die 
Uberlegenheit der Ionenrakete iiher die thermische Atomrakete keineswegs sicher. Dazu muB 
noch beachtet werden, da® die thermische Atomrakete sowohl fiir starkbeschleunigte (Luft- 
fahrt) als auch fiir schwachbeschleunigte (Interplanetare Raumfahrt) Fahrten geeignet er- 
scheint, wihrend die Ionenrakete sich sicher nur fiir letztere eignen kénnte. 

Anderseits diirfte die vorausgesetzte Auspuffgeschwindigkeit von w —20000m/sec fiir 
thermische Wasserstoff-Atomraketen ebenfalls eine allerauBerste Grenze darstellen, die durch 
die starke Ionisation des Wasserstoffplasmas bei sehr hohen Temperaturen und das damit 
unter Eins sinkende Atomgewicht des Plasmas begiinstigt wird. Die Vorstufe dieser starken 
Tonisation ist aber eine sehr hohe Elektronenanregung der Wasserstoffatome, d.h. sehr inten- 
sive Strahlung des Plasmas, die zu sehr hohen Warmeiibergingen an die Ofenwande fithren 
kann. Eine Folge der starken Ionisation gerade bei Wasserstoff ist anderseits die rasche 
Zunahme der freien Weglangen des Protonen-Elektronen-Gases, die allmahlich mit den 
Ofenabmessungen vergleichbar werden und dann zu den sehr ungiinstigen gaskinetischen 
Strémungs- und Konvektions-Verhaltnissen fiihren!?*. 

In den Bereich der Energieausbeuten von Kernreaktionen kann man mit Totalraketen 
aus den bekannten Griinden nicht vorstoBen, méglicherweise aber mit partiellen Photonen- 
raketen, die nach Abb.4 sowohl den thermischen Atomraketen, als auch den elektrischen 
Tonenraketen sehr tiberlegen sein kénnten, trotz der aus der Abbildung deutlich erkennbaren 
erundsatzlichen Unterlegenheit der Photonenraketen gegeniiber den Totalraketen bei gleicher 
Energieausbeute, wie die Kurvenschar der Partialraketen mit 7=1 zeigt. 

Mit der partiellen Photonenrakete gelangen charakteristische Fluggeschwindigkeiten von 
v/e = 10-8 oder v = 300000 m/sec in den technischen Méglichkeitsbereich. 

SchlieBlich ist zu Vergleichszwecken in Abb.4 noch eine dritte Kurvenschar der partiellen 
Photonenraketen (7 =1) ohne Massenauswurf eingezeichnet. 

Die verbrauchten, aber an Bord behaltenen Treibstoffmassen machen sich natiirlich in 
der erreichten Endfluggeschwindigkeit umso starker bemerkbar, je gréBer die Treibstoff- 
zuladung m,/M, an sich ist, und je geringer deren Umwandlungsgrad « bleibt. 

Die Linien gleicher Treibstoffladung m,/M, jeder Kurvenschar konvergieren naturgema} 
mit wachsender Energieausbeute ¢, bis schlieBlich alle Raketensysteme in die totale Photonen- 
rakete mit vélliger Treibmassenzerstrahlung einmiinden, deren charakteristische Fluggeschwin- 
digkeiten sich auch bei maBigen Treibstoffladungen der Lichtgeschwindigkeit nahern. 

Weitere Kurvenscharen fiir beliebige a-Werte zwischen Null und Eins lassen sich mit 
Hilfe der Gleichungen (32) bis (42) ohne weiteres in Abb.4 eintragen. 


4, Partielle Photonenraketen. Die Photonenrakete als Spezialfall 7 =1 der Partial- 
raketen, mit dem massenspezifischen Impuls J’/m) = ec, dem energiespezifischen Impuls 
J'/E’ =1/e und der Auspuffgeschwindigkeit wc kann gegeniiber der Totalrakete gleicher 
Energieausbeute ¢ dadurch interessant werden, daB sich die Photonen leichter in eine ge- 
wiinschte Richtung lenken lassen, als die korpuskularen Produkte der Atomreaktion selbst. 

Abb.5 gibt vorerst die bekannte Ubersicht itber das Gesamtgebiet der elektromagnetischen 
Strahlung, geordnet nach Wellenlangen in cm, also ttber die Arten der Photonen und ihre 
wesentlichsten Higenschaften. 

a) Reflektor. Die hier zunachst interessierende Lenkung von Photonen in eine ge- 
wiinschte Richtung fallt in das Gebiet der Optik, gegebenenfalls der Réntgen- und Gamma- 
Optik, wobei die optischen Reflexionskoeffizienten in unserem Zusammenhang eine aus- 
schlaggebende Rolle spielen. 


1 E. Sanger, Gaskinetik sehr groBer Fluggeschwindigkeiten, Deutsche Luftfahrtforschung F. B. 972, 
Berlin 1938, Douglas Aircraft Comp., Transl. FR 369, Santa Monica 1946, N.A.C.A. Techn, Mem 
1270, Washington 1950. } 

* FE. Sdnger, Schweizer Arch. angew. Wissensch. Technik 16 (1950), S. 43. 

3 KE, Sanger, J. Métaux 22 (1952), S. 44. 
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Bei chemischen Raketen betragt der Warmeiibergang an die feuerberiihrten Wande rund 
ein Hundertstel der umgesetzten Energien, also etwa 100 keal/secton, der in der Regel von 
den Treibstoffen selbst, vor deren Einbringung in den Ofen, aufgenommen wird. 

Bei den Projekten der Ionenraketen sehr geringer Beschleunigung wird mitunter ange- 
nommen, dafs bis zu ?/; der umgesetzten Energie als Warmeiibergang an das Fahrzeug zu- 


lassig und durch Strahlungsflachen nach augen abgebbar waren, d.h. iiber 240000 keal/secton, 


\ 


secton. 


bei einem niitzlichen Energieaufwand je Antriebseinheit nach Abb.3 von etwa 120000 keal/ 


Wenn dieser letzte Wert sich bei Photonraketen nach Abb.3 um 7.108 keal/secton bewegt, 
also fast 6000fach hoher ist, und bei schwach beschleunigten Photonenraketen derselbe ab- 


solute Warmeiibergang je Antriebseinheit wie bei Ionenraketen als zulassig angesehen wird 


2 


dann diirften bei der Photonenrakete also etwa 0,34 Promille der umgesetzten Energie an das 
Fahrzeug, d.h. also praktisch an die Photonenoptik iibergehen. 
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Abb. 5. Photonenarten, Erzeugung, Lenkung, Frequenz », Energie hv, Masse hy/c? und Impuls hy/e iitber der Wellenlinge 1 = c/p. 


Sobald die optischen Reflexionskoeffizienten also héher als 0,9997 werden, ist in dieser 
Hinsicht Wettbewerbsfahigkeit mit der Ionenrakete méglich, wenn man an ein gekiihltes 
Spiegelsystem nach Art von Scheinwerfern denkt. Bei Totalreflexion nahert sich der optische 
Reflexionskoeffizient tatsachlich dem Wert Eins, es wird dem zweiten Medium keine Energie 
zugefiihrt. Es ist dies der in der Grenzflachenphysik einzigartige Fall eines Impulsaustausches 
fast ohne Energieaustausch, mit Reflex'onsgraden der GréSenordnung (1—10~%). Bei senk- 
rechtem Strahleinfall auf glatte Metalloberflachen nimmt das Reflexionsvermégen mit ab- 
nehmender Wellenlange von Werten nahe an Eins im Ultrarot und teilweise auch noch im 
Sichtbaren bis zu kleinen Werten im Ultraviolett ab. An Metallspiegeln mit Verstarker- 
lamellen hat man durch Ausniitzung der entstehenden Interferenz und des Phasensprunges in 
der Lamelle Vorderflachen-Silberspiegel mit einem Reflexionsgrad von 0,995 fiir weiBes Licht 
erreicht, wahrend normale Oberflachen-Silberspiegel nur bis 0,95 reflektieren. 

Im Réntgen- und Gamma-Gebiet sind die Verhaltnisse wesentlich ungiinstiger, dort 
bleiben fiir Photonenlenkung nur die Braggsche Kristallreflexion und die Totalreflexion tibrig. 
Bei Kristallreflexion im Gamma-Gebiet werden erfahrungsgemaf weniger als 1% der ein- 


fallenden Intensitat reflektiert. 
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Totalreflexion derart kurzer Wellen tritt nur bei sehr streifendem Kinfall auf, da die 
Brechungszahlen der kurzwelligen Photonen nur sehr wenig von Eins abweichen. Es scheinen 
sich also am ehesten im langwelligeren Bereich Méglichkeiten fiir die erforderlichen hohen 
Reflexionszahlen zu eréffnen. Dabei wird noch von Bedeutung sein, daB die hohe Reflexion 
bei ganz auBerordentlich hohen Strahlungsintensitaten und gegebenenfalls bei hohen Tem- 
peraturen der reflektierenden Oberflichen erreicht werden muf. 

Hier tut sich ein weitgehend neuer Forschungsbereich der physikalischen Optik auf, 
Welche Reflexionszahlen man aber auch immer erreichen wird, geht ein kleiner Teil der Energie 
an das Reflexionssystem iiber, der auf geeignete Weise aus dem Fahrzeug entfernt werden mul. 
Man kénnte beispielsweise an eine stufenformige Reflexion der Art denken, da8B die Ursprungs- 
intensitét I, des Photonengenerators mit der Reflexionszahl r,; von einem Spiegel zurtick- 
geworfen wird, der die aufgenommene Intensitat (1—r,) J, von einer a-fach gréBeren Flache 
der Temperatur T einem zweiten Spiegel schwarz zustrahlt, welch letzterer diese Restinten- 
sitit in einer zweiten Stufe in die gewiinschte Richtung reflektiert, falls die Gleichrichtung 
dieser Abwirme noch fiir lohnend gehalten wird. Zwischen den einzéelnen GriBen besteht 
offenbar die Beziehung - 

(1—r,) =aoT’, 
worin ¢ =5,77-:10-12 Wattem*grad~4 die Strahlzahl im Stefan- Boltzmannschen Strahlungs- 
gesetz schwarzer Kérper ist. 

Fiir Photonenraketen interessieren Ursprungsintensitaten bis in die GréSenordnung 
1019 Watt/em? =2,4-108 keal/em’sec. Bei den praktisch méglichen Strahlungstemperaturen T 
der Kithlflachen wird man meist zu betrachtlichen Werten fitr « kommen, die den bekannten 
Kithlvorschlagen fiir Ionenraketen entsprechen. 

Mit « =1 wiirde man besondere Einrichtungen zum Warmetransport ersparen und kénnte 
mit T daher bis an die Grenzen der Gliihtemperaturen fester Kérper gehen, doch setzt dies 
wieder sehr kleine (1—r,) voraus, iiber deren Realisierungsméglichkeiten noch nicht geniigend 
Kenntnisse vorliegen. 

Wahrend also fiir das unmittelbare Gleichrichten der Reaktionsprodukte von Kernreak- 
tionen zunachst iiberhaupt keine technischen Lésungen zu erkennen sind, und fiir das Gleich- 
richten solcher Impulse iiber den Umweg der Ionenbeschleunigung durch elektrische Felder 
mit ?/, Energieverlust an das Fahrzeug gerechnet wird, kénnte ein Gleichrichtvorgang iiber 
den Umweg der Photonen trotz seiner prinzipiellen Energievergeudung technische Vorteile 
bieten, u. U. sogar geringere Warmeiibergange an das Fahrzeug je Antriebseinheit ergeben. 

Jedenfalls erweist sich, daB das erste Teilproblem der Photonenraketen in der Schaffung 
von Reflektoren allerhéchsten Reflexionsvermigens besteht. 


b)-Lampe. Wenn das erste Problem einer méglichst verlustlosen Photonenlenkung uns 
zunaichst auf das Gebiet langerer Wellen verweist, so deckt sich diese Tendenz nach Abb.5 
auch mit dem zweiten Hauptproblem der Photonenerzeugung selbst, wo zweifellos im Gebiet 
der sichtbaren und ihnen benachbarten Wellenlangen die umfangreichsten Erfahrungen und 
Méglichkeiten vorliegen. Wegen der viel zu geringen zulassigen Temperaturen und Strahlungs- 
intensitaten fester Kérper kommt in erster Linie Gasstrahlung in Frage. 

Bei der Anregung von Gasmolekiilen entstehen Rotationsspektren im fernen Ultrarot, 
Rotations --Schwingungs-Spektren im nahen Ultrarot und im Sichtbaren, Rotations +-Schwin- 
gungs +Elektronen-Spektren im Sichtbaren und im Ultraviolett. Atomgase liefern Spektren 
der Valenzelektronen und Ionisations-Spektren im Sichtbaren und beiden Nachbargebieten. 
Die Méglichkeiten der Molekil- und Atomanregung von Gasen sind in erster Linie: 


a) thermisch (Flammen, Lichtbogen, Sterne, Entladungslampen, Funkenentladung, Blitz, 
starke VerdichtungsstéBe, DetonationsstéBe, Kernreaktionen usw.), 

b) chemisch (Chemilumineszenz, Flammen, Nachthimmelleuchten, Phosphor, kalte Na- 
trium-Chlor-Reaktion usw.), 

ce) elektrisch (Glimmentladung, elektrische Entladungsréhren, Nordlicht usw.), 

d) optisch (Fluoreszenz, Kometen, Ionosphare usw.). 
Hier interessiert zunachst besonders die thermische Anregungsmiglichkeit von Gasen un- 
mittelbar durch Kernreaktionen. 

Ersetzt man in gliihenden Gasplasmen, wie sie beispielsweise in Hochstrombogen, Queck- 
silberdampflampen, Xenonlampen usw. auftreten, die elektrische Energiezufuhr zum Plasma 
durch eine stationdre Kernreaktion, d.h. den Elektronensto®prozeB® durch den ProzeB der 


XXI. Band 19538. Singer: Zur Theorie der Photonenraketen. 225 


StéBe von Kernreaktionsprodukten, so erhalt man in dieser »Kern-Lampe“‘ ein Gerat zur 
direkten Umsetzung «von Kernenergie in Photonen, wobei die Art der erzeugten Photonen 
durch die Auswahl der Gasfiillung beeinflu&t werden kann. 

Hier eréffnet sich ein zweiter, weiter Forschungsbereich eines Grenzgebietes zwischen 
physikalischer Optik und Atomphysik, mit dem Ziel, den Wirkungsgrad dieser Umsetzung 
in die niitzlichen, mittelwelligen Photonen méglichst hoch zu treiben und kurzwellige Kern- 
- strahlung ebenso wie Gasionisation weitgehend zu unterdriicken, 

Man ist in elektrisch stationiir beheizten Plasmen anscheinend noch nicht iiber stationdre 
Plasmentemperaturen von 50000°K gegangen!, wihrend man sich in diskontinuierlichen 
Hochstrombigen bei den Versuchen der chemischen Ziindung thermischer Kernreaktionen 
sehr kurzzeitig den zehnfachen Werten nihert. 

Jedenfalls erlaubt die Kernheizung eines Plasmas ohne weiteres dauernd die in Detonations- 
stéBen vermuteten? Plasmentemperaturen der Gréenordnung 100000° K aufrechtzuerhalten, 
bei denen praktisch alle Atome des Plasmas zur Strahlung angeregt werden, und die Strahlungs- 
intensitat I bei optisch unendlicher Schichtdicke der Gasmasse auf die GréSenordnung 
10 Watt/cm? steigen kann, entsprechend Lichtdriicken von etwa I/¢ =3,4 Atmosphiren. 

Die in Photonenraketen technisch erforderlichen Lichtdriicke liegen bei den schwach 
beschleunigten Fahrzeugen der interplanetaren Raumfahrt etwa um 0,01 atm., bei den stark 
beschleunigten Fahrzeugen der Luftfahrt und der interstellaren Raumfahrt bis gegen 10 atm., 
also im wesentlichen im Bereich der durch Gasstrahlung wahrscheinlich erzielbaren Werte. 
Dem genannten oberen Lichtdruckbereich entsprechen natiirlich auch die héchsten An- 
forderungen an den Reflektor. 

Die Schichtdicke der strahlenden Gasschicht wird etwa dann als optisch unendlich gelten 
kénnen, wenn sie 


n = 4,6 d2/Q 


mittlere Atomabstande d betragt®, wobei Q der optische Wirkungsquerschnitt eines Atoms ist. 
Zwischen dem Durchmesser D einer kugelférmigen Kernlampe und dem in ihr herrschenden 


Gasdruck p muf demnach die Beziehung bestehen 
pD =4,6kT/Q, 


wenn T die Plasmatemperatur und k die Boltzmann-Konstante sind. Mit beispielsweise 
T=10°°K und Q=10-% cm? folgt pD =6,5-10% atm.cm, d.h. man hat Interesse, eine 
einzige groBe Lampe und nicht viele kleine zur Herstellung des gewiinschten Schubes zu 
beniitzen. 

Derselbe Gesichtspunkt wird wahrscheinlich auch fiir den Vorgang der Kernreaktion selbst 
mafgebend sein, wo ein unzulassiges Abdiffundieren der wirksamen Reaktionspartner durch 
entsprechend groBe Werte von pD verhindert werden muB. 

Z.B. finden Shepherd-Cleaver‘, allerdings auf Grund etwas veralteter Unterlagen, fiir 
einen gasférmigen Wasserstoff-Uranmeiler ein notwendiges pD=1,2-10° atm. cm, also un- 
ausfiihrbar groBe Meilerdurchmesser. Fiir thermische Kernaufbau-Reaktionen scheinen diese 
Werte jedoch ganz erheblich kleiner zu sein, und sich denen der optischen Durchmesser- 
bedingung zu nahern. 

Die bisher ausschlieBlich betrachtete rein thermische Anregung des leuchtenden Gases 
braucht indes nicht die einzige Lésungsméglichkeit zu sein. 

Im Interesse der sehr schwierigen Kithlungsverhaltnisse der Kernlampen liegt es, in Analo- 
gie zur Chemilumineszenz auch eine direkte Ubertragung der Kernreaktionsenergie an Atom- 
anregung ohne den Umweg iiber die iibrigen thermische Freiheitsgrade in Erwagung zu 
ziehen, wobei also z.B. die Translationstemperatur des Gases relativ niedrig bleiben kénnte. 
Damit wiirde natiirlich die Erwarmung der Lampenwande durch Konvektion, aber noch nicht 
jene durch Strahlungsabsor ption infolge mangelhafter Strahlungsdurchlassigkeit vermindert. 


5. Totale Photonenrakete. Die totale Photonenrakete mit einem spezifischen Impuls von 
30500 secton/kg ware auf weitere Zukunft gesehen zweifellos die Endlésung des Strahlantriebs- 
problemes fiir Luftfahrt und Raumfahrt. Derzeit fehlen fiir ihre technische Inangriffnahme 
jedoch noch die wesentlichsten physikalischen Voraussetzungen. 


1 W.H. Westphal, Physikalisches Worterbuch, 5. 467, Berlin 1952. 
2 E, Sanger, Z. Sprengtechnik 1 (1952), S. 151. 

3 FE. Sénger, Z. Brennstoff-Chemie 32 (1951), S. 36. 

4 ZL, R. Shepherd u. A. V. Cleaver, J. Brit. Interpl. Soc. 81 (1949), 3. 36. 
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Die vollstindige Massenzerstrahlung ist in Form der Elektron-Positron-Reaktion vor kaum 
zwanzig Jahren als Laboratoriumsexperiment erstmals gelungen, fiir eine Ubertragbarkeit in 
technischen Mafstab sind noch keinerlei Handhaben bekannt geworden, insbesondere ist die 
offenbar im Kosmos auftretende Zerstrahlung gréBerer Partikel, wie Atome, auch im physi- 
kalischen Experiment noch nicht darstellbar. 

Die Wellenlangen der bei solchen Umwandlungen: entstehenden Photonen liegen im Be- 
reich der Gamma- und kosmischen Strahlung, z.B. fiir die Elektron-Zerstrahlung bei 
2,43-10-1 em, fiir die Proton-Zerstrahlung bei 1,33-10-' cm, fiir die Uranatom-Zerstrahlung 
bei 1,45-10- em. 

Um die nach unseren Betrachtungen zum Reflektor-Problem erwiinschten Photonen im 
Bereich des sichtbaren Lichtes durch Massenzerstrahlung unmittelbar zu erhalten, miiften 
Partikel von etwa einem Millionstel der Elektronenmasse zur Verfiigung stehen, wie sie sich 
vielleicht in den Neutrinos andeuten. Anderseits scheint eine Wellenlangen-Transformation 
im technischen Mafstab, etwa mit Vorrichtungen nach Art der erwogenen Kernlampe, bei 
derart kurzen Primdrwellenlangen nicht sehr erfolgversprechend, wie die auBerordentliche 
Durchdringungskraft der Photonen kosmischer Strahlen zeigt. 

Insgesamt scheinen somit die physikalischen Voraussetzungen fiir eine Inangriffnahme 
technischer Forschungsarbeiten in Richtung der totalen Photonenrakete noch nicht auszu- 
reichen, wahrend diese Voraussetzungen fiir die partiellen Photonenraketen vielleicht schon 
in naher Zukunft bestehen werden, so daB® deren vorbereitende Bearbeitung gerechtfertigt 
erscheint. 


(Eingegangen am 13. November 1952.) 
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Herausgegeben von Reichsbahndirektor i. R. Dr.-Ing. R. KUhnel, Minden, Zweite, 
verbesserte Auflage. Mit 323 Abbildungen. VII, 454 Seiten. 1952. Ganzleinen DM 69,— 


Inhaltsibersicht: Erster Teil: Priifung und Bewertung der Gleitlagerwerkstoffe, Von Reichsbahndirektor i, R. 


Dr.-Ing. R. Ktihnel,- Minden. — Geschichtliche Entwicklung, — Die Lagerpriifverfahren und die Auswertung ~ 


ihrer Ergebnisse fiir die Werkstoffwahl. — Formgebung in Wechselwirkung zum Werkstoff, — Gleitbedingungen 
in Wechselwirkung mit dem Werkstoff. — Bedeutung der einzelnen Werkstoffeigenschaften fiir den Lagerlauf. — 
Zweckmi&Bige Darstellung von Priifungs- und Erfahrungswerten der Lagerwerkstoffe. — Zweiter Teil: Metallische 


\ -Gleitlagerwerkstoffe: Legierungen mit Blei oder Zinn als Hauptbestandteil. Von Professor Dr.-Ing. W. Bungardt, 


Essen. — Gehiirtete Bleilagermetalle. Von Dr.-Ing. R. Weber, Frankfurt a. M. — Legierungen mit Kadmium 
als Hauptbestandteil. Von Dr.-Ing. R. Weber, Frankfurt a. M. — Legierungen mit Aluminium oder Magnesium 
als Hauptbestandteil. Von Professdr Dr.-Ing. W. Bungardt, Essen, — Legierungen mit Silber als Hauptbe- 
standteil. Von Dr. H. Mann, Osnabriick. — Legierungen mit Zink als Hauptbestandteil. Von Reichsbahnrat 
‘Dr. E. Martin, Minden, — Legierungen mit Kupfer als Hauptbestandteil. Von Dr. H. Mann, Osnabriick, — 
Sintermetalle. Von Dr. H. Wiemer, Mehlem (Rhid.). — GuBeisen. Von Dr.-Ing. W. Meboldt, Mannheim, und 
Reichsbahndirektor i. R. Dr.-Ing, R. Kiihnel, Minden. — Dritter Teil: Nichtmetallische Werkstoffe: Holz. Von 
Ober-Ing. j. Arens, Saarbriicken. — Kunststoffe. Von Professor Dr. A.Thum, Ziirich, und Dr.-Ing. C. M. Frh. v. 
Meysenbug, Darmstadt. — Kohle. Von Dr. H. Wiemer, Mehlem (Rhid.). — Sachverzeichnis. 


Die erste Auflage dieses Buches — 1939 — war nach zwei Jahren bereits vergriffem. Die Verhiltnisse der Kriegs- 
und Nachkriegszeit brachten es mit sich, daB erst 1949 mit der Bearbeitung der zweiten Auflage begonnen werden 
konnte. Die Rohstoffnot des vergangenen Jahrzehnts zwang Deutschland, die Forschungen auf dem Gebiet der 
Lagerwerkstoffe besonders vorwirts zu treiben. Die reichen Ergebnisse dieser Untersuchungen bedingten eine 
véllige Neubearbeitung der einzelnen Beitraige und eine Umstellung, Die Abschnitte Grundlagen der Konstruktion 
und Priifung der Lager sind entfallen. Das, was fiir die Lagerwerkstoffe davon noch in Frage kam, ist mit dem 
Abschnitt Priifung und Bewertung der Gleitlagerwerkstoffe vereinigt’ worden, der den ersten Teil des Buches 
bildet. Der dadurch entstandene Raumgewinn kam dem zweiten und dritten Teil zugute. Der zweite Teil enthalt 
eon die metallischen Gleitstoffe und wurde um die Abschnitte Silber und GuSeisen vermebrt. In den dritten 
eil sind die nichtmetallischen Gleitstoffe aufgenommen, Zusitzlich wurde noch Holz und Kohble eingefigt, 


+ 
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Hochwertiges GuBeisen (GrauguB), seine Eigewschaften und die physikalische 
_- Metallurgie seiner Herstellung. Von Dr.-Ing. habil. Eugen Piwowarsky, o. Professor 

der Eisenhiittenkunde, Direktor des Instituts fiir allgemeine Metallkunde und das ge- 

samte GieBereiwesen an der Rheinisch-Westfilischen Technischen Hochschule Aachen. 
Zweite, verbesserte Auflage. Mit 1063 Abbildungen. XII, 1070 Seiten. 1951. 

Ganzleinen DM 135,— 

Aus den Besprechungen:... Die Erkenntnisse tiber die technologischen, mechanischen, elastischen und andere 

Eigenschaften geben bedeutsame Hinweise fiir die Praxis. Auch die Kapitel iiber legiertes GuBeisen, den Ober- 

flichenschutz und die Eigenarten der mechanischen sowie chemischen Widerstandsfihigkeit sind fiir die neu- 


zeitliche Fabrikation wichtig. Von grundlegender Bedeutung sind auch die thermischen Nachbehandlungs-Ver- 
fahren. Mit Interesse wird man auch die Neuerungen auf dem Gebiete der Schmelzapparate verfolgen und aus 


' dem Kapitel iiber die vielseitigen Anwendungsgebiete besten Rat holen, Daf das Werk noch dazu mit hervor- 


ragenden Bildern versehen ist, soll besonders hervorgehoben werden... - 


5 


» Zeitschrift des Vereins Deutscher Ingenicure“ 


Lehrbuch der allgemeinen Metallkunde. von Dr. seorg Masing, o. 6. Pro- 


.fessor an der Universitat Gottingen, Direktor des Instituts fiir Allgemeine Metallkunde 
Géttingen, Unter Mitwirkung von Dr. Kurt Licke, Assistent am Institut fiir Allgemeine 


Metallkunde Gittingen. Mit 495 Abbildungen. XV, 620 Seiten. 1950. 
DM 56,—; Ganzleinen DM 59,60 


Aus den Besprechungen: Das Buch gibt zunichst die allgemeinen chemischen und physikalischen Grundlagen 
des Stoffgebietes an. Hiernach behandelt es die speziellen Grundlagen, wie die Konstitution der Legierungen 
(Heterogene Gleichgewichte vom Einstoffsystem bis zum Vielstoffsystem). Danach folgt eine Besprechung des 
atomistischen Aufbaus des metallischen Kristalles. Nach Ausfiihrungen iiber die Diffusion und die Entstehung 
des kristallinischen Metallkérpers werden die physikalischen Eigenschaften der Metalle, die plastische Verformung, 
die Eigenspannungen, die sog. Erholung und Rekristallisation, die Vorginge im kristallisierten Zustand und Zun- 
dernng und Korrosion theoretisch und praktisch ausfiihrlich behandelt. SchlieBlich folgt in einem Anhang noch 
die Besprechung der wichtigsten Legierungen, wie der des -Hisens, des Kupfers, der Leichtmetalle, des Zinks und 
des Nickels. Eine Zusammenstellung der wichtigsten physikalischen Konstanten der Metalle ist dem Werk voraus 


geschickt... 
Das Werk ist wegen der Fiille des Stoffes und wegen der Beriicksichtigung auch der neuesten Forschung wohl als 
Standardwerk der allgemeinen Metallkunde zu bezeichnen und kann sowohl den Forschern und Lernenden als 


i fa beste empfohlen werden. 
ere eee » Zeitschrift des Vereirs Deutscher Ingenieure“ 


- 


‘des Schnelliufers. Vereinfachte Bestimmung des Spaltverlustes bei Axialrddern. — F. Die Leitvorrichtung: Das 


Strémungsmaschinen. Von Dr.-Ing:, Dr.-Ing.eh. Carl Pfleiderer, Professor an der 


Technischen Hochschule Braunschweig. Mit 200 Abbildungen. XII, 383 Seiten, 1952. , 
Ganzleinen DM 36,— : 


Inhaltsiibersicht: Einleitung: A. Der Stromungsmechanismus im Laufrad; Absolute und relative Bewegung, 
stoBfreier Eintritt. Die Fallhéhe (Férderhéhe) H. Verluste und Wirkungsgrade in der Strémungsmaschine. Be- = 
sondere Formen des Wirkungsgrades bei Gasférderung. Das Moment der Schaufelkrifte und die anf 1 kg bezogene ; a 
Schaufelarbeit H,},. Haufige Anwendungsformen der Hauptgleichung. Der EinfluS der Endlichkeit der Schaufel-. 
zahl. Rechnerische Berticksichtigung des Einflusses der Endlichkeit der Schaufelzahl. Bedeutung des Spaltdruckes, 
Erlauterung von Reaktionsgrad r, Laufzahl u/C und Druckziffer Y. Elementare Erlauterung der Kraftewirkung ‘ 
zwischen Strémung und Schaufel. Wahl des Schaufelwinkels f,, Verwendungsgebiete der verschiedenen Schaufel-. ae 
formen, Die spezifische Drehzahl. Schnelliufigkeit und Wirkungsgrad. Kavitation. Die Saughéhe der Wasserpumpe. ; 
Die Saughéhe der Wasserturbine. Uberschallgrenze bei Verdichtern. Einlaufziffere (bzw. AuslaSwert) 6? und Stré- 
mungswinkel §, an der Saugkante. — B. Entwurt des Laufrades: Festigkeit und Formgebung von Welle und Lauf- 
rad. Der allgemeine Gang der Berechnung der Laufradschaufel, erlautert am Langsamlanfer. Schanfelentwurf bei 
einfacher Schaufelkriimmung. Die doppelt gekrimmte Radialschaufel, Mittelldufer. Schnelliufer, gekennzeichnet 
durch die Schraglage der Druckkante der Schaufel im Meridianschnitt. — C, Austiihrungsbeispiele von Lautradern: 
Kreiselpumpe mit einfach gekriimmter Radialschaufel. Einstufiges Geblise. Schnellaufige Kreiselpumpe. Francis- ~ 
Turbine. Peltonrad. Einstufige Dampfturbine. Einstufige Gasturbine. — D. Die Kennlinien einstutiger Maschinen 
bei vernachlassigbarer Dichtednderung: Allgemeine Regeln. Kenulinien der radialen Kreiselpumpen. Kennlinien der 
Turbinen. Besonderheiten der Schnelliufer. Modellgesetze fiir alle Strsmungsmaschinen. — E, Spaltverlust und . “a 
Achsschub: Der Spaltverlust. Der Achsschub des Langsamlaufers und sein Ausgleich. Spaltverlust und Achsschub 
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Leitrad. Der schaufellose Ringraum (Leitring). Das Spiralgehause. Leitvorrichtung am Saugmund des Laufrades, — a 
G. Besonderheiten.der Axialschaufel: Allgemeines. Axialraider mit enggestellten Schaufeln. Axialrider mit weit aus-_ 
einandergestellten Schaufeln. Die Leitvorrichtung der Axialpumpe. Rechnungsbeispiele fir Axialmaschinen. 
Berechnung der axialen Zwischenstufe einer Dampf- oder Gasturbine. Gleichdruck mit geringer Uberdrack- 
wirkung. Die Schaufelverluste in der axialen Turbinenstufe. — H. Die Mehrstufigkeit bei Turbinen und Verdichtern: 
Die beiden Arten der Abstufung. Einflu8 der Druckabstufung auf den Wirkungsgrad. Ausfiihrungsformen der 
Denes eearens- Ausfihrungsformen und Eigenschaften der Geschwindigkeitsabstufung mit geringer Uberdrack= 
wirkung. Auswirkung der Reibungswarme bei Druckstufung. Die Gitezahl. Mehrgehauseturbinen. Drehzahl und 
Grenzleistung. Wahl des Reaktionsgrades. Gang der Berechnung einer Druckstufenturbine. Radial beaufschlagte — 
Dampfturbinen. Der mehrstufige Verdichter. — 1. Die Kennlinien ein= und mehrstufiger Verdichter und Dampf- — <3 
turbinen: DerVerdichter bei Abweichung vom Berechnungspunkt. Die Dampf- oder-Gasturbine bei Abweichung 
vom Berechnungspunkt. Besondere Bedingungen bei den einzelnen Arten von Strémungsmaschinen, — Lites cm 
verzeichnis. — Sachverzeichnis, : 
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Die Kreiselpumpen fir Flissigkeiten und Gase. Wasserpumpen Ventilatoren 
Turbogeblase / Turbokompressoren. Von Dr.-Ing. C.Pfleiderer, Professor an-der | 
Technischen Hochschule, Braunschweig. Dritte, neubearbeitete Auflage. Mit 353 Text- i 
abbildungen. XI, 518 Seiten. 1949. DM 51,—; Ganzleinen DM 54,60 
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Ventilatoren. Entwurf und Betrieb der Schleuder- und Schraubengeblase. Von Dr.-Ing. 
BrunoEck. Zweite, verbesserte und erweiterte Auflage. Mit 344 Abbildungen, XI, 
304 Seiten. 1952. Ganzleinen DM 36,— 


Inhaltstibersicht: A. Radialgeblase: Elementare Stromfadentheorie. Genauere rechnerische Behandlung der 

Schaufelstrémung. Einflu8 der endlichen Schaufelzahl. Gestaltung der Schaufelenden. Ahnlichkeitsbeziechungen. 
Verluste. Giinstigste Gestaltung des Laufrades. Betriebseigenschaften von _ Radialgeblasen, Leitvorrichtungen, 
Ermittlung der Hauptdimensionen eines Geblises. Sonderausfiihrungen und Sonderprobleme. — B. Theorie und 
Berechnung von Axialgeblasen: Berechnung nach der Gittertheorie. Ermittlung der Hauptdimensionen eines 
Axialgeblises. Betriebseigenschaften von Axialgeblisen, Meridianbeschleunigte Axialgeblise. — ©. Konstruk- 


tive Gestaltung der Geblase, Sonderausfiihrungen. — D. Festigkeitsberechnungen. — E. Experimentette Erprohung =. 
von Ventilatoren. — Namen- und Sachverzeichnis. 


Technische Stromungslehre. von D:-Ing. Bruno Eck. Vierte Auflag. ‘| 


In Vorbereitung, 
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